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The  results  of  the  two  models  are  combined  to  obtain  the  dynamic 
motions  of  the  gun  about  the  initial  aim^position.  Finally,  the  statistical 
mean  and  the  standard  deviation  errors  in  gun  pointing  are  derived  from  the 
dynamic  gun  vibration  patterns.  - • 

These  analytical  techniques  have  been  used  to  predict  the  gun-pointing 
errors  for  several  types  of  weapons  and  to  investigate  the  dynamic  effects  of 
certain  proposed  modifications  of  the  servocontrols.  Some  specific  weapons  that 
have  been  investigated  include  the  20mm  XM197  gun,  the  Bushmaster  25mm  gun,  and 
the  AMCAWS  30mm  gun.  Both  the  standard  spring  adapter  mounts  and  the  constant 
recoil  mounts  were  considered  in  the  analysis;  the  weapons  were  mounted  on  a 
modified  version  of  the  XM97  turret  that  was  attached  to  the  Bell  AH-1  heli- 
copter airframe.  The  analysis  of  each  weapon  system  is  performed  for  several 
initial  orientations  of  the  turret  to  determine  the  sensitivity  of  the  error  to 
the  gun-pointing  direction. 

The  major  findings  from  these  investigations  are  that  the  constant  recoil 
device  improves  the  gun-pointing  accuracy;  that  a simple  rate-gyro-feedback  in 
the  turret  servocontrols  does  not  significantly  improve  the  dynamic  response  of 
the  weapon;  and,  finally,  that  the  dynamic  response  of  the  weapon  is  rather  sen- 
sitive to  the  initial  turret  orientation  and  to  the  weapon  firing  rate. 
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I.  INTRODUCTION 


One  of  the  factors  that  contributes  to  the  gun-pointing  errors  of  the 
helicopter  automatic  cannon  system  is  the  dynamic  vibration  of  the  system 
components.  The  vibrating  components  of  the  automatic  cannon  system  con- 
sist of  the  helicopter  structure,  the  gun  turret,  and  the  gun  itself.  When 
this  system  is  excited  by  a dynamic  force,  the  resulting  motion  will  alter 
the  gun  orientation,  which  in  turn  will  affect  the  accuracy  of  the  weapon 
system. 

To  analyze  and  predict  the  gun  errors  due  to  this  effect  necessitated 
extensive  analytical  and  computer  studies.  The  purpose  of  the  present  re- 
port is  to  describe  these  studies  and  present  results  that  have  been  obtained 
in  the  course  of  this  investigation. 

For  the  dynamic  analysis,  the  helicopter- weapon  system  was  divided  into 
two  parts.  The  first  part  of  the  system,  called  the  helicopter  model,  is  a 
detailed  representation  of  the  helicopter  mass  and  structural  components, 
with  the  turret  and  the  gun  represented  by  a concentrated  mass.  The  second 
part  of  the  system,  called  the  gun-turret  model,  la  a dynamic  model  used 
to  describe  the  moving  components  of  the  turret  and  the  gun.  The  dynamic 
motion  of  the  helicopter  is  analyzed  by  use  of  a finite-element  model  of 
the  Bell  AH-1  structural  airfr-oe.  A prescribed  forcing  function  is  applied 
to  induce  the  structural  vibrations.  Since  this  model  contains  a relatively 
large  number  of  degrees  of  freedom  (approximately  1800),  a considerable  re- 
duction in  the  order  of  the  system  is  firs:  achieved  by  a matrix  condensation 
procedure.  Then  the  characteristic  frequencies  and  the  mode  shapes  of  the 
structure  are  determined,  and  the  governing  dynamic  equations  of  motion  are 
reformulated  in  terms  of  a set  of  uncoupled  modal  equations.  The  modal 
equations  are  solved  numerically  to  obtain  the  translational  and  the  angular 
accelerations  of  the  concentrated  mass,  representing  the  turret  and  the 
gun.  These  accelerations  are  then  used  as  an  input  into  the  gun-turret  dy- 
namic model,  which  is  a two-degree-of-f reedom  model  representing  the  azimuth 
and  the  elevation  motions  ot  one  turret. 

The  response  of  the  turret  is  affected  by  the  inertial  properties  of  the 
turret  and  the  gun,  the  characteristics  of  the  servocontrols , and  the  stiff- 
ness of  the  drive  systems  usee,  for  the  azimuth  and  elevation  controls. 


1 


The  dynamic  equations  for  this  system  are  solved  by  numerical  methods,  and 
the  solutions  for  the  azimuth  and  elevation  dynamic  motions  are  presented 
in  a form  which  is  useful  for  the  statistical  analysis  of  gun  errors.  In 
particular,  the  mean  and  the  standard  deviation  errors  of  the  azimuth  and 
the  elevation  motions  are  determined  from  the  vibrations  about  the  initial 

aim-position.  These  statistical  parameters  are  evaluated  for  the  time  

interval  during  which  the  gun  is  firing. 

Two.  different  finite-element  structural  models  are  used  to  represent 
the  AH-U  and  AH-1G  helicopter  structures.  These  two  models  differ  primarily 
in  the  way  that  the  structure  of  the  helicopter  is  modeled.  The  AH-U  model, 
described  in  Reference  1,  is  an  elastic  axis  NASTRAN  representation  of  the 
structural  airframe.  The  AH-1G  NASTRAN  model,  described  in  Reference  2, 
contains  a much  more  detailed  representation  of  the  fuselage  and  wingstub 
structures.  Both  models  have  been  validated  by  comparisons  of  model  pre- 
dictions with  results  from  static  deflection  tests  and  dynamic  shake  tests. 
Structural  plots  of  the  AH-U  and  AH-1G  models  are  shown  ir.  Figures  1 and 
2,  respectively.  

A number  of  different  forcing  functions  can  initiate  the  dynamic 
vibration  of  the  total  system.  The  main  unsteady  forces  are  the  main  rotor 
loads  and  the  recoil  force  of  the  gun.  Experimental  data  from  a recent 
series  of  flight  tests  of  the  AH-lG  helicopter  show  that  the  recoil  forces 
are,  by  far,  the  most  predominant  in  producing  airframe  and  turret  vibra- 
tions. Consequently,  numerical  results  are  obtained  and  presented  for  the 
recoil  force  excitation  only.  Weapons  of  different  caliber  have  been  ex- 
amined, including  20,  25,  and  30  mm  automatic  cannons.  Different  recoil 
forces  for  each  of  these  weapons  were  obtained  from  experimental  tests  and 
from  computer  simulations  of  the  assembled  weapon-and-mount  mechanisms. 


^ronkhite,  J.D.  and  Wilson,  W.F. , "Dynamic  Analysis  of  Two-Per-Rev 
Vibrations  in  the  Model  AH-1J  Helicopter",  Report  No.  299-100-021,  Bell 
Helicopter  Company,  Fort  Worth,  Texas  (Feb  1972) 

2Cronkhite,  J.D.,  and  3erry,  V.7.  and  Brunken,  J.E.,  "A  NASTRAN  Vibra- 
tion Model  of  the  AH-lG  Helicopcer  Airframe",  Volumes  I and  II,  Bell 
Helicopter  Company,  AD  R-TR-74-045.  Fort  Worth,  Texas  (April  1974) 
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II.  ANALYTICAL  PROCEDURE 


A.  Factors  Affecting  the  Gun-Pointing  Error  Analysis 
Dynamic  gun-pointing  errors,  in  a turreted  attack-helicopter,  are 
caused  by  the  dynamic  effects  of  weapon  recoil.  Essentially,  two  primary 
errors  are  attributed  to  the  weapon  recoil  loads  - namely,  the  helicopter 
airframe  vibrations  and  the  turret  servocontrol  oscillations.  The  airframe 
vibrations  that  affect  the  gun-pointing  accuracy  are  the  three  angular 
vibration  components  at  the  turret  attachment  platform.  Since  the  turret 
consists  of  a stiff  linkage  of  rigid  body  segments,  the  oscillations  of  the 
turret  platform  are  a direct  source  of  gun-pointing  error.  In  addition, 
the  relative  oscillations  between  the  interconnected  rigid  body  segments 
of  the  turret  (i.e.  motions  of  the  azimuth  and  elevation  drive  systems) 
constitute  a second  source  of  gun-pointing  error.  The  turret  servocontrol 
system  responds  dynamically  to  the  applied  recoil  torques  and  also  to  the 
accelerations  of  the  turret  attachment  platform. 

A two-step  analytical  procedure  has  been  developed  to  determine  the 
cumulative  gun-pointing  error  due  to  weapon  recoil  loads.  First,  the  air- 
frame flexural  vibrations  are  determined  with  a structural  vibration  model. 
The  location  of  the  gun-barrel  axis  is  defined  in  the  model  by  azimuth  and 
elevation  angles,  and  the  predicted  turret  platform  vibrations  are  resolved 
in  terms  of  the  azimuth  and  elevation  gun-line  components.  Furthermore,  the 
accelerations  of  the  turret  attachment  platform  are  determined  by  the 
structural  vibration  model.  The  dynamics  of  the  turret  servocontrol  system 
are  then  predicted  with  a model  that  utilizes  the  platform  accelerations 
and  recoil  data  generated  by  the  structural  vibration  model.  In  the  turret 
model,  the  turret  servocontrol  oscillations  are  calculated  and  the  resulting 
gun-line  turret  vibrations  are  combined  with  the  platform  vibrations. 

This  procedure  leads  to  predictions  of  the  time  variations  in  the  gun- 
line  azimuth  and  elevation  angles  about  the  initial  aim-position.  Statis- 
tical error  parameters  are  then  calculated  for  the  burst  interval.  In 
particular,  the  means  and  the  standard  deviations  of  the  dynamic  azimuth 
and  elevation  angle  errors  are  determine',  and  these  parameters  are  used  to 
represent  the  gun-poir.ting  error  components  statistically. 

Since  the  magnitude  of  each  statistical  error  parameter  is  dependent 
on  the  initial  aim-position,  a matrix  of  statistical  error  parameters  is 
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calculated  for  several  discrete  turret  orientations.  In  this  manner,  the 
sensitivity  of  the  gun-pointing  error  to  the  initial  aim-position  is  de- 
termined over  the  normal  operating  range  of  the  turret.  The  turret  con- 
sidered in  this  analysis  operates  approximately  between  0°  and  + 110°  in 
azimuth,  and  between  0°  and  -50°  in  elevation.  However,  in  practical  ap- 
plications, the  azimuth  angle  position  normally  does  not  exceed  + 70°. 

To  show  how  the  present  investigation  fits  into  the  overall  error 
analysis  of  a helicopter  mounted  automatic  cannon,  a schematic  diagram  of 
the  Helicopter  Weapon  System  Model  is  presented  in  Figure  3.  This  figure 
shows  the  flow  chart  representing  the  order  of  the  various  calculations, 
which  are  performed  to  determine  the  hit  probability  for  a specified  set 
of  input  parameters. 

With  reference  to  Figure  3,  note  that,  with  this  model,  four  different 
types  of  input  parameters  are  used  for  the  evaluation  of  the  final  hit 
probability.  The  four  inputs,  shown  on  the  left  side  of  the  flow  chart, 
consist  of  the  line  of  sight  and  aircraft  sensor  data,  the  weapon  recoil 
force,  wind  and  ballistic  data,  and  the  target  data.  The  present  investi- 
gation is  designed  to  develop  a portion  of  this  system  model  that  comprises 
the  Turret  Servocontrol  Model,  the  Turret  Dynamic  Model,  and  the  Helicopter 
Vibration  Model.  This  portion  of  the  system  model  responds  to  two  different 
input  functions:  the  first  is  the  gun  aim  angle  commands  passed  from  the 
fire  control  computer  model  and  the  second  is  the  weapon  recoil  force. 
During  the  target  tracking  process,  the  command  aim-angles  produced  by  the 
fire  control  computer  are  relatively  slowly  varying  input  signals  into  the 
turret  control  system  and,  therefore,  do  not,  in  general,  significantly 
affect  the  dynamic  vibration  of  the  system.  Consequently,  the  main  exci- 
tation for  the  dynamic  motion  is  the  recoil  force  of  the  weapon. 

B.  Finite-Element  Structural  Models 

The  helicopter  vibration  analysis,  to  be  described  in  the  following 
section  of  this  report,  is  used  to  predict  the  dynamic  structural  response 
to  a general,  time-dependent  recoil  force.  In  this  analysis,  the  natural 
frequencies,  the  generalized  masses,  and  the  f ree-vibration  modes  shapes 
of  the  helicopter  structure  are  used  in  the  governing  dvnamic  equations  to 
represent  the  dynamic  characteristics  of  the  airframe.  Before  a 
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forced-vibration  problem  can  be-  solved,  the  free-vibration  properties 
of  the  structure  have  to  be  obtained.  In  the  present  investigation,  these 
properties  were  supplied  from  the  finite-element  NASTRAN  computer  program. 
This  program  is  a general-purpose  structural  analysis  program  developed  by 
NASA  and  by  several  aerospace  companies.  One  of  the  capabilities  of  the 
NASTRAN  program  is  to  solve  for  the  natural  frequencies  and  mode  shapes 
of  the  multi-degree-of-freedom  systems. 

In  the  modeling  of  the  AH-1  helicopter  structure  in  terms  of  finite- 
elements,  'two  different  structural  models  were  used.  These  two  models  cor- 
respond to  the  AH-1J  and  the  AH-1G  helicopter  structures.  Structural  plots 
of  the  AH-lJ  and  the  AH-1G  models  are  shown  in  Figures  1 and  2,  respectively, 
iu  spite  of  the  minor  structural  differences  in  the  AH-lJ  and  in  the  AH-1G 
helicopter  configurations,  the  two  models  differ  substantially  in 
the  detail  of  the  structural  representations.  In  the  case  of  the  AH-lJ 
model,  the  helicopter  structure  was  represented  by  the  so-called  "elastic 
axis"  configuration.  in  this  approach,  the  structure  is  modeled  entirely 
by  beam  finite  elements  that  can  carry  axial  tension,  bending,  and  torsion. 
The  mass  of  the  system  was  lumped  at  the  nodes.  This  is  a relatively  simple 
representation,  as  is  evident  from  the  structural  plot  of  the  AH-lJ  model 
in  Figure  1.  For  example,  the  main  fuselage  of  the  helicopter  ^.s  modeled 
by  a row  of  elements  in  series.  In  the  case  of  the  AH-1G  structure,  a more 
detailed  finite-element  model  representation  was  developed.  In  this  case, 
the  total  structure  was  modeled  by  many  different  types  of  elements  in- 
cluding beams,  rods,  plates,  and  shear  panels.  Consequently,  this  led  to  a 
larger  degree-of-freedom  system.  As  in  the  first  model,  the  mass  of  the 
helicopter  was  lumped  at  the  codes  of  the  finite  elements.  Both  in  the 
AH-lJ  and  in  the  AH-1G  models,  one  concentrated  mass  was  chosen  to  repre- 
sent the  mass  of  the  gun-turret  system.  This  mass  was  placed  as  close  as 
possible  to  the  center  of  gravity  of  the  turret  and  weapon. 

C.  Helicopter  Vibration  Model  Description 

In  the  finite-element  NASTRAN  formulation  of  the  structural  model,  the 
helicopter  structure  is  idealized  by  a set  of  discrete  structural  elements 
connected  at  suitable  nodal  or  grid  points.  For  this  type  of  model,  the 
dynamic  equations  for  the  forced  vibration  of  the  helicopter  can  be  written 
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(i) 


in  the  matrix  form 

[m]  {»]  + Ik]  {x)  - {F} 

where 

[m]  is  the  mass  matrix, 

[k]  is  the  stiffness  matrix, 

(F)  is  the  force  matrix, 

and  (x)  is  the  vector  of  the  grid  point  displacements. 

Note  that  equation  (1)  does  not  contain  any  damping  terms.  With  re- 
spect to  the  present  problem,  the  incorporation  of  terms  to  represent  struc- 
tural damping  is  desirible;  but,  because  of  its  frictional  nature,  damping 
can  only  be  represented  empirically.  Such  a correction  can  be  done  most 
conveniently  after  the  modal  formulation  has  been  introduced  and  the  equa- 
tions (1)  have  been  decoupled. 

To  transform  the  matrix  equation  (1)  to  the  modal  formulation  necessi- 
tated the  examination  of  the  free-vibration  characteristics  of  the  homo- 
geneous  problem  corresponding  to  equation  (1) , 

[m]  {X}  + [k]  (x)  * 0.  (2) 

The  eigenvalues  for  equation  (2)  will  be  denoted  by  w2,  and  these  are 
obtained  from  a characteristic  equation 

| - w2  [m]  + [k]  | = 0 (3) 

where  the  vertical  bars  indicate  a determinant.  This  procedure  results 
in  the  determination  of  the  characteristic  frequencies,  where  i varies 
over  the  number  of  degrees  of  freedom  of  the  system.  Once  the  frequencies 
are  known,  the  corresponding  eigenvectors  or  mode-shape  vectors  are 

determined  by  a standard  procedure.  Both  of  these  steps  were  carried  out 
by  the  NASTRAN  computer  program.  The  modal  transformation  matrix  is  defined 
by  the  relation 

(4) 


[$]  “ [♦j.t  $2*  • • • • 

where  N is  the  number  of  degrees  of  freedom.  In  the  definition  (4) , the 
modal  matrix  consists  of  column  vectors  of  the  mode  shapes. 

Consider  now  a modal  coordinate  transformation 

{x}  - [*]  {q}.  (5) 
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By  the  substitution  of  equation  (5)  in  equation  (1)  and  the  multiplication 

T 


of  both  sides  by  the  transposed  matrix  [$]  , the  dynamic  equations  can  be 
written  in  the  new  form 


[*]T  In]  I*]  £*}  + U]T  Ik]  M {q}  - WT  {F}. 


(6) 


The  normalization  of  the  modal  matrix  U]and  the  use  of  the  orthogonal 
properties  of  the  natural  modes  lead  to  the  following  simplifications: 


M1  Im]  !♦]  - HI 


(7) 

and  M*  Ik]  [+]  - [V]  (8) 

where  [I]  is  a unit  diagonal  matrix  and  [''u?]  is  a diagonal  matrix  contain- 


ing the  squares  of  the  natural  frequencies.  From  equations  (6),  (7),  and 
(8)  a set  of  uncoupled,  second-order  differential  equations  is  obtained 


in  the  variables  q^. 


^ + “i  qi  " Q 


(i-  1,  2 , 3,  . 


N) 


(9) 


where  no  summation  on  i is  implied  and  is  a generalized  force  defined  by 


Qi  - IF}. 


The  consequence  of  equations  (9)  is  that  the  helicopter  vibration  problem 
has  been  reduced  to  a solution  of  a large  number  of  uncoupled  differential 
equations. 

Normally,  however,  an  accurate  solution  to  the  forced-vibration  problem 


can  be  obtained  by  the  solution  to  only  those  equations  for  which  is 


below  some  cutoff  frequency.  For  the  helicopter  vibration  problenv  only  the 
first  30  of  equations  (9)  are  solved  and  the  higher  frequency  equations  are 
neglected.  In  general,  the  number  of  equations  that  need  to  be  retained  to 
obtain  an  accurate  solution  is  dependent  upon  the  nature  of  the  excitation 
force  vector  {F}. 

Before  the  equations  (9)  are  solved,  they  are  modified  to  Include  struc- 
tural damping  effects.  Any  correction  for  this  type  of  damping  has  to 
satisfy  the  experimental  observations  that  the  damping  force  is  in  phase 
with  velocity  and  is  Independent  of  frequency.  When  the  structure  is  subject 
to  sinusoidal  vibration,  these  facts  can  be  represented  mathematically  by 
the  introduction  of  the  complex  stiffness  and  the  replacement  of  the  frequency 
term  w*  in  equations  (9)  by  the  complex  quantity  (1  + jy),  where  j rep- 
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resents  the  pure  imaginary  number  /^T  and  y is  defined  as  the  damping  coeffi- 
cient. The  complex  representation  for  the  structural  damping  is  not~conven- 
ient  when  it  is  applied  to  general  time-dependent  problems.  Consequently, 
in  this  case,  a more  general  correction  is  applied  that  include  the  first 
time-derivative . In  this  approach,  the  term  in  equations  (9)  is  replaced 


by  ui£(l  + y/u^  );  this  substitution  leads  to  the  modified  dynamic  equations 


^ q±  + wj  q±  = 


a = 1.  2,  3, 


.,  N). 


(10) 


Note  that  the  behavior  of  the  damping  term  in  these  equations  reduces 
to  the  desired  form  when  the  system  is  excited  with  a sinusoidal  force 


with  frequency  to^.  In  general,  the  value  of  the  damping  coefficient  y has 


to  be  obtained  experimentally;  however,  for  the  present  analysis,  a value 
of  0.02  was  used  because  this  value  has  been  found  by  previous  investigators 
to  be  quite  representative  of  airframe  construction.  _ _ 

Before  the  solutions  to  equations  (10)  are  obtained  for  general  forcing 


functions,  these  equations  are  first  solved  for  impulsive  loading.  Dropping 
the  subscript  i for  convenience  results  in  the  following  equation  for  the  im- 
pulsive load  problem: 

q + y o>  q + ui2q  - <5  (t)  (n) 

where  6 (t)  is  a unit  delta  function.  The  solution  to  equation  (n)  is  de- 
noted by  q^  and  can  be  obtained  by  the  expression  of  equation  (11)  as  a homo- 
geneous equation  with  a nonzero  initial  velocity  condition.  This  solution 
is  given  by 


<6 

where  a 
and  8 


1 at  . 
g e sin  8t 


yu/2 


(y/2)4 


(12) 

(13a) 

(13b) 


The  solution  for  the  general  forcing  function  Q(t)  can  be  obtained  by  use  of 
equation  (12)  and  the  superposition  integral.  The  general  solution  to 
equation  (10)  is  then  given  by 

q(t)  - f Q(x)  q^  (t  - r)  d t. 

' n 


(14) 


By  the  substitution  of  equation  (12)  in  (14),  the  solution  for  q(t)  is  then 


f 


■ 


3 


P 

1 


A 


ii 


given  in  the  following  form: 

q(t)  - | eot  jsin  8t  jQ  Q(t)  e~GT  cos  St  d t 

Q(t)  e °T  sin  0t  d t > . (15) 

o ) 

In  the  present  analysis,  equation  (15)  was  used  for  the  solution  to  each 
modal  equation,  and  the  integration  with  respect  to  time  was  carried  out 
numerically.  Altogether,  30  modal  equations  were  used  in  the  analysis  for 
both  the  AH-U  and  the  AH- 10  helicopter  models. 

After  the  set  of  modal  equations  is  solved  to  obtain  the  modal  coordi- 
nate vector  (q),  the  grid-point  displacement  vector  {x}  is  determined  from 
equation  (5).  The  results  of  this  analysis  are  used  as  the  input  to  the 
Gun-Turret  Dynamic  Model,  which  is  described  in  the  next  section  of  this 
report.  This  input  comprises  the  angular  displacements,  velocities,  and 
accelerations,  and  the  linear  accelerations  of  the  grid  point  corresponding 
to  the  turret  attachment  platform.  Therefore,  only  the  six  displacement 
terms  that  correspond  to  the  turret  attachment  grid  point  in  the  vector  {x} 
need  to  be  evaluated.  The  velocities  and  accelerations  of  the  turret  attach- 
ment grid  point  are  then  calculated  by  differentiation. 

The  FORTRAN  IV  source  program  for  the  Helicopter  Vibration  hodel  is 
listed  in  the  Appendix  A.  The  output  data  that  is  required,  in  the  analysis, 
by  the  Gun-Turret  Dynamic  Model  is  written  on  a sequential  data  set  (file 
12).  Two  temporary  I/O  files  (files  13  and  14)  are  also  used  in  the  program, 
for  the  storage  of  the  intermediate  results.  These  temporary  files  are 
provided  primarilv  to  reduce  the  central  core  storage  requirements  and  to 
Increase  the  execut  on  efficiency. 

D.  Gun-Turret  Dynamic  Model 

In  this  section,  the  dynamic  model  for  the  gun-turret  configuration, 
including  the  effects  of  the  servocontrols  in  the  turret,  is  described. 

The  forcing  function  for  the  dynamic  motion  is  provided  by  the  gun  recoil 
force  and  by  the  accelerations  of  the  helicopter  structure  to  which  the 
turret  is  attached.  The  dynamic  motion  of  the  turret  can  be  described  by 
two  rotational  equations  of  motion:  one  tor  the  azimuth  motion,  and  the 

other  for  the  elevation  motion.  The  three  main  components  of  the  turret 
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system,  shewn  in  Figure  4,  consist  of  the  stationary  ring,  the  rotating 
ring,  and  the  gun.  The  stationary  ring  is  attached  to  the  helicopter  plat- 
form and  will  therefore  move  with  the  vibration  of  the  helicopter  structure. 
The  azimuth  motion  will  involve  the  rotation  of  the  rotating  ring,  the  saddle 
mechanism,  and  the  gun  about  the  turret  azimuth  axis.  The  elevation  motion 
) will  involve  the  rotation  of  the  gun  and  saddle  mechanism  about  the  eleva- 

tion trunnion  axis. 


Axis 

Figure  4.  Schematic  Diagram  of  the  Turret  and  Gun  Configuration 

In  the  derivation  of  the  dynamic  equations,  various  coordinate  systems 
are  used.  The  three  different  coordinates  systems  x^,  and  z^  (i  = 1,  2, 
3)  used  in  the  analysis  are  shown  in  Figure  5.  The  coordinates  are  the 
inertial  coordinates,  and  the  coordinates  x^  are  moving  coordinates  that 
rotate  with  the  helicopter  vibration  and,  therefore,  differ  from  by  small 
angular  deflections.  The  coordinates  z^  are  oriented  with  respect  to  the 
coordinates  x±  by  the  initial  azimuth  and  elevation  angles  that  are  denoted 
by  and  a,  respectively.  The  azimuth  dynamic  equation  refers  to  the  ro- 
tation about  the  axis  X2»  and  the  elevation  dynamic  equation  refers  to  the 


rotation  about  the  axis  z^.  The  origin  of  the  coordinates  is  assumed 
to  be  situated  at  the  concentrated  mass  corresponding  to  the  gun-turret 
configuration  in  the  NASTRAN  finite-element  model. 


Figure  5.  Coordinate  Systems  Used  in  the  Analysis 


To  obtain  the  rotational  equations  of  motion  about  the  axes  and  z 
involves  the  use  of  the  general  vector  dynamic  equations  for  a moving  co- 
ordinate system.  The  moment  vector  equation  is  expressed  as 


M = m p xa  + H 
c p 


M is  the  vector  moment. 


Pc  is  the  vector  position  of  the  center  of  gravity  of  the  rotating 
mass  relative  to  the  £ origin, 

a is  the  linear  acceleration  of  the  x^  origin, 

m is  the  mass  of  the  rotating  parts, 


H is  the  portion  of  the  rate  of  change  of  angular  momentum  that  is 
dependent  upon  the  moments  of  inertia  and  the  angular  accelerations 
of  the  rotating  parts. 

The  relative  positions  of  the  various  pivot  points  and  centers  of  mass  are 
defined  in  Figure  4 by  the  parameters  A^,  k^,  ky  B^,  and  Bj.  The  center 
of  gravity  CG^  is  denoted  by  the  length  parameters  and  B^  defined  with 
reference  to  zero  azimuth  and  elevation  angles.  The  center  of  gravity  CG^ 
applies  to  the  total  mass  involved  in  the  azimuth  vibration,  which  includes 
the  rotating  ring  and  the  saddle-gun  combination.  The  center  of  gravity 
CG2  applies  to  the  saddle-gun  combination  only,  and  it  is  defined  by  the 
lengths  ky  ky  and  B^.  The  parameter  A^  defines  the  pivot  position,  and 
A^  and  B^  define  the  postion  of  CG2. 

Since  the  dynamic  equations  for  the  azimuth  and  the  elevation  motions 
are  similar  in  form,  a general  set  of  equations  can  be  derived  that  are 
valid  for  either  of  these  motions.  Consider  the  schematic  representation 
in  Figure  6 in  which  the  three  components  involved  in  the  vibration  are 
shown.  This  schematic  applies  both  to  the  azimuth  and  to  the  elevation 
motions.  The  three  components  in  Figure  6 are  the  drive  motor,  the  gear 
transmission  system,  and  the  turret  mass.  For  convenience,  the  moments  of 
inertia  of  the  motor  and  the  turret  are  referred  to  the  turret  coordinates. 
The  moment  of  inertia  of  the  gear  box  and  the  drive  can  be  included  in  the 
motor  and  the  turret  inertias.  The  torques  transmitted  through  the  total 
system  are  also  shown  in  Figure  6.  The  torque  T^  is  the  torque  applied 
electrically  to  the  motor,  and  Tfc  is  the  torque  transmitted  to  the  turret. 
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Figure  6.  Schematic  Diagram  of  the  Moving  Parte  in  the  Turret-Gun 
Dynamic  Analysis 


The  effect  of  friction  in  the  system  is  represented  by  the  torque  T^  act- 


ing on  the  turret.  The  torque  Tr  represents  the  contribution  of  the  recoil 


force.  The  rates  of  change  of  the  angular  momentum  for  the  motor  and  for 
thr  turret  are  denoted  by  and  P , respectively.  With  the  use  of  the 
torques  shown  in  Figure  6,  the  dynamic  equations  of  equilibrium  are  written 


as  follows : 


H - T - T 
m m t 


H “ T + T + T, 
t r t f 


The  transmitted  torque  Tfc  is  related  to  the  elastic  deformation  of 


the  system.  To  express  equations  (17)  and  (18)  in  terms  of  angular  dis- 
placements, let  and  Aa  define  dynamic  rotations  in  the  azimuth  and 


elevation  mechanisms,  respectively;  and  furthermore  let 

a*  t ♦„  + K + ♦„»  <«> 

o ffl  e 

and 

An  ^ a + a + a , (20) 

o m e’ 

where  and  aQ  are  the  dynamic  components  due  to  helicopter  vibration, 
and  are  the  dynamic  motions  for  rigid  gear  system  drives,  and 
and  cie  are  the  contributions  from  the  elastic  effects  in  the  system.  For 
convenience,  let  ^ and  define  the  relative  displacements  given  by  the 
expressions 

♦r  " + ♦*»  (21) 

r m e 

and 

a - a + a . (22) 

r m e 

The  angles  $ and  are  relative  angles  between  various  parts  of  the 

turret  and  are  detected  by  the  angular  position  resolvers  in  the  servocontrol 
system.  The  motor  torque  in  equations  (17)  and  (18)  can  be  expressed  in 
terms  of  the  elastic  displacements.  For  the  azimuth  motion: 


and  for  the  elevation  motion: 


where  k represents  the  elastic  stiffness  of  the  gear  box  and  the  transmission. 

Consider  now  the  torque  T^  produced  by  the  friction  in  the  system. 

Since  this  term  represents  a combination  of  all  frictional  forces  in  the 
system,  an  exact  analysis  of  this  term  is  extremely  difficult,  and  an  em- 
pirical relation  has  to  be  used.  In  this  analysis,  the  assumption  is  that 
this  torque  acts  to  oppose  the  relative  motion  of  the  various  parts  of  the 
turret;  but  since  it  is  a frictional  force,  it  is  independent  of  the  magni- 
tude of  the  angular  velocities.  Consequently,  for  the  azimuth  motion,  thi3 
torque  is  represented  by: 


Tf  - -kg(<frr)/|$r| 
and  for  the  elevation  motion 

Tf  - -kg<*r)/|Arl 
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where  g Is  the  damping  coefficient. 

Consider  now  the  rates  of  change  of  angular  momentum  and  in  equa- 
tions (17)  and  (18).  Before  writing  the  expression  for  Ht>  let  a^,  and 
a^  denote  the  linear  acceleration  components  of  the  turret  attachment  plat- 
form in  the  axis  system  and  let  9^,  and  6^  denote  the  angular  acceler- 
ation components.  These  are  the  hull  accelerations  that  are  calculated  by 
the  Helicopter  Vibration  Model.  With  the  use  of  the  azimuth  and  elevation 
angles  ^ and  a,  and  the  geometrical  parameters  shown  in  Figure  4,  the  com- 
ponent of  the  rate  of  change  of  angular  momentum  for  the  turret  about  the 
axis  of  rotation  is  given  by 

&t  " -m^A^(-a^  sin  <j>  + Sj  cos  <j>)  + I (27) 

where  m^  is  the  turret  mass  involved  in  the  azimuth  motion  and  It  is  the 
turret  moment  of  inertia  about  the  turret  azimuth  axis.  Consider  now  the 
angular  acceleration  term  in  equation  (27).  Differentiating  equations  (19) 
and  (21)  twice  result  in  the  expression 


■ <f>  + $ 

To  Yr 

But  $ is  related  to  the  hull  vibration  in  the  azimuth  direction 
o 


(28) 

by 


0. 


(29) 


o 3 

For  the  motor,  the  predominant  contribution  to  the  angular  momentum  rate  is 
from  the  relative  motor  rotations  and,  therefore,  to  a good  approximation 

(30) 


H 


m 


I <f> 

m m 


where  I is  the  moment  of  inertia  of  the  motor  referred  to  the  turret  coordi- 
m 

nates. 

For  the  case  of  the  elevdtion  vibration,  the  rotational  equilibrium  equa- 
tions are  referred  to  the  saddle  pivot  axis.  The  rate  of  change  ot  the  angular 
momentum  can  be  shown  to  be 

Ht  “ m2  1(^2  + A3  coe  °0  a3  + ( cos  <t>  +a2  sin  (_B2  + A3  sin  °01 

+ It  A a (3D 

where  m^  is  the  mass  of  the  saddle-gun  combination  and  I is  the  moment  of 
inertia  about  the  saddle  pivot  axis.  For  the  elevation  motor,  all  terms 
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except  the  motor  rotational  effects  car.  be  neglected  and  therefore, 


H - 1 a 
m mm 


(32) 


Equations  (20)  and  (22)  can  be  differentiated  to  get 


Ad  = a +o 
o r 


(33) 


and,  as  is  shown  in  the  next  section,  the  acceleration  a is  related  to 

o 

the  hull  accelerations  by 
••  ••  •• 

°*o  * “e1  sin  d>  + e2  cos  . (34) 

Consider  now  the  azimuth  vibration.  By  use  of  equations  (17)  and  (18) 
and  substitution  from  equations  (23),  (25),  (27),  (28),  (29),  and  (30),  the 
equations  of  equilibrium  can  be  written  in  the  form: 


I ♦ - k<(> 
mm  e 


m 


+ M>e  + k84>r/|$rl 


m^A^(-a^  sin  <J>  + a^  cos  <f>)  - Tr 


(35) 


(36) 


Equation  (21)  can  be  used  to  write  equations  (35)  and  (36)  in  terms  of  two 
unknowns : 


~ M<t>  - <t>  ) 

mm  r m 


m 


It(*’r  + k^r  * ^ + 


T_  + T 
I r 


(37) 


(38) 


where,  for  convenience,  an  effective  hull  inertia  torque  is  defined  as  follows: 


Tj  " m^  (-a^  sin  <J>  + &2  cos  d>)  - It  63  . 


(39) 


For  the  elevation  vibration,  the  coriesponding  equations  have  the  same 
form  a?  equations  (37)  and  (38)  except  that  the  variable  <J>  is  replaced  by  a. 
The  inertial  torque  for  this  case  will  have  the  form: 

A 

Tj.  * -n^KAj  + A^  cos  cOa^  + ( a^  cos  <J> +a2  sin  <f>)  (-B2  + A^  sin  a)] 

- It  (-61  sin  d>  + 02  cos  <#») . (40) 

Before  the  solution  to  equations  (37)  and  (38)  is  determined,  an  appro- 
priate expression  for  the  motor  torque  must  be  obtained.  This  expression 
will  be  developed  next.  Both  the  azimuth  and  the  elevation  motions  are 
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governed  by  the  same  type  of  servosyetem;  the  only  differences  are  in  the 
Inertial  and  geometrical  parameters,  and  in  the  gear  ratio  used  in  each  of 
these  motions.  A simplified  block  diagram  for  the  servocontrol  system  is 
shown  In  Figure  7. 


I 1 

A 

1 

i SSS 

T 

.75 

Inertia 

n 

n 

Transmission 

Figure  7.  AH-1J  Turret  Servosystem  Diagram 


In  the  simplification  of  the  actual  servodiagram,  certain  very  high-frequency 

effects  were  neglected.  Basically,  two  feedback  systems  are  present.  The 

first  feedback  signal  is  the  relative  motion  <j>r  or  that  is  measured  by 

use  of  a gear  position  resolver.  This  relative  motion  is  then  subtracted 

from  the  desired  control  signal  <j>  or  a . The  second  feedback  is  the  speed 

c c 

of  the  motor  that  is  measured  by  a tachometer  and  modified  by  a transfer 

function.  The  transfer  function  is  defined  in  terms  of  the  motor  angular 

displacements  $ or  a , and  the  resulting  signal  is  subtracted  from  the  dif- 
m m 

ference  '•£  the  control  ana  the  resolver  signals.  Actually,  before  the  sub- 
traction takes  place,  the  control  signal  is  amplified  by  the  factor  5.15, 
as  is  shown  in  Figure  7.  The  transfer  function  A,  defined  in  terms  of  the 


transformed  variable  s,  has  the  form: 
f 1 0.1s 

km  [T^  X 1+  .Is  X -067s 


. .0192s 

067  s + -j  n 


where  n is  the  gear  ratio.  For  the  azimuth  motion,  n = 620;  for  the  eleva- 
tion motion,  n - Bin.  Consequently,  in  terms  of  the  transformed  variable, 
the  motor  torque  for  the  azimuth  motion  is  given  by: 
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where  6 is  defined  as 


(A/ 535)* 


Writing  the  transfer  function  A over  a common  denominator  and  using  n 
620  results  in  the  following  equation: 

. _ .0029667s  + .0062692  s2  . 


Equation  (44)  19  given  in  terms  of  the  transformed  variable  s;  but,  for 
convenience,  this  equation  will  be  expressed  in  differential  form  by  the 
replacement  of  the  parameter  s with  the  time-differential  operator: 


+ .0062692 


0029667 


For  the  case  of  the  elevation  motion,  equation  (45)  is  modified  by  the 

replacement  of  d>  and  <J>  with  a and  a , and  by  the  multiplication  of  each 
s m s m 

term  on  the  right  side  with  the  factor  810/620.  The  motor  torque  at  any 
given  time  is  defined  from  equations  (42)  and  (45),  although  it  is  expressed 
In  the  form  of  a differential  equation. 

The  final  form  of  the  governing  equations  is  given  here  in  terms  of  the 
azimuth  deflection  variable  4> , and  similar  types  of  equations  are  valid  for 
the  elevation  parameter  a.  Consider  the  equation  of  motion  (37)  and  sub- 
stitute from  equation  (42)  to  get 


Equation  (46) , together  with  equations  (38)  and  (45) , represents  a system  of 
three  governing  equations  in  the  varinhlee  4-^,  and  *g.  For  the  case  of 
the  elevation  motion,  the  variable  $ is  changed  to  u and  the  right  side  of 
equation  (45)  is  altered  by  a multiplying  constant,  as  Indicated  previously. 

The  solution  to  the  three  governing  equations  will  be  obtained  numer- 
ically on  a computer.  To  do  this  calls  for  the  transformation  of  the  three 
second-order  equations  to  five  first-order  equations  by  the  definition  of 
the  following  five  variables: 


(47a) 


“ *r 


(47b) 


*3**. 

d<|i 

a,  * Jj “ 

^4  dt 


(47c) 

(47d) 


& <p 


8 • 


(47e) 


By  the  substitution  of  the  above 
(38),  (45),  and  (46),  the  set  of 


definitions  in  the  governing  equations 
five  first-order  equations  follows: 

(48a) 


ir-T  [-«*!  - +3>  - k8+2/|i2l  +TI  + Tr 

— — ■ A 
at 


(48b) 

(48c) 


d$4  1 

' *3>  + 401 -24n  <*c  - *i  " *5)]  (48d) 

m 

d* 

- - 104.  + 0.029667  4, 
dt  d 4 

+ 0.062692  [ y k(41  - 43)  + 401. 25n  (<J>c  - ♦1  - $5)].  (4Re) 

m 

Note  that,  in  the  last  of  equations  (48),  the  second  derivative  of 

<p  was  replaced  by  the  use  of  equation  (46).  Equations  (48)  are  solved 
n 

numerically,  and  the  solution  has  been  programmed  on  a digital  computer. 

The  FORTRAN  IV  source  program  for  the  Gun-Turret  Dynamic  Model  is 
listed  in  the  Appendix  ..  In  this  program,  file  1 is  the  input  file  that 
contains  the  data  passed  from  the  Helicopter  Vibration  Model.  A plot  rou- 
tine is  provided,  in  the  Gun-Turret  Dynamic  Model,  to  graph  the  dynamic 
responses  in  the  azimuth  and  elevation  directions  on  the  CALCOMP  Model  914 
plotter. 


22 


E.  Evaluation  of  the  Gun  Pointing  Errors 

The  results  of  the  Gun-Turret  Dynamic  Model,  described  in  the  previous 
section,  are  used  to  calculate  the  gun-pointing  errors.  These  errors  are 
calculated  from  the  statistical  properties  of  the  angular  dynamic  responses 
of  the  weapon.  The  statistical  properties  Include  the  mean  and  the  standard 
deviation  of  the  angular  motions  of  the  gun  barrel  about  two  mutually  per- 
pendicular axes.  With  the  U8e  of  the  initial  gun  barrel  direction  as  the 
reference  line,  the  gun-pointing  errors  can  be  defined  by  the  deviation  of 
the  gun  from  this  orientation.  This  deviation  can  be  calculated  by  the  ro- 
tation of  the  barrel  about  two  mutually  perpendicular  axes,  which,  in  this 
analysis,  are  chosen  to  be  z 3 and  z2,  The  choice  of  the  z 3 axis  is  a natural 
one  since  it  corresponds  to  the  elevation  motion  previously  defined  by 
equation  (20)  and  denoted  by  A a.  This  quantity  contains  two  parts,  the 
helicopter  rotation  ctQ  and  the  relative  rotation  a . The  helicopter  vibra- 
tion term  oq  can  be  obtained  from  the  data  generated  by  the  Helicopter 
Vibration  Model  and  it  is  related  to  the  angular  rotations  of  the  helicopter 
at  the  attachment  point  of  the  turret.  This  relation  can  be  obtained  by 
the  consideration  of  the  coordinate  transformation  matrices  for  the  three 
coordinate  systems,  shown  in  Figure  5.  On  the  assumption  that  the  coordi- 
nates x^  and-C^*  *2  ^3»  an<^  x3  ^2  by  small  angles  and  that 

the  second-order  terms  in  these  angles  are  negligible,  the  coordinate  trans- 
formations are  then  given  by 

[*.]  f-l  0 ol  ft.) 


0J 


cosa  cos<(>,  sina,  -cosot  sin<j» 
-sina  cos$,  cosa,  sina  simj> 


and 


4 J 


' ■■ 

r 

-cosa  cos  $ , -cosa  sin<t>,  slna 

*1 

sina  cos$,  sina  sin$  , cosa 

?2 

-Bin4  , cos$  , 0 

*3 

. . . ___ 

(51) 


Equations  (49)  to  (51)  represent  coordinate  transformations,  but  they  are 
also  valid  for  the  transformation  of  small  angular  rotations  about  the  various 
axes  of  the  three  coordinate  systems. 

Consider  now  the  rotation  about  the  z^  axis  that  has  been  defined  as 
the  elevation  motion  and  Is  given  by  equations  (20)  and  (22) . The  contri- 
bution c*r  Is  calculated  directly  by  the  Gun-Turret  Dynamic  Model,  and  the 
contribution  ctQ  is  related  to  the  helicopter  motion.  Using  the  third  of  the 
transformation  relations  given  in  equation  (51)  and  applying  it  to  the 
small  angular  rotation  vector  {6}  yields  the  result 


ao  - -01  sin  <j>  + 02  cos  $ , (52) 

since  0^,  6^,  and  0^  have  been  defined  as  the  helicopter  attachment  point 
rotations  about  the  axes.  Therefore,  the  total  angle  used  in  the  error 
analysis,  in  the  elevation  directidn,  is  given  by 


A a - -0^  sin$  + 0 ^ cos  + a . (53) 

Consider  now  the  rotation  about  the  z ^ axis,  which  is  denoted  by 
oud  is  divided  into  two  parts 

i 

An  - nc  + nr  , (54) 

where  no  is  the  contribution  from  the  helicopter  structural  vibration  and 

n is  the  contribution  from  the  turret  motion.  The  term  n is  obtained  by 
r o J 

use  of  the  second  transformation  in  equation  (51) , and  is  given  by 

no  " (6^  cos$  + ©j  sin$)  sina  + 0^  cosa.  (55) 

The  second  term  or  is  related  to  the  relative  azimuth  motion  <j>r  that  is 
defined  as  the  rotation  about  the  x2  or  the  axis.  Therefore,  the  contri- 
bution of  this  rotation  is  again  obtained  by  use  of  the  second  of  the  trans- 
formations in  equation  (51).  Consequently,  the  total  angle  to  be  used  in 
the  calculation  of  the  error  analysis  is  given  by 
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An  - (0^  cosij)  + 02  siniti)  sina  + (63  + (J>r)cosa.  (56) 

Equations  (53)  and  (56)  are  used  by  the  Gun-Turret  Dynamic  Model  to  compute 
the  mean  and  the  standard  deviation  of  the  gun  barrel  rotations  about  the 
two  mutually  perpendicular  axes.  Equation  (53)  represents  the  dynamic 

I 

error  in  the  elevation  direction,  and  equation  (56)  represents  the  dynamic 
error  about  an  axis  orthogonal  to  both  the  elevation  axis  and  the  barrel 
axis.  For  convenience,  the  dynamic  error  Aa  will  be  designated  the  turret 
elevation  error,  and  the  dynamic  error  An  will  be  designated  the  turret 
azimuth  error. 


F.  Recoil  Forces 

During  the  present  investigation,  a number  of  different  types  of  gun 
recoil  forces  were  utilized.  These  forces  correspond  to  different  caliber 
weapons  and  have  been  generated  by  different  methods  including  experimental 
and  simulated  approaches.  The  first  set  of  data  was  measured  experimentally 
(for  the  20mm  XM197  weapon)  on  the  six-degree-of-freedom  simulator,  at  the 
GEN  Keith  L.  Ware  Simulation  Center.  These  data  are  summarized  in  Table  1; 
and  the  plots  of  three  of  these  recoil  forces,  as  a function  of  time,  are 
given  in  Figures  8 to  10.  The  data  presented  in  these  figures  are  the  total 
recoil  force  transmitted  to  the  rigid  turret  structure  for  a twenty-round 
burst  sr.d  correspond  to  different  firing  rates  and  to  different  turret  posi- 
tions. In  each  case  the  data  points  were  recorded  at  one  millisecond  time 
intervals.  The  rates  of  fire  were  varied  from  500  to  817  shots  per  minute. 


TABLE  1 


List  of  Experimentally  Measured 
Standard  Recoil  Data 

Firing  Rate 


No. 

Test  No. 

Date 

spin 

No.  SI 

1 

2 

3/26/74 

670 

20 

2 

3 

3/26/74 

500 

20 

3 

1 

3/28/74 

669 

20 

4 

2 

3/28/74 

817 

20 

5 

3 

3/28/74 

500 

20 

6 

4 

3/28/74 

660 

6 

I 


I 


\ 


i 
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The  second  set  of  data  used  in  the  present  study  corresponds  to  con- 
stant recoil  data  that  were  provided  by  Honeywell  Incorporated.  These  data 
contain  eight  different  records  and  are  summarized  in  Table  2.  The  data 
again  represent  the  total  recoil  force  transmitted  to  the  turret  for  the 
20mm  gun  mounted  on  a constant  recoil  mechanism.  The  plots  of  these  data 
as  a function  of  time  are  shown  for  three  representative  cases  in  Figures 
11  to  13.  The  first  six  records  of  these  data  in  Table  2 correspond  to 
computer  simulated  results,  and  the  last  two  of  these  simulated  records 
contain  misfires  in  the  rounds  2 and  4 out  of  25.  The  rates  of  fire  for 
these  six  records  were  varied  from  675  to  820  shots  per  minute.  The  data 
records  7 and  8 correspond  to  experimentally  measured  constant  recoil 
forces  at  the  rate  of  approximately  750  shots  per  minute;  about  10  rounds 
were  fired.  The  simulated  and  the  measured  constant  recoil  data  are  notice- 
ably different;  nevertheless,  the  peak-to-peak  variation  in  the  recorded 
constant  recoil  data  is  well  below  that  of  the  test  data  recorded  from 
standard  spring  adapter  mounts. 


TABLE  2 

List  of  Honeywell  XM197  20mm 
Constant  Recoil  Data 


No. 

Code  Name 

Round  Impulse 
lb-sec 

Firing  Rate 
spm 

No.  Shots 

Type  of 
Data 

1 

RAANM 

34.56 

750 

25 

Simulated 

2 

RALN11 

34.56 

675 

25 

Simulated 

3 

RLLNM 

31.10 

675 

25 

Simulated 

4 

RHHNM 

38.02 

820 

25 

Simulated 

5 

RAAM2 

34.56 

750 

25 

Simulated 
with  misfire 
Rds.  2 and  4 

6 

RHHM2 

38.02 

820 

25 

Simulated 
with  misfire 
Rds.  2 and  4 

7 

DDGGT 

- 

750 

10 

Experimental 

8 

DDFFTA 

- 

750 

10 

Experimental 

29 
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HONEYWELL  CONSTANT  RECOIL  FORCE 


For  the  Bushmaster  25mm  and  the  AMCAWS  30mm  weapons,  which  were  also 
examined,  the  detailed  time  records  for  the  recoil  forces  were  unavailable. 
The  only  available  data  for  these  weapons  were  limited  to  the  peak  recoil 
force,  the  total  round  impulse,  and  the  rate  of  fire.  With  the  use  of  this 
information,  the  recoil  time  records  for  the  larger  caliber  weapons  were 
generated  by  the  modification  of  the  20mm  recoil  data.  The  modifications 
involved  the  adjustment  of  the  rate  of  fire  and  the  magnitude  of  the  peak 
recoil  forces.  The  rate  of  fire  was  adjusted  by  the  changing  of  the  time 
scale  of  the  20mm  records,  and  the  magnitude  was  altered  by  the  multipli- 
cation of  the  force  by  a conversion  factor.  The  actual  numbers  used  for 
these  changes  will  be  given  in  the  separate  sections  of  this  report  that 
deal  with  the  results  for  these  particular  weapons.  In  the  case  of  the 
25mm  gun,  both  the  regular  and  the  constant  recoil  forces  were  used  in  the 
analysis.  For  the  30mm  gun,  only  the  constant  recoil  force  was  used  in  the 
analysis,  because,  without  the  use  of  a constant  recoil  device,  the  peak  re- 
coil force  for  this  weapon  would  be  too  large  for  the  practical  operation 
on  the  AH-1G  helicopter. 


III.  ANALYSIS  OF  THE  20mm  XM197  WEAPON 


A.  Weapon  Deacriptlon 

The  XM197  weapon  Is  a three-barrel  20mm  automatic  gun.  The  rate  of 
fire  of  this  weapon  can  be  varied;  however,  for  helicopter  applications, 
the  weapon  firing  rate  is  normally  selected  at  about  750  shots  per  minute. 
The  gun  is  supported  by  the  XM97  turret,  which  controls  the  motion  of  the 
gun  in  the  azimuth  and  in  the  elevation  directions.  The  azimuth  orienta- 
tion can  be  varied  from  0°  (front  of  helicopter)  to  + 110°.  The  elevation 
motion,  measured  positive  up  from  the  horizontal  plane  of  the  helicopter, 
can  be  varied  between  5°  and  about  -50°.  The  various  servocontrol  para- 
meters have  been  defined  in  the  Section  II-D  in  which  the  Gun-Turret  Dynamic 
Model  is  treated.  Additional  parameters,  by  which  the  physical  properties 
of  the  weapon  are  defined,  consist  of  the  inertial  parameters  and  the 
positions  of  the  centers  of  mass  of  the  moving  parts  of  the  turret  and  the 
gun.  These  positions  are  defined  in  Figure  4 by  the  parameters  A^,  A^,  A^, 
B^,  and  By  an  examination  of  the  specifications  of  the  gun  and  the 

turret,  these  parameters  were  chosen  to  have  the  following  values: 

A^  ■ 2.4",  * 2.4",  A^  “ 10",  ■ 6",  and  * 6".  The  two  masses 

that  are  considered  in  the  analysis  are  the  mass  of  the  moving  part  of  the 
turret  and  the  mass  of  the  gun.  These  masses  were  estimated  as  130  and 
170  pounds,  respectively.  The  moments  of  inertia  of  the  rotating  compon- 
ents of  the  assembled  gun  and  turret  about  the  azimuth  and  the  elevation 

2 

pivots  are  20.7  and  13.7  slug-ft  , respectively.  The  axial  moments  of  in- 

-4 

ertla  for  the  two  electric  drive  motors  are  identical  and  are  2.4  x 10 
2 

slug-ft  . The  gear  ratios  for  the  azimuth  and  the  elevation  motions  are 
620  and  810,  respectively;  the  transmission  stiffness  in  each  drive  was 
chosen  as  5 x 10^  lb-ft/rad.  The  structural  damping  factor  of  the  drive 
system  was  taken  as  0.0001. 

The  various  weapon  parameters  defined  in  this  section  were  used  in  the 
numerical  analysis  of  the  dynamic  motion  of  the  gun  mounted  on  the  AH-1J  or 
on  the  AH-1G  helicopter  frame.  The  results  of  these  calculations  are  dis- 
cussed in  the  next  two  sections  of  this  report. 

B.  Analysis  of  the  AH-1J  Helicopter 

The  first  phase  of  investigation  involved  the  analysis  of  the  pointing 
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errors  of  the  20mm  XM197  gun  attached  to  the  AH-lJ  helicopter  frame.  The 
analysis  was  performed  for  five  different  recoil  force  histories.  Three 
of  these  histories  correspond  to  experimentally  measured  forces  and  are 
listed  as  forces  No.  3,4,  and  5 in  Table  1.  These  forces  were  chosen  to 
cover  a broad  range  of  firing  rates;  these  rates  correspond  to  669,  817, 
and  500  spm.  The  remaining  two  forces  are  the  Honeywell  simulated  constant 
recoil  forces  obtained  by  a mechanism  model  of  that  device.  These  forces 
are  identified  by  the  codes  RAANM  and  RAAM2;  the  first  code  refers  to  the 
data  with  no  misfires  and  the  second  refers  to  the  data  with  two  misfires. 
These  are  two  of  the  forces  described  in  Table  2 and  plots  of  these  two 
forces  are  provided  in  Figures  11  and  12.  Both  of  these  forces  had  a 
firing  rate  of  750  spm. 

For  each  of  the  five  forces,  the  dynamic  helicopter  and  turret  models 
were  analyzed  for  12  different  angular  orientations  of  the  gun  barrel.  The 
orientations  include  combinations  of  the  four  different  azimuth  angles  of 
0,  22.5,  45,  and  67.5  degrees,  and  the  three  different  elevation  angles  of 
-5,  -25,  and  -45  degrees.  In  each  case,  the  results  of  the  calculations 
are  presented  in  terms  of  the  mean  and  the  standard  deviation  of  the  an- 
gular motion  error  in  both  the  azimuth  and  the  elevation  directions.  The 
results  are  presented  in  Tables  3 to  7,  in  which  the  values  of  the  statisti- 
cal error  parameters  are  giver,  in  units  of  milliradians.  The  most  impor- 
tant statistical  parameter  is  the  standard  deviation  because  it  represents 
the  theoretical  gun  scatter.  The  maximum  standard  deviation  er-or,  from 
the  recoil  force  excitation,  occurs  when  the  gun  is  fired  from  the  side  of 
the  helicopter,  that  is,  at  the  maximum  azimuth  angle.  This  result  was 
expected  since,  for  large  azimuth  orientations,  the  recoil  torque  imparts 
a predominant  torsional  mode  of  vibraticn  in  the  helicopter  frame.  In 
the  comparison  of  the  responses  of  the  standard  recoil  forces  (Tablos  3 to 
5)  with  the  corresponding  responses  of  the  constant  recoil  forces  (Tables 
6 and  7),  the  conclusion  can  be  made  that  the  constant  recoil  mechanism  is 
effective  in  the  reduction  of  the  dynamic  motion,  even  when  misfires  occur. 
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TABLE  3 

AH-1J  Analysis  Results  for  Teat  3 
3/28/74  (300sr>m)  Ftnndufd  Recoil 


AH-1J  Analysis  Results  for  Honeywell 
Constant  Recoil  Force  (Coae  RAANM) 


AH-1J  Analytic  Retultt  for  Hontyvtll 
Con* tart  Recoil  Force  (Code  RAAM2) 


'•  Analysis  of  the  AH-1G  Helicopter 


The  AH-1J  helicopter  model  is  sufficiently  accurate  to  predict  the 
dynamic  response  to  weapon  recoil  forces.  The  model,  however,  contains 
some  deficiencies  that  needed  to  be  corrected  to  improve  the  accuracy  of 
the  results.  In  particular,  the  dynamic  predictions  of  the  AH-1J  model  are 
valid  only  for  frequencies  below  20  Hz;  however,  the  weapon  recoil  force 
data  usually  include  large  amplitude  Fourier  series  components  that  corre- 
spond to  frequencies  above  20  Hz.  To  improve  the  accuracy  of  the  model, 
a more  detailed  representation  of  the  AH-1  helicopter  structure  was  there- 
fore required,  which  would  yield  accurate  solutions  in  the  frequency  range 
from  zero  to  about  40  Hz.  Another  problem  encountered  with  the  AH-1J 
model  is  that  the  elastic  axis  model  of  the  forward  fuselage  strucutre  does 
not  represent  the  torsional  stiffness  of  the  actual  helicopter  structure  in 
sufficient  detail.  For  these  reasons.  Bell  Helicopter  Company  was  contracted 
to  provide  a NASTRAN  model  of  the  AK-1G  helicopter  that  includes  a detailed 
representation  of  the  fuselage  structure.  The  resulting  model  contains 


over  1700  independent  degrees  of  freedom  compared  with  about  300  for  the 
AH-1J  model.  The  free  vibration  solution  for  the  AH-1G  model  is  provided 
in  Reference  2.  This  solution  includes  the  generalized  mass  and  stiffness 
matrices,  the  natural  frequencies,  and  the  matrix  of  mode  shapes.  Again, 
only  the  30  lowest  natural  frequencies  and  the  corresponding  mode  shapes 
were  ised  in  the  forced-vibration  analysis  with  the  Helicopter  Vibration 
Model. 

In  the  subsequent  analyses  with  the  Helicopter  Vibration  Model,  the 
data  from  the  AH-1G  NASTRAN  model  were  used  in  place  of  the  data  from  the 
AH-1J  model.  Before  any  new  problems  were  solved,  however,  some  comparisons 
of  results  obtained  with  the  AH-1J  and  with  the  AH-1G  models  were  made  for 
the  dynamic  response  of  the  20mm  XM197  gun.  The  calculations  were  per- 
formed for  several  of  the  20mm  recoil  forces  listed  in  Tabies  1 and  2;  and 
the  results  obtained  from  the  AH-1G  model  version  were  compared  with  the 
corresponding  results  obtained  previously  from  the  AH-1J  model  version. 

The  analysis  was  performed  for  each  combination  of  the  four  initial  azimuth 
angles  (0°,  22.5°,  45°,  and  67.5°)  and  for  the  three  initial  elevation  an- 
gles (-5°,  -25°,  and  -45°).  Again,  the  dynamic  responses  of  the  gun  bar- 
rel in  the  azimuth  and  elevation  directions  were  determined  during  the 
burst  interval;  and,  from  these  results,  the  statistical  means  and  standard 
deviations  were  calculated. 

An  example  of  the  dynamic  response  of  the  20mm  weapon  mounted  on  the 
AH-1G  airframe  is  provided  in  Figures  14  and  15  for  one  set  of  conditions. 
The  case  presented  is  that  of  the  standard  20mm  recoil  force  with  a 670  spm 
firing  rate  and  an  initial  gun  position  of  45°  in  azimuth  and  -5°  in  eleva- 
tion. The  recoil  force  excitation  is  plotted  in  Figure  8,  and  the  azimuth 
and  elevation  vibration  patterns  are  plotted  in  Figures  14  and  15,  respec- 
tively. Also,  the  statistical  mean  and  standard  deviation  parameters  de- 
termined from  the  vibration  patterns  are  given  in  Figures  14  and  15.  Note 
that,  in  this  particular  cas^  the  peak  amplitudes  of  the  elevation  vibration 
are  significantly  larger  than  those  of  the  azimuth  vibration. 


Figure  15.  Dynamic  Elevation  Response  to  the  670  sum  Standard  20mm  XM197  Recoil 
Force  at  the  45°  Azimuth  and  -5°  Elevation  Position 
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The  dynamics  of  the  elevation  mechanism  are  affected  to  a large  ex- 
tent by  both  the  azimuth  orientation  of  the  turret  and  the  mass  unbalance 
of  the  weapon,  as  defined  by  the  parameter  in  Figure  4.  At  large  azi- 
muth orientations,  the  recoil  force  induces  predominantly  large  torsional 
vibrations  of  the  helicopter  structure.  The  torsional  vibrations  of  the 
turret  platform  directly  affect  the  dynamic  motions  of  the  elevation  mech- 
anism; and,  furthermore,  the  platform  accelerations  induce  large  inertial 
torques  about  the  elevation  trunnion.  Ideally,  the  turret  should  be  de- 
signed in  such  a manner  that  the  center  of  mass  of  the  weapon  lies  along 
the  intersection  of  the  azimuth  axis  and  the  elevation  trunnion  axis.  For 
practical  reasons,  however,  mass  unbalances  cannot  always  be  eliminated. 

A second  example  of  the  dynamic  azimuth  and  elevation  responses  is 
provided  in  Figures  16  and  17,  respectively.  In  this  example,  the  applied 
excitation  force  is  the  Honeywell  constant  recoil  force  that  is  specified 
by  the  code  DDFFTA  and  is  shown  in  Figure  13.  The  initial  turret  orienta- 
tion was  the  same  as  that  of  the  previous  example,  namely,  45°  in  azimuth 
and  -5°  in  elevation.  Again,  the  predominant  vibration  pattern  occurs  in 
the  elevation  mechanism;  however,  the  magnitudes  of  the  statistical  errors 
are  much  less  than  the  corresponding  results  in  the  previous  example.  For 
instance,  the  standard  deviation  error  in  the  elevation  motion  is  reduced 
from  a value  of  9.28  to  a value  of  1.73,  by  the  use  of  the  constant  recoil 
mount . 

The  sensitivity  of  the  gun-pointing  error  to  the  initial  turret  azi- 
muth orientation  is  illustrated  by  a set  of  results  presented  in  Figures 
18,  19,  and  20  that  correspond  to  three  types  of  weapon  recoil  forces.  In 
these  figures,  each  of  the  four  statistical  errors  is  plotted  versus  the 
turret  azimuth  orientation  and  the  initial  turret  elevation  angle  is  set 
at  -25°.  The  results  for  the  817  spm  standard  recoil  force  are  shown  in 
Figure  18;  the  results  for  the  Honeywell  constant  recoil  mechanism  model 
data  (code  RAANM)  are  shown  in  Figure  19;  and  the  results  for  the  Honey- 
well constant  recoil  test  data  (code  DDFFTA)  are  shown  in  Figure  20.  In 
each  case,  the  mean  and  standard  deviation  errors  associated  with  the  ele- 
vation vibration  are  the  predominant  errors,  especially  when  the  initial 
turret  azimuth  position  is  greater  than  about  30°.  At  the  -25°  elevation 
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Figure  17.  Dynamic  Elevation  Response  to  the  Honeywell  Constant  Recoil  Force  (Code 
DDFFTA)  at  the  45°  Azimuth  and  -5  Elevation  Position 


Figure  18.  Statistical  Gun  Vibration  Errors  vs.  Turret  Azimuth  Angle  for  the  Test 
2-3/28/74  Recoil  Force,  -25°  Turret  Elevation  Angle 


setting,  the  elevation  error  parameters  generally  become  larger  as  the 
turret  azimuth  angle  increases.  The  azimuth  errors,  however,  are  rela- 
tively insensitive  to  the  turret  azimuth  position.  As  expected,  the  con- 
stant recoil  device  had  a noticeable  effect  in  the  reduction  of  the  magni- 
tudes of  the  statistical  errors. 

The  improvement  in  the  gun-pointing  accuracy  as  a result  of  the  con- 
stant recoil  device  is  further  illustrated  in  Figure  21.  In  this  figure, 
the  standard  deviation  error  in  elevation  is  plotted  versus  the  initial 
turret  azimuth  angle  for  the  standard  recoil  force  (Test  2-3/26/74),  for 
the  Honeywell  corutant  recoil  terc  data  (code  DDFFTA) , and  for  the  Honey- 
well constant  recoil  mechanism  model  data  (code  RAANM).  The  case  repre- 
sented agaia  corresponds  to  an  initial  elevation  angle  of  -25°.  The 
standard  deviation  errcr  function  for  the  constant  recoil  test  data  (code 
DDFFTA)  lies  slightly  above  the  error  function  for  the  idealized  mechanism 
model  data  (code  RAANM);  nevertheless,  both  graphs  for  the  constant  recoil 
forces  lie  well  below  the  graph  for  the  standard  recoil  force.  In  this 
particular  example,  the  firing  rate  for  the  standard  recoil  force  was  670 
shots  per  minute,  and  the  firing  rate  for  both  sets  of  constant  recoil 
data  was  750  shots  per  minute. 

The  sensitivity  of  the  gun-pointing  error  to  the  weapon  firing  rate 
was  also  investigated  for  the  20mm  XM197  weapon.  The  weapon  recoil  data 
listed  in  Tables  1 and  2 were  recorded  at  several  different  firing  rates 
that  ranged  from  about  500  spm  to  820  spm.  In  addition,  the  constant  re- 
coil mechanism  model  data  were  recorded  with  slightly  different  round 
impulses.  Each  of  these  recoil  forces  was  used  in  the  sensitivity  investi- 
gations. For  each  particular  turret  orientation,  the  statistical  error 
predictions  obtained  for  each  of  the  different  constant  recoil  forces  did 
not  vary  appreciably,  regardless  of  the  firing  rate  or  the  round  impulse. 
Only  the  occurrence  of  misfires  in  the  constant  recoil  force  data  affected 
the  results  by  a noticeable  degree. 

An  examination  of  the  results  obtained  for  the  standard  recoil  forces 
reveals,  however,  that  the  weapon  firing  rate  significantly  affects  the 
vibration  patterns.  In  Figure  22,  the  standard  deviation  error  in  elevation 
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is  plotted  versus  the  turret  azimuth  angle  for  three  standard  recoil  force 
excitations.  The  weapon  firing  rates  correspond  to  500,  670,  and  817  shots 
per  minute.  The  weapon  firing  rate  not  only  affects  the  magnitudes  of  the 
statistical  errors  at  each  particular  turret  orientation,  but  also  affects 
the  functional  dependency  of  the  statistical  errors  on  the  turret  azimuth 
position  angle.  These  observations  were  anticipated  because  the  vibration 
amplitudes  of  each  natural  mode  of  the  helicopter  are  dependent  upon  both 
the  frequency  spectrum  and  the  direction  of  the  applied  force.  Based  on 
these  results,  the  conclusion  is  made  that  the  dynamic  response  of  a given 
weapon  system  to  the  standard  weapon  recoil  loads  can  be  significantly 
improved  by  the  adjustment  of  the  firing  rate.  If,  however,  the  weapon  system 
includes  a constant  recoil  mount,  then  the  weapon  firing  rate  will  not  appre- 
ciably affect  the  dynamic  response. 

The  results  of  the  AH-1G  structural  sensitivity  analysis  for  the 
20mm  XM197  gun  are  summarized  in  Tables  8 to  12.  Tables  8,  9,  and  10 
contain  the  statistical  error  predictions  for  three  of  the  standard  recoil 
forces  that  are  listed  in  Table  1.  Tables  11  and  12  contain  the  results 
for  the  Honeywell  constant  recoil  forces  that  are  designated  by  the  codes 
RAANM  and  DDFFTA.  A comparison  of  the  results  summarized  in  these  five 
tables  reveals  again  that  the  Honeywell  constant  recoil  mechanism  is  very 
effective  in  the  reduction  of  the  angular  gun  vibrations.  The  constant 
recoil  mount  is  particularly  effective  in  the  reduction  of  the  standard 
deviation  error  in  the  elevation  mechanism.  For  example,  the  standard  de- 
viation error  in  elevation  is  greater  than  13  milliradians , for  the  case 
in  which  the  initial  elevation  angle  is  -5°  and  a standard  recoil  force 
is  applied  at  a 500  spin  firing  rate.  These  results  are  shown  in  Table  10. 

The  corresponding  results  for  the  Honeywell  constant  recoil  force  DDFFTA, 
which  are  shown  in  Table  12,  range  from  0.84  to  2.90  milliradians.  In 
this  case,  the  use  of  the  constant  recoil  mount  resulted  in  a significant 
improvement  in  the  gun-pointing  accuracy.  Moreover,  in  all  of  the  other 
cases  analyzed,  large  improvements  in  the  standard  deviation  errors  re- 
sulted from  the  use  of  the  constant  recoil  mount.  The  constant  recoil 
mount,  however,  is  not  very  effective  in  the  reduction  of  the  mean  errors 
in  gun  pointing.  The  mean  error  appears  to  depend  primarily  on  the  mean 
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level  of  the  weapon  recoil  force  instead  of  the  peak-to-peak  variation 
in  the  force. 

The  results  presented  in  Tables  8 to  12,  which  were  obtained  using  the 
AH-1G  model  version,  can  be  compared  directly  with  the  results  shown  in 
Tables  3 to  7,  which  were  obtained  using  the  AH-lJ  model  version.  In  both 
the  AH-U  and  the  AH-1G  analyses,  the  applied  excitation  forces  consisted 
of  a selection  of  three  standard  recoil  forces,  with  firing  rates  of  approx- 
imately 500,  670,  and  817  shots  per  minute,  and  two  Honeywell  constant  re- 
coil forces.  An  examination  of  the  AH-lJ  and  the  An-IG  model  results  for 
the  same  type  of  recoil  force  shows  that  both  models  give  results  of  the 
same  order  of  magnitude.  The  differences  in  the  predictions  of  the  two 
models  are  most  pronounced  in  the  variations  of  the  standard  deviation  er- 
rors with  the  turret  azimuth  position.  With  the  AH-1G  model,  these  errors 
are  relatively  large  for  all  turret  positions;  whereas,  with  the  AH-lJ 
model,  these  errors  generally  increase  in  magnitude  from  a relatively  small 
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azimuth  position  angle  increases  from  0°  to  67.5°. 

This  discrepancy  is  a consequence  of  the  much  more  accurate  torsional 
stiffness  representation  in  the  AH-1G  model.  In  this  model,  the  torsional 
and  the  bending  mode  shapes  for  the  forward  fuselage  structure  are  strongly 
coupled;  however,  in  the  AH-lJ  model,  the  bending  and  the  torsional  modes 
are  primarily  uncoupled.  Hence,  the  AH-lJ  model  exhibits  a greater  sensi- 
tivity to  the  turret  azimuth  position  than  the  AH-1G  model. 

D.  The  XM97  Turret  Servocontrol  Modifications 

The  results  presented  in  the  two  previous  sections  of  this  report  showed 
that,  under  certain  conditions,  the  dynamic  motion  of  the  gun  would  produce 
relatively  large  gun  errors,  especially  the  error  in  the  standard  deviation 
of  the  elevation  motion.  Closer  examination  cf  the  computer  results  revealed 
that  the  dynamic  motion  was  predominantly  of  one  frequency,  which  is  closely 
related  to  the  natural  frequency  of  the  turret  elevation  mechanism.  Con- 
sequently, any  reduction  of  this  motion  would  involve  a reduction  of  the 
natural  frequency  response  of  the  turret.  One  objective  of  the  overall 
study  was  to  determine  whether  such  a reduction  is  possible  by  the  modifi- 
cation of  the  existing  XM97  servocontrols,  with  the  addition  of  a rate  gyro 
sensor  attached  at  one  of  several  possible  locations  in  the  svstem.  Al- 
together, four  different  modifications  of  the  controls  have  been  tried; 
these  are  labeled  as  Methods  1 to  4. 

Method  1 

In  Method  1,  the  rate  gyro  is  placed  on  the  gun  and  the  signal  from 
this  gyro  is  fed  back  into  the  servoloop.  The  corresponding  modification 
of  the  controlsystem  is  illustrated  in  Figure  23  and  the  unmodified  system 
is  shown  in  Figure  7.  The  rate  gyro  modification  is  represented  by  a loop 
containing  the  total  angular  displacement  of  the  gun,  or  its  time  derivative. 
This  displacement  is  denoted  by  <f>r  + <j>Q  for  the  azimuth  motion.  The  ex- 
pressions for  the  elevation  motion  are  similar;  therefore,  the  development 
of  duplicate  relations  will  not  be  necessary.  The  term  ^ is  the  feedback 
signal  from  the  azimuth  gear  position  resolver.  This  term  represents  the 
relative  displacement  between  the  turret  attachment  platform  and  the  ro- 
tating azimuth  mechanism.  The  term  <j>Q  is  the  absolute  displacement  of  the 
turret  attachment  platform,  in  the  azimuth  direction.  Some  additional 
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Figure  23.  Modified  AH-1J  Turret  Servosystem  Diagram 


a^rcra^  instrumentation,  such  as  a rate  gyro  or  some  other  tvoe  of  iner- 
tial sensor,  would  be  required  to  measure  the  displacement  <}> q . 

In  the  feedback  loop,  the  total  angular  displacement  signal  is 
processed  through  a transfer  function,  Q.  To  choose  an  expression  for 
the  transfer  function,  Q,  requires  a consideration  of  the  basic  purpose 
of  this  function.  As  mentioned  in  the  beginning  of  this  section,  the 
main  angular  response  of  the  gun  and  turret  occurs  at  a predominant  fre- 
quency related  to  the  natural  frequency  of  the  elevation  drive  system. 
Consequently,  in  the  modification  of  the  control  system,  a transfer  func- 
tion is  desirable  that  would  have  certain  frequency  characteristics. 
Specifically,  the  transfer  function  would  amplify  the  signal  in  the  con- 
trol loop  by  the  maximum  amount,  when  the  predominant  frequency  of  this 
signal  is  at  or  near  the  natural  frequency  of  the  turret.  Furthermore, 
the  effect  of  the  transfer  function,  Q,  should  decrease  at  frequencies 
that  are  above  and  below  this  natural  frequency.  The  Laplace  transforma- 
tion of  the  expression  for  the  feedback  signal  is  given  by 

<t>g(s)  = Q(s)  [<f>r(s)  + 4>0(s)]  (57a) 

and  a transfer  function  Q(g)  that  would  have  these  desired  characteristics 
might  be  represented  by 
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where  s is  the  transform  variable,  G is  an  amplifying  constant,  and 
is  a time  parameter  used  to  define  the  maximum  amplification  frequency. 
The  parameter  is  chosen  to  be  very  close  to  the  natural  period  of  vi- 
bration of  the  turret  drive  system.  An  illustration  of  the  transfer 
function,  Q,  is  provided  in  Figure  24,  in  which  the  graph  of  Q (S)  is 
plotted  on  a semi-log  scale  for  a particular  choice  of  the  parameters 
and  G.  The  actual  values  that  were  chosen  for  the  parameters  T^  and 
G,  in  the  modified  gun-turret  analysis,  are  presented  later  in  this 
section. 

The  additional  feedback  signal  in  the  modified  servocontrol  loop, 

namely  <f)  , is  used  to  modify  the  drive  motor  torque.  With  reference  to 
S 

equation  (42),  the  new  expression  for  the  motor  torque  T^  becomes: 

T = 401. 25n  (*  - <J>  - <J>  - <p  ) (58) 
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For  the  incorporation  of  this  modification  into  the  gun-turret,  equation 
(57a)  is  first  transformed  in  the  time  domain  and  the  resulting  differen- 
tial equation  is  given  by 


dt2 


2T. 


d<j> 

dt 


d(j> 
= c 

dt 


d<j> 

i 

dt 


(59) 


This  equation  is  then  written  in  terms  of  two  first-order  equations  by  the 
definition  of  two  new  variables  <j>,  and  <(>.,  where 

D / 


d»6 

dt 


(60) 

(61) 


With  the  use  of  these  definitions,  equation  (59)  can  be  written: 


^7  = 1 d*o  + — 2T-.  4> 

dt  2 [_dt 


*6] 


(62) 


The  two  first-order  equations  (61)  and  (62)  are  added  to  the  set  of  equa- 
tions (48)  to  modify  the  dynamic  model.  Accordingly,  the  Gun-Turret  Model 
computer  program  was  modified  and  then  an  analysis  was  performed  to  check 
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TRRNSFER  FUNCTION  QfS)  VS. 


Figure  24.  Transfer  Function  O(S)  For  The  Modified  Turret  Servocontrol  System 


out  the  effect  of  this  modification.  An  examination  of  the  response  of 

the  modified  system  revealed  that,  for  this  study,  the  term  ^o  is  small 

d a dt 

compared  with  the  term  vr  in  equation  (59) , and  therefore  this  term  was 

neglected.  dt 

A specific  angular  orientation  was  chosen  to  check  the  servocontrol 
modification  and  this  orientation  corresponds  to  an  azimuth  angle  of  67.5° 
and  an  elevation  angle  of  -5°.  Usually,  the  largest  gun-pointing  error 
was  obtained  ac  this  particular  turret  orientation.  In  Section  III-B  the 
results  of  the  analysis  with  the  AH-1J  Helicopter  Vibration  Model  and  the 
unmodified  Turret  Dynamic  Model  are  presented.  One  specific  set  of  results 
is  presented  in  Table  3 for  the  standard  20mm  recoil  force  with  a firing 
rate  of  669  shots  per  minute.  The  predominant  turret  vibration,  in  this 
analysis,  is  characterized  by  che  standard  deviation  of  the  elevation  mo- 
tion t and  the  value  obtained  for  this  parameter  at  the  chosen  turret  orien- 
tation is  17.5  milliradians.  The  amount  of  improvement  obtained  in  this 
standard  deviation  is,  therefore,  one  way  to  estimate  the  effectiveness  of 
the  present  modification  to  the  control  system.  Consequently,  in  this 
section,  a ratio  of  the  new  value  to  the  previous  value  of  this  standard 
deviation  is  presented  as  a measure  of  the  effectiveness  of  the  servocontrol 
modification  . This  ratio  is  denoted  by  the  term  "modified  ratio". 

In  the  present  modification,  the  two  additional  servocontrol  parameters 
that  are  defined  in  equation  (57b)  are  G and  T^.  So  that  a reasonable  value 
for  T^  could  be  obtained,  the  natural  frequency  of  the  drive  system  was 

examined  in  the  azimuth  and  in  the  elevation  modes.  These  frequencies  were 

2 

calculated  from  the  servoparameters  and  have  been  found  to  be  1.91  x 10 
2 

and  1.55  x 10  radians  per  second  for  the  elevation  and  for  the  azimuth 

motions,  respectively.  Consequently,  the  parameter  T1  for  the  elevation 

motion  was  chosen  to  be  5.6  x 10  , by  the  requirement  that  ioT^  = 1,  in 

which  a)  is  the  natural  frequency  of  the  drive  system.  However,  since  this 

is  only  an  approximate  estimate,  several  different  values  of  the  parameter 

\ 

T^  were  examined.  The  numerical  results  are  shown  in  Cases  a.  and  b.  of 

Table  13.  The  numerical  results  in  Case  a.  show  the  modified  ratio  for 
_3 

T^  * 5.6  x 10  and  for  different  values  of  the  parameter  C defined  in 
equation  (57b).  Note  that,  in  Cas.;  a^  G is  varied  over  a range  of  0.1  to 
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Case  e. 


G 

Modified  Ratio 

2 * 

10  4 

1.25 

0.5  x 

10 A. 

1.04 

-0.5  x 

10-4 

1.14 

-2  x 

10  4 

1.45 

-0.1  and  that  a minimum  value  for  the  modified  ratio  of  1.18  is  obtained 
for  G = 0.05.  Increasing  or  decreasing  G from  this  value  produced  larger 
motions.  Since  none  of  the  values  of  G produced  a value  below  1.0  for  the 
modified  ratio,  the  modified  control  system,  therefore,  did  not  improve  the 
dynamic  motion.  The  sensitivity  cf  the  modified  ratio  to  the  variable  was 
then  examined  and  the  results  cf  this  investigation  are  presented  in  Case  b. 
of  Table  13  for  G = 0.05  and  for  different  values  of  T^.  Note  that  in  Case 
b.  the  modified  ratio  reaches  approximately  a constant  value  of  one  for 
large  values  of  T^.  Similar  results  were  also  obtained  for  G = 0.1.  On 
the  basis  of  these  results,  the  conclusion  is  made  that  the  present  modifi- 
cation would  not,  in  most  instances,  improve  the  dynamic  response  of  the 
control  system. 
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Method  2 

The  second  type  of  modification  of  the  control  system,  which  was  ex- 
amined, consisted  of  a similar  approach  to  Method  1,  but  with  a slight  dif- 
ference in  the  transfer  function  Q.  The  new  transfer  function  is  given  by 

Q(s)  - G/(l  + (63) 

Note  that  this  function  is  similar  to  the  one  in  equation  (57b)  except  that 
the  parameter  s is  omitted  from  the  numerator.  The  effect  of  this  change 
is  to  produce  a response  that,  for  sinusoidal  inputs,  will  differ  by  a 
90-degree  phase  angle  from  that  of  equation  (57b) . 

With  the  use  of  equation  (63)  instead  of  equation  (57b) , the  gover- 
ning equation  (62)  is  modified  to  read 

d*7  - — i-  (<J>  +<!>,-  2T  - ♦,).  (64) 

dt  T^  ° X 1/0 

Equations  (61)  and  (64)  are  now  the  two  additional  first-order  equations 
that  are  used  to  modify  the  existing  Gun-Turret  Model. 

As  in  the  case  of  Method  1,  the  new  model  was  applied  to  the  case  in 
which  the  azimuth  and  the  elevation  angles  are  given  by  67.5  and  -5°,  re- 
spectively. The  modified  ratio  was  determined  for  values  of  G equal  to 

_3 

10  and  -10  and  for  the  parameter  T^  equal  to  5.6  x 10  seconds.  The  re- 
sults of  this  investigation  are  shown  in  Case  c.  of  Table  13.  Note  from 
these  results  that  choosing  either  value  for  G increases  the  modified  ratio. 
In  fact,  the  large  result  that  occurs  for  G = -10  indicates  an  unstable 
system.  Other  values  of  were  also  examined,  but  were  found  to  give  simi- 
lar results  $ that  is,  no  reduction  in  the  dynamic  motion  was  obtained.  On 
the  basis  of  these  results,  the  conclusion  is  made  that  this  modification 
does  not  improve  the  response  of  the  control  system. 

Method  3 

In  this  approach  and  in  Method  4 to  be  described  next,  a slightly  dif- 
ferent technique  was  used  to  modify  the  control  system.  Specifically,  the 
motor  torque  is  modified  by  a function  proportional  to  angular  velocity  of 
the  gun.  Consequently,  equation  (42)  is  modified  to  read: 
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T - 401.25  n(*c  — <|>r  — 4>g  — G 4>r>  . (65) 

The  gun-turret  model  was  modified  accordingly  and  the  numerical  results 

for  the  modified  model  are  shown  in  Case  d.  of  Table  13  for  different 

values  of  the  parameter  G.  These  results  show  that  a minimum  value  of  the 

_2 

modified  ratio  is  obtained  if  G is  set  equal  to  -1  x 10  and,  on  either 
side  of  this  value,  the  modified  ratio  increases.  Again,  a modified  ratio 
of  less  than  one  was  not  obtained}  consequently  this  approach  also  does 
not  bring  about  a reduction  in  the  gun  motion. 

Method  4 

In  Method  4,  the  motor  torque  is  modified  by  the  gun-barrel  accelera- 
tion feedback  rather  than  by  the  velocity  feedback.  Using  a rate  gyro 
with  this  approach  would  require  a differentiating  circuit.  With  this  ap- 
proach, the  motor  torque  is  given  by 

T = 401. 25n  (<|>  - $ - <J>  - G$  ) . (66) 
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The  results  of  this  modification  are  summarized  in  Case  e.  of  Table  13. 

-4 

A minimum  value  for  the  modified  ratio  was  obtained  for  G = 0.5  x 10  , but 

again  none  of  the  values  of  G were  found  to  give  a modified  ratio  of  less 
than  one. 

On  the  basis  of  the  results  presented  in  this  section,  the  conclusion 
is  made  that  none  of  the  simple  modifications  of  the  existing  control  sys- 
tem provide  any  effective  reduction  ir.  the  motion  of  the  gun.  What  the 
numerical  results  do  indicate  is  that  a aifficulty  exists  in  the  generating 
of  enough  motor  torque  to  correct  the  elastic  motion  in  the  transmission. 
This  difficulty  may  be  caused  by  an  upper  bound  on  the  motor  torque  due  to 
amplifier  saturation  and,  in  all  the  modified  model  versions  described  here, 
this  saturation  effect  was  included.  All  the  modifications  described  up 
to  this  point  can  be  considered  as  external  to  the  control  system  in  the 
sense  that  the  only  change  is  an  addition  of  an  extra  feedback  signal  to  the 
control  system.  Hc.*aver,  more  fundamental  changes  could  also  be  made  to  the 
controls  and,  therefore,  the  question  arises  about  the  effectiveness  of  such 
changes  in  reducing  the  dynamic  motion.  One  method  that  has  been  suggested 
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to  Improve  the  response  of  the  control  system  Is  the  use  of  optimal  control 
techniques  In  the  design  of  the  servocontrol  system.  Additional  studies 
are  therefore  being  planned  to  investigate  the  dynamic  response  of  an  op- 
timally controlled  turret. 
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IV.  ANALYSIS  OF  THE  25mm  WEAPON 


A.  Weapon  Description 

In  the  previous  chapter,  the  analysis  of  gun  pointing  errors  was 
presented  for  the  20mm  XM197  weapon.  As  a continuation  of  this  study, 
the  helicopter  model  and  the  gun-turret  model  were  used  to  evaluate  the 
dynamic  responses  of  th<"  larger  caliber  weapons,  including  the  25mm  and 
30mm  automatic  cannons.  This  section  of  the  report  deals  with  the  25mm 
gun;  the  30mm  gun  analysis  is  described  in  the  following  section. 

The  weapon  used  in  this  study  is  the  25mm  Bushmaster  cannon  that  is 
an  automatic,  single-barrel  weapon  with  the  rate  of  fire  in  the  vicinity 
of  500  spm.  This  gun  has  a round  impulse  of  about  57  lb-sec  which  can  be 
compared  with  33.15  lb-sec  for  the  20mm  XM197  weapon. 

B.  Recoil  Force 

Two  different  recoil  forces  were  used  in  this  study.  The  first  of 
these  is  a standard  recoil  force  and  the  second  is  a constant  recoil  force. 
Since  no  experimentally  measured  recoil  data  were  available  for  the  25mm 
gun,  the  necessary  data  were  generated  by  use  of  the  measured  recoil  force 
for  the  XM197  20mm  weapon,  which  is  described  as  recoil  force  No.  3 in 
Table  1.  The  original  data  were  generated  for  a firing  rate  of  669  spm, 
but  these  data  were  conversed  to  the  500  spm  firing  rate,  which  was  assumed 
to  be  representative  of  the  25mm  gun.  This  conversion  was  made  by  an  in- 
crease in  the  time  scale  of  the  weapon  recoil  data  by  a factor  of  669/500. 

In  addition  to  the  adjustment  for  the  different  firing  rates,  an  adjustment 
was  made  on  the  amplitude  of  the  recoil  force.  If  the  assumption  is  made 
that  the  recoil  traces  for  the  20mm  and  25iam  weapons  are  similar  except 
for  the  amplitude,  then  the  maximum  force  of  each  weapon  will  be  propor- 
tional to  the  impulse.  On  the  basis  of  a round  impulse  of  33.15  lb-sec  for 
the  XM197  gun  and  57  lb-sec  for  the  Bushmaster  25mm  gun,  c'ne  correction  for 
the  maximum  force  was  obtained  by  the  application  of  the  factor  57/33.15  to 
the  recoil  force.  This  constant  was  included  in  the  computer  program  of  the 
AH-1G  Helicopter  Vibration  Model,  as  a multiplying  factor  to  the  vector 
B(I)  in  subroutine  REDREC. 

The  25mm  constant  recoil  force  data  were  generated  from  the  Honeywell 
data  with  the  code  name  DDFFTA,  which  was  again  originally  obtained  for  the 
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20mm  weapon.  Since  these  data  correspond  to  a 750  spm  firing  rate,  the 
time  scale  was  modified  by  a factor  of  750/500.  The  modification  of  the 
peak  force  was  the  same  as  that  used  for  the  standard  recoil  force;  that 
is,  the  force  was  increased  by  a factor  of  57/33.15. 

C.  Turret  Controls 

Since  the  characteristics  of  the  turret  for  the  25mm  weapon  are  not 
as  yet  determined,  the  decision  was  made  to  use  the  XM97  turret  with  a few 
design  modifications.  These  modifications  were  performed,  in  order  to  im- 
prove the  dynamic  response  of  the  turret  drive  systems  to  the  transmitted 
recoil  force  of  the  higher  impulse  weapon.  They  include  the  increasing  of 
the  motor  torque  response  from  0.10  to  0.13  lb-ft  per  volt  for  both  the 
azimuth  and  the  elevation  drives,  the  increasing  of  the  motor  saturation 
voltage  from  26  to  35  volts,  and  the  increasing  of  the  transmission  stiff- 
ness from  5 x 10^  to  10  x 10^  lb-ft /rad. 

D.  Mass  and  Inertia  Properties 

The  mass  of  the  modified  XM97  turret  was  taken  to  be  220  pounds  and 
that  of  the  gun  to  be  250  pounds.  The  position  of  the  center  of  gravity 
of  the  gun  relative  to  the  elevation  trunnion  axis  plays  an  important  role 
in  the  dynamic  analysis  because  a large  mass  unbalance  is  produced  by  the 
location  of  the  gun  center  of  gravity  at  a large  distance  from  the  pivot. 
Since  the  final  gun-mounting  position,  relative  to  trunnion  location,  is 
not  firm  at  this  stage,  four  different  weapon-mounting  positions  were  ex- 
amined. These  four  posicons  are  listed  in  Table  14,  in  which  the  parameters 
and  are  the  distances  that  were  defined  previously  in  Figure  4.  The 
trunnion  elevation  axis  was  assumed  to  intersect  the  turret  azimuth  axis 
and  to  be  perpendicular  to  the  azimuth  axis.  The  center  of  gravity  of  the 
combined  gun  and  turret,  denoted  by  the  parameter  A^,  will  vary  depending 
on  the  position  of  the  gun  from  the  elevation  trunnion,  represented  by  the 
parameter  A^.  Table  14  contains  a listing  of  the  parameters  A^  and  A^  for 
each  of  the  selected  gun  positions.  For  the  25mm  turret  configuration,  the 
distances  B^,  and  B ^ were  assumed  to  be  zero. 

The  maximum  value  of  A^  was  chosen  to  be  2.67  feet  (32  inches)  and  the 
minimum  value  was  taken  to  be  zero,  in  which  case  the  gun  center  of  gravity 
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TABLE  14 

Position  of  Gun  Center  of  Gravity 


Gun  Center 
of  Grevlty 
Position 

Distance 
At  (ft) 

for  the  25mm 
Distance 

*3  (ft) 

Analysis 

Azimuth  Moment  of 
Inertia 
(slug-ft*) 

Elevation  Moment  of 
Inertia 
(slug-ft2) 

1 

1.2 

2.67 

65.0 

63.0 

2 

0.89 

2.0 

41.0 

39.0 

3 

0.59 

1.33 

23.7 

21.7 

4 

0.0 

0.0 

10.0 

8.0 

would  be  directly  on  the  turret  elevation  axis.  The  distance  was 

chosen  to  be  about  45  percent  of  A^.  This  choice  is  somewhat  arbitrary, 

but  should  be  quite  close  to  the  actual  value.  The  largest  effect  of  these 

distances  is  on  the  rotational  moments  of  inertia.  These  rotational  moments 

of  inertia,  both  for  the  azimuth  and  for  the  elevation  motions,  are  given 

in  Table  14.  Both  of  these  inertias  are  expressed  in  relation  to  the  turret 

and  gun  pivot  points  and  were  calculated  on  the  assumption  that  the  turret 

2 

moment  of  inertia  about  its  azimuth  axis  is  2.0  slug-ft  and  that  the  gun- 

saddle  moment  of  inertia  about  any  axis  normal  to  the  barrel  is  approximately 
2 

8.0  slug-ft  . Consequently,  the  total  inertia  for  the  azimuth  motion  was 
calculated  from  the  relation 

I = 2.0  + 8.0  + 250(A«)2/32.2 
a j 

where  the  last  term  represents  the  gun-mass  unbalance  and  is  proportional 
to  the  square  of  the  distance  A^.  The  corresponding  expression  for  the 
elevation  motion  was  calculated  from 

I = 8.0  + 250(A_)2/32.2. 
e 3 

E.  Results 

The  25mm  analysis  was  performed  with  the  AH-lG  version  of  th«  Helicop- 
ter Vibration  Model  and  with  the  modified  Gun  Turret  Model.  Calculations 
were  performed  for  the  feur  different  center  of  gravity  positions  given 
in  Table  14,  and  for  the  two  different  recoil  forces  described  previously. 
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The  results  of  these  calculations  are  the  statistical  means  and  the 
standard  deviations  of  the  dynamic  errors  in  the  azimuth  and  elevation 
gun-position  angles.  For  each  gun-denter-of-gravity  position  in  Table  14, 
calculations  were  performed  for  twelve  different  angular  positions  of  the 
weapon.  These  positions  involve  three  different  elevation  angles  of  -5°, 
-25°,  and  -45°;  and  four  azimuth  angles  of  0°,  22.5°,  45°,  and  67.5°. 

The  results  of  these  calculations  are  given  in  Tables  15  to  22.  Tables 
15  to  18  contain  the  results  for  the  standard  reccil  force,  and  Tables  19 
to  22  contain  the  results  for  the  constant  recoil  force.  Each  table  cor- 
responds to  a different  position  of  the  gun  center  of  gravity.  The  statis- 
tical errors  are  in  units  of  milliradians.  These  results  reveal  that  the 
maximum  values  of  the  standard  deviation  parameters  were  obtained  for  the 
gun-center-of-gravity  position  No.  1,  which  corresponds  to  the  largest  mass 
unbalance. 

The  maximum  standard  deviation  occurs  in  the  elevation  motion  and 
corresponds  to  the  largest  azimuth  angle.  In  that  respect,  the  helicopter 
response  to  the  25mm  standard  recoil  force  is  similar  to  that  of  the  20mm 

weapon.  Some  quantitative  comparisons  of  the  results  for  the  20mm  and  the 

tabu  is 

Statistical  Error  ?ariMt«ri  for  tha 
Standard  25am  focza  and  Gun  CG  Position  Wo.  1 


Tur 

Orlant 

E location 

■AlO  « 

Hnsi 

jg| 

-5* 

0! 

.56 

— 

1.26 

4.28 

Hi 

.28 

1.08 

0.25 

4.98 

p BP 

.37 

1.43 

-2.90 

7.74 

67.5* 

.36 

1.10 

-4.90 

11.29 

-M* 

H 

.29 

.27. 

.81 

1.84 

mm 

-.65 

.94 

.04 

2.91 

45.V 

-i.38 

2.41 

-1.76 

11.09 

o 

o 

1.80 

-3.62 

19.34 

-«• 

HI 

.13 

.42 

3.16 

3.84 

HI 

-1.99 

1.64 

2.03 

4.11 

45.0* 

-2.B9 

2.50 

-.83 

8.24 

67.5* 

-2.29 

L12 

-2-37 

14.36 
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25mm  weapon  systems  are  now  given.  The  results  were  obtained  for  the  20mm 
XM197  weapon  mounted  on  the  unmodified  XM97  turret  and  for  the  25mm  Bush- 

master  weapon  mounted  on  the  modified  XM97  turret.  The  firing  rate  was 
500  spm  for  both  weapons,  but  the  two  weapon  systems  used  in  the  com- 
parisons differed  in  several  respects  besides  the  weapon  caliber.  In 
both  analyses,  the  AH-lG  version  of  the  helicopter  dynamic  model  was  used. 

Consider  first  the  azimuth  and  the  elevation  angles  of  67.5°  and 
-5°,  respectively.  For  this  orientation  of  the  gun,  the  standard  de- 
viation of  the  elevation  vibration  was  found,  in  Table  10  for  tiif  20mm 
weapon,  to  be  13.25  milliradians  and  the  corresponding  value  for  the  25mm 
weapon,  given  in  Table  15,  is  11.29  milliradians.  This  value  corresponds 
to  the  largest  mass  unbalance,  whereas  the  corresponding  value  for  the 
smallest  unbalance,  given  in  Table  18,  is  5.17  milliradians.  Consider  now 
the  gun  orientation  of  67.5°  and  -45°,  for  the  azimuth  and  for  the  elevation, 
respectively.  At  that  orientation,  the  20mm  analysis  gave  the  result  for 
the  elevation  standard  deviation  of  7.77  milliradians;  the  corresponding 
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values  for  the  25mm  analysis  vary  from  14.36  to  5.00  milliradians . The 
largest  standard  deviation  found  in  these  studies  is  19.3  milliradians 
in  the  elevation  motion  and  it  corresponds  to  an  initial  angular  orien- 
tation of  67.5°  in  azimuth  and  -25°  in  elevation.  This  value  was  obtained 
for  the  largest  mass  unbalance  from  Table  15.  From  these  results,  the 
conclusion  is  made  that,  if  the  mass  unbalance  of  the  weapon  is  large,  the 
25mm  weapon  standard  recoil  force  can  lead  to  gun  vibration  errors  of 
approximately  twice  those  obtained  for  the  20irm  weapon. 

The  significance  of  the  gun-mass-unbalance  parameter,  A^,  is  il- 
lustrated in  Figure  25.  This  figure  contains  plots  of  the  standard  de- 
viation error  for  the  elevation  motion  versus  the  turret  azimuth  orien- 
tation. The  plots  correspond  to  three  separate  values  of  the  gun  position 
parameter  and  the  initial  turret  elevation  orientation  was  -25°  in 
each  case.  As  expected,  the  smallest  standard  deviation  error  is  genera- 
ted when  the  parameter  is  zero.  For  larger  values  of  the  parameter  A^ , 
the  standard  deviation  error  in  elevation  usually  increases,  although  this 
does  not  always  happen.  Usually,  however,  the  dynamic  response  of  the 
turret  would  be  improved  by  the  reduction  of  the  distance  A^  to  as  small 
a value  as  is  possible. 

Consider  now  the  comparison  of  the  response  to  the  standard  recoil 
force  with  the  response  to  the  constant  recoil  force.  The  effect  of  the 
constant  recoil  mechanism  can  be  seen  by  a comparison  of  the  results  in 
Tables  15  to  18  with  the  corresponding  results  in  Tables  19  to  22.  For 
each  position  of  the  gun  center  of  gravity,  the  constant  recoil  force 
yields  appreciable  smaller  values  of  the  standard  deviations.  This  is 
especially  apparent  for  the  standard  deviation  error  in  the  elevation  an- 
gle. The  constant  recoil  force  yields  a reduction  in  this  error  that 
ranges  from  about  30  to  60  percent  of  the  corresponding  result  for  the 
standard  recoil  force.  A constant  recoil  mechanism  would,  therefore,  bring 
about  a significant  reduction  in  the  dynamic  gun-pointing  error. 


V.  ANALYSIS  OF  THE  30mm  WEAPON 


A.  Weapon  Description 

The  dynamic  response  of  the  AH-1G  helicopter  and  the  gun-turret 
system  was  next  examined  with  the  application  of  the  weapon  recoil  for- 
cing function  for  the  AMCAWS  30mm  weapon.  This  weapon  is  a single- 
barrel  automatic  cannon  and  has  an  approximate  firing  rate  of  450  spm 
and  a very  large  round  impulse  of  approximately  150  lb-sec. 

In  the  performance  of  this  analysis,  the  decision  was  made  to  use 
only  the  constant  recoil  force  because  the  magnitude  of  the  peak  recoil 
of  this  weapon  is  very  large.  A peak  recoil  force  of  about  15,000  pounds 
was  measured  for  the  AMCAWS  30ram  weapon  on  a hard  stand.  Consequently, 
on  the  basis  of  the  results  for  the  20mm  gun,  extremely  large  motions  of 
the  helicopter  and  the  gun-turret  configuration  would  probably  result 
without  the  use  of  the  constant  recoil  mechanism.  In  addition,  the  AH- 
1G  helicopter  is  not  designed  to  withstand  peak  recoil  forces  of  this 
magnitude.  Such  large  motions  could  not  be  tolerated  by  the  system. 

Since  no  turret  has  yet  been  designed  to  house  the  AMCAWS  30mm  gun, 
the  present  analysis  makes  use  of  a modified  XM97  turret.  The  computer 
program  for  the  analysis  of  the  original  XM97  turret  has  been  modified 
and  used  in  the  present  investigation. 

B.  Recoil  Force 

Two  different  types  of  weapon  recoil  data  were  used  to  provide  the 
dynamic  excitation  force.  The  first  of  these  is  a simulated  constant  re- 
coil force  with  no  misfires  and  the  second  is  a simulated  constant  recoil 
force  with  two  misfired  rounds. 

No  experimental  data  are  available  at  the  present  time  for  the  recoil 
history  of  the  AMCAWS  30mm  weapon  and,  consequently,  the  appropriate  data 
for  this  analysis  were  generated  from  existing  constant  recoil  data  sup- 
plied by  the  Honeywell,  Inc.,  for  the  20mm  weapon.  The  two  data  sets 
tha:  were  used  are  identified  by  codes  RAANM  and  RHHM2  and  are  described 
in  Table  2.  To  use  these  data  for  the  AMCAWS  30mm  gun  required  the  modi- 
fication of  both  the  round  impulse  and  the  firing  rate.  The  projected 
firing  rate  of  the  AMCAWS  30mm  gun  is  450  shots  per  minute  and,  therefore, 
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the  time  scale  of  the  RAANM  data  was  increased  by  a factor  of  750/450 
and  the  time  scale  of  the  RHHM2  data  was  increased  by  820/450.  To 
account  for  the  higher  impulse  of  the  AMCAWS  30mm  weapon  an  increase 
was  made  to  the  constant  recoil  amplitude  of  the  20mm  recoi’  force  by 
a factor  of  3. 

C.  Turret  Controls 

Some  design  modifications  were  performed  on  the  XM97  turret  servo- 
control  system  so  that  the  turret  would  respond  adequately  to  the  higher 
impulse  loading  of  the  AMCAWS  30mm  weapon.  The  first  modification  in- 
volved a change  in  the  servocontrol  drive  motor  torque  response,  which 
was  increased  from  0.1  ft-lb  per  volt  to  0.13  ft- lb  per  volt.  The  second 
modification  concerned  the  motor  saturation  voltage  that  was  raised  from 

26  to  35  volts.  The  third  change  involved  the  motor  inertia  that  was 

_4  _4  2 

increased  from  2.4  x 10  to  4.0  x 10  slug-ft  . The  final  change  was 
an  increase  in  the  stiffness  of  the  gear  transmission.  The  original  stiff- 
ness was  5 x 10^  lb-ft/rad  and  it  was  increased  to  10  x 10^  and  15  x 10^ 
lb-ft/rad.  doth  of  these  larger  values  of  the  transmission  stiffness  were 
used  in  tne  analysis. 

D.  Mass  and  Inertia  Properties 

before  the  dynamic  analysis  of  the  turret  could  be  performed,  some 

estimates  had  to  be  made  as  regards  the  mass  and  the  moments  of  inertia 

of  the  AMCAWS  30mm  gun  and  turret.  The  total  mass  of  the  gun  and  the 

turret  was  estimated  at  400  pounds,  and  the  moments  of  inertia  about  both 

2 

the  azimuth  and  the  elevation  axes  were  chosen  as  50  slug-ft  . Further- 
more, the  gun  and  the  turret  centers  of  gravity  were  assumed  to  coincide, 
and  the  position  of  the  gun  center  of  gravity  was  chosen  to  be  at  the 
gun  elevation  pivot  point.  The  pivot  was  chosen  to  be  on  the  turret  azi- 
muth axis  and  29.5  inches  below  the  turret  attachment  point. 

These  modifications  to  the  turret  geometric  parameters  were  performed, 
ir.  spite  of  practical  considerations,  to  minimize  the  vibration  amplitude 
of  the  turret  drive  mechanisms.  Although,  in  this  regard,  the  modified 
30mm  turret  is  an  idealized  design  concept,  the  results  of  the  AMCAWS  30mm 
gun-pointing  error  analysis  should  provide  an  indication  of  the  practicality 
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of  mounting  a large  impulse  weapon  on  a helicopter. 

E.  Results 

Two  sets  of  numerical  results  were  obtained  corresponding  to  the  two 
different  values  of  the  transmission  stiffness,  k.  For  each  case,  the 
results  are  given  by  the  mean  value  and  by  the  standard  deviation  of  the 
dynamic  motion  in  both  the  azimuth  and  the  elevation  directions.  These 
results  are  obtained  for  different  orientations  of  the  gun.  The  orienta- 
tions of  the  gun  correspond  to  two  elevation  angles  of  -5°  and  -45°, 
and  four  azimuth  angles  of  0°,  22.5°,  45°,  and  67.5°.  These  results  are 
summarized  in  Tables  23  and  24.  The  results  in  Table  23  correspond  to 
the  transmission  stiffness  k = 10  x 10^  lb-ft/rad;  the  results  in  Table  24 
are  for  k * 15  x 10^  lb-ft/rad.  Each  table  contains  results  for  the  constant 
recoil  forces  RAANM  and  RHKM2.  As  expected,  the  recoil  force  with  misfires 
produces  the  largest  standard  deviation,  which  occurs  in  the  elevation  mode 
for  large  azimuth  angles.  These  results  are  similar  to  the  results  for  the 
20mm  gun,  but  are  larger  in  magnitude.  A comparison  of  the  results  in 
Tabxes  23  and  24  shows  that  a large  increase  in  the  transmission  stiffness 
appears  to  cause  a slight  decrease  in  the  dynamic  response.  In  most  cases, 
the  statistical  errors  for  the  idealized  AMCAWS  30mm  constant  recoil  force 
(code  RAANM  modified)  are  less  than  2 milliradians.  The  corresponding 
results  for  the  iaealized  constant  recoil  force  with  misfires  (code  RHHM2 
modified)  are  much  larger;  however,  the  statistical  errors  did  not  exceed 
10  milliradians,  for  any  turret  orientation.  Before  a complete  assessment 
of  the  dynamic  errors  can  be  made  for  the  AMCAWS  30mm  weapon,  this  analysis 
should  be  repeated  with  experimental  data  measured  from  an  actual  constant 
recoil  mechanism.  To  date,  no  constant  recoil  mount  design  has  been  de- 
veloped for  the  AMCAWS  30mm  weapon.  Such  a device  might,  however,  be  ex- 
pected to  yield  results  which  are  slightly  larger  than  the  results  that 
are  presented  in  this  section  for  the  idealized  constant  recoil  force. 
Nevertheless,  from  the  results  of  this  analysis  the  feasibility  of  mount- 
ing the  AMCAWS  30mm  gun  on  an  attack  helicooter  has  been  demonstrated. 
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VI.  EXTRAPOLATION  MODEL 

During  the  present  investigation,  large  amounts  of  data  were  obtained 
for  the  statistical  gun-pointing  errors  from  various  helicopter,  gun, 
and  turret  configurations.  For  each  type  of  helicopter  model,  recoil 
force,  and  weapon  configuration,  the  statistical  errors  were  evaluated 
for  specific  azimuth  and  elevation  angles  of  the  gun.  However,  in  an 
analysis  with  the  complete  Helicopter  Weapon  System  Model,  shown  schemat- 
ically in  Figure  3,  gun  orientations  different  from  those  that  have  been 
examined  may  frequently  arise.  In  this  situation,  the  helicopter  vibration 
analysis  and  the  turret  dynamic  analysis  could  be  repeated  for  the  desired 
orientation.  However,  this  approach  is  not  practical  since  these  two 
models  are  relatively  complex  and  require  an  appreciable  amount  of  com- 
puter storage  and  execution  time.  To  remedy  this  difficulty,  an  extra- 
polation model  was  developed. 

With  this  approach,  the  statistical  mean  and  standard  deviation  quan- 
tities for  any  desired  angular  gun  orientation  are  extrapolated  from  a 
table  of  results,  obtained  fcr  a specific  set  of  initial  turret  orientations. 
The  model  incorporates  a linear  extrapolation  procedure  to  calculate  the 
desired  statistical  quantities.  One  requirement  of  this  approach  is  that 
the  desired  orientation  of  the  gun  shouli  be  within,  or  close  to,  the  range 
of  the  angular  orientations  that  have  previously  been  analyzed.  The  com- 
puter program  is  relatively  simple  and  requires  little  storage  or  execution 
time.  The  input  for  the  extrapolation  model  computer  program  consists  of 
matrices  containing  the  statistical  results  for  a particular  type  of  recoil 
force.  Each  element  of  this  matrix  corresponds  to  one  combination  of  the 
azimuth  and  elevation  angles.  The  computer  program  allows  for  the  storage 
of  more  than  one  set  of  data,  and  each  set  of  data  corresponds  to  a par- 
ticular recoil  force. 
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VII.  SUMMARY  AND  CONCLUSIONS 

The  gun-pointing  errors  induced  by  weapon  recoil  forces,  in  a heli- 
copter weapon  system,  were  analyzed  for  several  weapons  of  different  cali- 
ber. In  these  investigations,  the  weapons  were  mounted  on  a biaxial 
turret  and  the  gun-pointing  accuracies  were  determined  for  various  tur- 
ret configurations.  Three  distinct  models  were  used  to  carry  out  the  ana- 
lysis. The  first  model  - a NASTRAN  finite-element  model  of  the  helicop- 
ter structural  airframe  - was  used  to  solve  the  free-vibratxon  problem, 
and  the  primary  outputs  of  this  model  are  the  natural  frequencies  and  mode 
shapes  of  the  helicopter  airframe.  These  data  are  required  by  the  Heli- 
copter Vibration  Model  to  solve  the  forced-vibration  problem  for  weapons 
recoil  force  excitations.  In  particular,  the  Helicopter  Vibration  Model 
is  used  to  predict  the  displacement,  velocity,  and  acceleration  components 
of  the  turret  attachnant  platform.  These  data  are  then  fed  into  the  Gun- 
Turret  Dynamic  model,  and  this  model  is  used  to  determine  the  dynamics  of 
the  turret  drive  systems  in  the  azimuth  and  in  the  elevation  directions. 

As  a final  step  in  the  analysis,  the  statistical  mean  and  the  standard  de- 
viation parameters  are  estimated  from  the  vibrations  of  the  azimuth  and 
the  elevation  drive  mechanisms.  These  statistical  parameters  provide  a 
convenient  measure  of  the  angular  errors  in  gun-pointing. 

A structural  sensitivity  analysis  is  carried  out,  with  the  use  of 
these  models,  by  the  evaluation  of  the  statistical  gun-pointing  error  com- 
ponents for  several  discrete  initial  orientations  of  the  turret.  Normally 
cne  analysis  is  carried  out  for  a selection  of  twelve  discrete  turret 
orientations,  which  correspond  to  a wide  range  of  azimuth  and  elevation 
angles.  The  results  of  this  analysis  reveal  the  sensitivity  of  each  gun- 
pointing error  component  to  the  initial  turret  orientation.  This  pro- 
cedure has  been  applied  to  the  investigation  of  the  20mm  XM197,  the  25mm 
Bushmaster,  and  the  30ram  AMCAWS  weapon  systems.  In  each  investigation, 
the  effects  of  different  turret  and  mount  configurations  on  the  gun-pointing 
error  components  vere  determined. 

For  the  20mm  XM197  weapon,  the  weapon  recoil  data  were  obtained  from 
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both  experimental  tests  and  computer  simulations.  These  data  consisted 
of  standard  recoil  data,  measured  from  recoil  adapter  mounts,  and  constant 
recoil  data,  obtained  from  a Honeywell  constant  recoil  mount  mechanism. 

Data  were  recorded  on  each  type  of  mount  for  several  weapon  firing  rates. 

The  basic  weapon  system  configuration,  which  was  used  for  the  analysis 
of  the  20mm  automatic  cannon,  consisted  of  the  XM197  gun,  the  XM97  turret, 
and  the  Bell  AH-1  helicopter  airframe.  Two  different  NASTFAN  models  were 
used  to  represent  the  free-vibration  characteristics  of  the  model  AH-lJ 
and  the  model  AH-1G  airframes.  The  models  differed  primarily  in  the  detail 
of  the  structural  representation.  For  example,  the  AH-1G  model  was  refined 
and  contained  about  1800  independent  degrees  of  freedom,  whereas  the  AH- 
1J  model  was  much  simpler  and  contained  only  about  300  degrees  of  freedom. 

The  results  of  the  structural  sensitivity  analysis  of  the  20mm  XM197 
weapon  are  presented  in  Chapter  III.  Analytical  predictions  of  the  gun- 
pointing  error  components  were  obtained  for  both  the  AH-lJ  and  the  AK-1C 
versions  of  the  structural  model.  In  general,  the  error  nredictions  for 
the  two  models  were  of  the  same  order  of  magnitude.  However,  the  results 
of  the  AH-lJ  model  revealed  a greater  sensitivity  to  the  turret  azimuth 
angle  than  the  corresponding  results  of  the  AH-1G  model.  The  discrepancies 
in  the  model  predictions  are  caused  by  the  differences  in  the  AH-lJ  and 
tne  AH-1G  model  representations.  Since  the  AH-1G  model  contains  a more 
refinad  representation  of  the  helicopter  structure,  the  analytical  results 
for  the  AH-iG  model  are  probably  more  accurate  than  the  corresponding  re- 
sults of  the  AH-lJ  model. 

In  almost  all  cases  analyzed,  the  largest  statistical  error  comnonent 
was  the  standard  deviation  of  the  elevation  motion.  In  general,  the  value 
obtained  for  this  error  was  dependent  upon  the  type  of  weapon  mount,  the 
weapon  firing  rate,  and  the  turret  orientation.  For  the  analysis  of  the 
AH-IG  helicopter  weapon  system  with  an  applied  20mm  standard  recoil  force, 
this  error  varied  between  about  1.2  milliradians  and  15.5  milliradians . 

In  this  analysis,  the  results  showed  that  the  gun-pointing  error  is  very 
sensitive  to  the  weapon  firing  rate  and  tc  the  turret  orientation.  The 
largest  errors  usually  occurred  at  the  largest  azimuth  orientation  of  the 
turret.  At  this  orientation,  large  torques  are  applied  to  the  helicopter 
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structure,  about  the  longitudinal  axis. 

The  magnitudes  of  the  error  components  obtained  for  the  20mm  constant 
recoil  forces  were  significantly  lower  than  the  magnitudes  obtained  for 
the  corresponding  analysis  of  the  standard  recoil  forces.  In  particular, 
the  statistical  error  components  were  almost  always  less  than  3 milli- 
radians;  and,  in  almost  all  instances,  the  errors  were  less  than  2 milli- 
radians.  These  results,  therefore,  show  that  the  constant  recoil  mount 
is  very  effective  in  the  reduction  of  the  gun-pointing  errors.  The  amount 
of  improvement  obtained  from  the  constant  recoil  mechanisms  was  dependent 
on  the  initial  turret  orientation  and  on  the  weapon  firing  rate.  However, 
the  constant  recoil  device  usually  resulted  in  the  reduction  of  the  gun- 
pointing errors  by  a factor  of  3 or  more.  In  addition,  the  gun-pointing 
errors  for  the  20mm  constant  recoil  force  excitations  were  somewhat  in- 
sensitive to  the  firing  rate,  whereas  the  results  for  the  standard  recoil 
force  excitation  showed  a much  greater  sensitivity  to  the  firing  rate. 

Other  studies  were  conducted  for  the  20mm  XM197  weapon  system  to  de- 
termine whether  the  XM197  turret  servocontrol  response  could  be  improved 
by  the  incorporation  of  an  additional  feedback  in  the  servocontrol  loop. 
This  additional  feedback  was  expressed  in  the  form  of  an  angular  rate  sig- 
nal or  an  acceleration  signal.  This  type  of  feedback  could  be  measured  by 
a rate  gyro  or  an  accelerometer  mounted  on  the  gun  or  turret.  None  of 
the  proposed  turret  servocontrol  modifications  affected  a noticeable  im- 
provement in  the  dynamic  response  of  the  gun.  Possibly,  a better  method 
for  improving  the  turret  servocontrol  system  would  be  to  redesign  the 
system  on  the  basis  of  optimal  control  techniques  rather  than  to  add  a 
feedback  to  the  classical  control  system.  Plans  are,  therefore,  being 
made  for  an  investigation  of  the  dynamic  response  of  an  optimally  con- 
trolled turret. 

After  the  completion  of  the  20mm  weapon  system  analysis,  an  investi- 
gation was  carried  out  for  the  .Larger  caliber  Bushmaster  25mm  and  the 
AMCAWS  30mm  weapons  and  the  results  of  the  investigation  are  presented 
in  Chapters  IV.  and  V.  The  analysis  of  the  25mm  gun  was  conducted  for 
both  the  standard  recoil  adapter  mount  and  the  constant  recoil  mount. 
However,  the  analysis  of  the  30mm  weapon  was  limited  to  the  constant  recoil 
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mount  only.  Since  no  experimental  recoil  data  were  available  for  either 
of  these  weapons,  the  recoil  traces  for  the  25mm  and  for  the  30mm  guns 
were  estimated  from  the  weapon  design  specifications  and  from  the  assump- 
tion that  the  25mm  and  30mm  recoil  traces  are  similar  to  the  recoil  traces 
for  the  20mm  weapon.  Because  the  larger  caliber  weapons  are  more  massive 
and  larger  in  size  than  the  20mm  XM197  gun,  certain  improvements  had  to 
be  made  to  the  physical  design  characteristics  of  the  XM97  turret.  Some 
of  these  improvements  involved  the  relocation  of  the  turret  trunnion  axis 
and  the  relocation  of  the  center  of  gravity  positions  of  the  gun  and  the 
major  turret  components.  As  a result  of  these  improvements,  the  mass 
unbalances  in  the  turret  were  reduced  considerably  and,  hence,  the  vi- 
brations of  the  turret  azimuth  and  elevation  drive  systems  were  also  re- 
duced. 

After  these  changes  had  been  made  to  the  turret  design,  the  struc- 
tural sensitivity  analysis  was  carried  out  for  the  Bushmaster  25mm  weapon 
and  for  the  AMCAWS  30mm  weapon.  In  the  analysis  of  the  25mm  weapon,  dif- 
ferent values  were  selected  for  the  parameter  A„  by  which  the  mass  un- 
balance of  the  gun  is  defined.  As  expected,  the  gun-pointing  errors  were 
generally  larger  for  large  mass  unbalances  than  for  small  mass  unbalances. 
The  largest  mass  unbalance  corresponded  to  a turret  mounting  configuration 
which  was  similar  to  that  used  for  the  20mm  XM197  weapon.  In  this  case, 
the  results  for  the  25mm  analysis  were  about  twice  as  large  as  the  cor- 
responding results  for  the  20mm  standard  recoil  analysis.  For  the  case 
in  which  the  mass  unbalance  parameter  was  equal  to  zero,  the  results  for 
the  25mm  standard  recoil  analysis  were  not  excessively  large.  The  mag- 
nitudes of  the  statistical  gun-pointing  errors  did  not  exceed  8.26  milli- 
radians  and,  in  most  cases,  these  errors  were  much  smaller.  An  even 
greater  improvement  in  the  gun-pointing  errors  occurred  when  the  gun  was 
mounted  on  a constant  recoil  mechanism.  In  this  case,  the  errors  were 
reduced  by  a factor  of  2 or  more. 

In  the  analysis  of  the  AMCAWS  30mm  weapon,  the  gun  mass  unbalance 
parameter  was  set  equal  to  zero,  and  different  values  were  selected  for 
the  transmission  stiffness.  The  choice  of  the  transmission  stiffness, 
however,  had  little  effect  on  the  numerical  results.  Almost  all  predicted 
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gun-pointing  errors  were  small  for  the  case  in  which  an  idealized  constant 
recoil  force  was  applied  to  the  modified  turret.  In  fact,  the  maximum 
errors  occurred  for  the  67.5°  azimuth  angle  position  and  these  errors  did 
not  exceed  3.2  milliradians.  The  idealized  30mm  constant  recoil  force 
that  contained  misfires,  however,  caused  larger  errors  to  occur  in  the 
elevation  mechanism,  especially  at  large  azimuth  angle  positions. 

From  the  AMCAWS  30mm  gun-pointing  error  analysis,  the  feasibility  of 
the  use  of  a high-performance  weapon  in  aircraft  applications  has  been 
demonstrated,  with  the  provision  that  the  gun  is  to  be  attached  to  a con- 
stand  recoil  mount.  These  results,  however,  are  somewhat  tentative  in 
that  no  actual  weapon  recoil  test  data  were  available  to  perform  the  ana- 
lysis. A continuation  of  the  structural  sensitivity  analysis  is  recommend- 
ed for  the  Bushmaster  25mm  weapon  and  for  the  AMCAWS  30mm  weapon,  in  which 
actual  experimental  recoil  force  data  are  used  in  place  of  the  computer 
simulated  data.  Before  this  analysis  can  be  undertaken,  a prototype  design 
of  a constant  recoil  mount  for  these  weapons  would  have  to  be  built  and 
tested.  Recently,  a theoretical  study  has  been  initiated  on  the  concept 
development  of  a constant  recoil  mount  for  the  AMCAWS  30mm  gun;  however, 
no  hardware  development  program  has  been  initiated  yet. 

Other  rypes  of  mount  designs  have  been  developed  for  various  air- 
craft weapons  that  may  also  improve  the  dynamic  responses  of  the  turret 
drive  systems.  For  example,  General  Electric  Company  has  developed  a new 
spring  mount  for  the  30mm  XM188  weapon  that  transmits  smaller  peak  recoil 
forces  than  the  standard  spring  adapter  mount.  The  force  transmitted  by 
the  improved  spring  mount  may  not  be  as  smooth  as  the  force  transmitted 
by  a constant  recoil  mount;  nevertheless,  the  improved  spring  mount  is 
more  compact  and  less  expensive  than  the  constant  recoil  mount.  An  ana- 
lysis is,  therefore,  recommended  to  quantify  the  amount  of  improvement  in 
gun-pointing  accuracy  that  could  be  achieved  by  the  replacement  of  the 

f 

standard  spring  adapter  mount  with  the  improved  mount. 

Besides  the  lack  of  suitable  recoil  data,  a problem  that  arose  in  the 
analysis  of  the  25mm  and  30mm  weapons  was  the  lack  of  a suitable  turret 
for  mounting  these  weapons.  A modified  version  of  the  XM97  turret  was, 
therefore,  used  in  the  25mm  and  30mm  investigations,  and,  in  most  instances, 
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the  resulting  gun-pointing  error  predictions  for  the  modified  turret  de- 
sign were  not  excessively  large.  A much  greater  improvement  in  the  dynamic 
response  of  a large  caliber  turret  would  probably  be  achieved  by  a com- 
plete design  optimization  of  the  turret.  An  extension  of  the  25mm  and 
30mm  recoil  sensitivity  investigations  is,  therefore,  recommended  that 
would  include  the  optimization  of  both  the  physical  design  characteristics 
and  the  servocontrol  characteristics  of  the  turret  for  precise  dynamic  gun 
pointing.  In  this  investigation,  the  weights,  center  of  gravity  locations, 
trunnion  locations,  and  other  important  design  parameters  of  the  turret  will 
be  chosen  to  minimize  the  dynamic  effects  of  weapon  recoil  forces.  Further- 
more, the  classical  servocontrol  system  will  be  replaced  by  a control  system 
design  developed  on  the  basis  of  an  optimal  control  technique.  Some  addi- 
tional turret  instrumentation  such  as  rate  gyros  or  accelerometers  may  be 
incorporated  in  the  optimal  control  system  design  concept. 

In  all  probability,  the  design  of  a high-performance  25mm  er  30mm 
weapon  system  would  be  intended  for  use  on  a next-generation  hell  copter 
such  as  the  Advanced  Attack  Helicopter  (AAH) . After  the  Army  has  made  a 
selection  between  the  Eell  and  the  Hughes  versions  of  the  AAH,  the  struc- 
tural sensitivity  analysis  should  be  repeated  for  the  selected  AAH.  The 
results  of  the  structural  sensitivity  investigations  for  the  AH-1G  heli- 
copter airframe  revealed  that  the  gun-pointing  error  is  very  sensitive  to 
the  natural  frequencies  of  vibration  and  to  the  mode  shapes  of  the  heli- 
copter structure.  Since  the  free-vibration  characteristics  of  the  AAH 
helicopter  might  be  quite  different  from  that  of  the  AF-lG  helicopter,  the 
results  of  a gun-pointing  error  analysis  for  the  AAH  helicopter  weapon 
system  may  differ  appreciably  from  the  AH-iG  results.  An  investigation 
of  the  gun-pointing  errors  is,  therefore,  recommended  for  the  AAH  weapon 
system. 

In  all  completed  analyses  of  the  various  factors  that  contribute  to 
the  total  dynamic  gun-pointing  error,  the  dynamic  gun  barrel  whin  had  not 
been  considered.  Often  the  flexural  vibration  amplitudes  of  the  gun  barrel 
are  less  than  one  milliradian  and,  therefore,  these  vibrations  can  be  ne- 
glected in  comparison  with  the  turret  and  helicopter  vibrations.  For  the 
larger  impulse  weapons  such  as  the  Bushmaster  25mm  and  the  A>*CAWS  30mm 
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weapons,  the  gun  barrels  are  rather  long  (85  inches  for  the  AMCAWS  30mm 
gun);  consequently,  the  flexural  vibrations  of  the  gun  barrels  of  these 
larger  caliber  weapons  may  be  a significant  source  of  gun-pointing  error. 

An  analysis  should,  therfore,  be  carried  out  for  these  weapons,  to  quan- 
tify the  magnitudes  of  the  dynamic  oscillations  of  the  gun  barrel  during 
a firing  event.  The  analytical  modeling  techniques  that  are  needed  to 
perform  the  barrel  vibration  calculations  have  already  been  developed, 
and  these  techniques  can  readily  be  applied  to  the  investigations  of  the 
Bushmaster  25mm  and  the  AMCAWS  30mm  weapons. 

In  conclusion,  the  structural  sensitivity  analysis  has  proven  to  be 
a valuable  tool  for  the  evaluation  and  quantification  of  the  dynamic  per- 
formance of  existing  or  conceptual  helioopter  weapon  systems.  Furthermore, 
the  dynamic  error  prediction  techniques  have  been  useful  for  the  identi- 
fication of  potential  problem  areas,  such  as  the  excitation  of  vibration 
resonance  at  certain  weapon  firing  rates,  and  these  techniques  have  also 
been  beneficial  for  the  evaluation  of  certain  concept  improvements,  such 
as  the  constant  recoil  mount. 
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APPENDIX 

A.  FORTRAN  IV  Source  Listing  of  the  Helicopter  Vibration  Model 

C 

C 

C HELICOPTER  VIRRATION  MOOEL:  MAIN  PROGRAM 

C 

C PURPOSE: 

C DETERMINES  THE  SOLUTION  TO  THE  FORCED  VIBRATION  PROBLEM  FOR  THE 
C HELICOPTER  STRUCTURAL  AIRFRAME.  THE  MODEL  GENERATES  THE  DISPLACEMENTS. 
C VELOCITIES.  AND  ACCELERATIONS  AT  THE  TURRET  ATTACHMENT  POINT. 

C 

C COMMON  BLOCKS  AND  THE  PROGRAMS  THAT  REQUIRE  THEM: 

C /ANGLE/  LOAD. MAIN 

C /DIST/  LOAD. MAIN 

C /MODE/  MRIN.OUPT.QMODE.QUAD 

C /RED/  MAIN  .OUTP  .REDREC 

C /WHY/  MAIN.  QMODE . QUAD 

C 

SUBROUTINES  CALLEO : 

DIME 
MODSEN 
C NOW 

C OUTP 

C QUAD 

C REDREC 

C RLOT 

C 

COMMON  /ANGLE/  ALPHA .PHI .SA .CR .SP .CP 
COMMON  /DIST/  B1.B2.D1.D2 

COMMON/MODE/A ( 180.30) .GM f 30 ) .GF( 30 ) .NGI 6 ) ,NM00 .NDF ,MI TER .WORK! 900 ) 
•.NiN.NOUT 

COMMON  /RED/  X0.H.CU003) 

COMMON  /WHY/  Y( 90 ) . IMOOEt 30 ) 

DIMENSION  IDISPC  30 ) , IVELt  30 ) .IRCCC30) 

DIMENSION  Y0UTC251  ) .TIME (251  ) 

DIMENSION  PRMT(S) 

DIMENSION  RPLTIS) .OONFIGC3) 

C 

C — INITIALIZE  TIMER 

C — TIMER  RESOLUTION  = 16440  MICROSECONDS 
C 

ITIMErC 

CALL  NOWC I TIME  ) 

C 

C READ  DISTANCES  B1  B2  B3  AND  D1  02  D3 
C 


A-l 


nnn  non  n n n non  non  non  non 


RERD  (5.1000)  Bl. 82. 01. 02 
WRITE  (6.2002)  81.B2.01.D2 

READ  CONSTRNT  RLPHR  AND  PHI 

RER0( 6 . 1000 )0RLPHR .0PHI 
PIr3. 1415927 
CPI=PI/180. 

RLPHR=CP I ■DRLPHA 
PHI=CPI«DPHI 
Sfl=SIN( RLPHfl ) 

CR=C0S( RLPHR ) 

SP=SIN( PHI ) 

CP=C0S( PHI ) 

RER0  NUMBER  OF  MODES  RND  DECREES  OF  FREEDOM 

RERD( 5.1001)  NM00.N0F 
IF(NM0D.GT .30)  NM00=30 
IF( NDF «0T • 180 ) N0F=180 
DO  40  1=1. NDF 
DO  40  Jrl.NMOD 
40  RC I .J)=0.0 

— RERD  THE  NUMBER  OF  MODAL  EQUATIONS  TO  BE  INTEGRATED 

REfiD  (5.1001)  NMODE 

RERD  (5.1001)  (IMODE(I). 1=1. NMODE) 

RERD  GLOBRL  POSITION  OF  TURRET  DISPLACEMENTS 

RERO(S.lOOl)  ( NOC I ) . 1=1 .6 ) 

READ  EIGENVALUES.  GENERALIZED  MASSES 

READ  (5.1008)  ( (OF(  I ) .QMU) ) . 1 = 1 .NMOO  ) 

N02=2«NM00E 

READ  INITIAL  DISPLACEMENTS  AND  VELOCITIES 

DO  60  1=1. N02 
60  Y(I)=0.0 

READ  LOWER  ANO  UPPER  BOUNDS  . I N I T I AL  INCREMENT .HMIN.HMRX 


I 


-j 

i 


I 

:! 
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READ  (5.1005)  PRMT 
WR I TE( 6 . 2003 ) 

00  62  1=1 .NM00E 

62  YOUT(I)  = SQRT(GF( IMODEI I) ) )/(2.«PI ) 

WR I TEC  6 .2004  ) (( IM00E( I) . YOUTH ))  .I=1.NM0DE) 

C 

C — DETERMINE  STEP  SIZE  FOR  INTEGRATION 
C 

h=PRMT ( 3 ) 

IF( I MODE ( NMODE ) .LE .6  ) GO  TO  63 
TRUMIN=2.0»3. 1415927/ SQRT(GF( IMODE( NM0OE ) ) ) 
H=0.05*TRUMIN 

63  SD=0 .02 

DELTR=PRMT( 2 )-PRMT( 1 ) 

HMIN  = PRMTC4) 

HMRX  = PRMTC5) 

IF(H.LT.HMIN)  H=HM I N 

IF(H.GT.HMHX)  H=HMRX 

L=H/HMIN 

rfX=H/HMIN 

NH= ( HX  +L+ 1 .5  )/2 

H=NH»HMIN 

PRMT( 3 )=H 

N I TER=1 +OELTR/H 

IFINITER.LE.4001)  GO  TO  65 

PRMT(2)=4000.»H+PRMT( 1 ) 

N I TER=4001 

65  WRITE  (6.2005)  (PRMTC  I ) . 1=1  .3  ) 

X0=PRMT( 1 ) 

WR I TE ( 6 .2006  JORLPHfl.OPHI 
C 

C *»*  RERO  TEST  NUMBER  AND  DATE 
C 

REROC 11  .2222  JRPLT 
WR I TE( 6 .2223 1RPLT 
C 

C RERO  HELICOPTER  CONFIGURATION 
C 

RERO ( 5 . 1 1 1 1 ) CONF I G 
WRITE(6.1112)C0NFIG 

*•»  RERO  RECOIL  FORCE  ORTR  FROM  UNIT  11 
C 

CALL  REORECCNH. NITER) 

C 


oou  cj  u u o u u <->  o o 


C — PLOT  RECOIL  FORCE 
C 

HOI VsDELTR/250 . 

DO  68  J=1 .251 
X=XO*(  J-n»HDIV 
CALL  RECOIL! X .RF ) 

YOUT! J)=RF 
68  TIME! J )=X 
WRITE16.2007) 

CALL  RL0T(2S1 .TIME. YOUT) 

CALL  DIME! RR  .SS . TT ) 

WRITE  (6.4000)  SS 

READ  MATRIX  OF  MOOESHAPES 

NS=N0(1) 

NE=N0!6) 

READ (5. 1005)  ( ( A! I . J ) . I=NS .NE ) .J=7.NM0D) 
CALL  DIME! RR .SS.TT ) 

WRITE  (6.4001)  SS 

---  COMPRESS  MOOESHAPE  MATRIX 

NMOO-NMOOE 
DO  160  Irl.NMOD 
00  160  Jrl.NDF 
160  A! J . I )=P( J . I MODE! I)) 

— - ZERO  THE  TEMPORARY  I/O  FILE  14 

NIN  = 13 
NOUT  = 14 
MITER  = NITER/100 

IF(NI TER.  NE-100* MITER)  MITER  = MITER  ♦ 1 
DO  170  1=1 .900 
170  WORK! I ) = 0. 

00  180  1=1. MITER 
WRITE! 14)  WORK 
180  CONTINUE 

SOLVE  MODAL  EQUATIONS 

200  DO  210  1=1 .NMOOE 
KMP=IMOOE( I ) 

0HEOfl=SQR  T! GF ( I MODE  ( I ) ) ) 
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DM 

? 


ALP=-0 -5«S0»0ttE0fl 

BETA=0 .5»0HEGA»SQRT(  4 .-SD»»2  ) 

CALL  QUAOC  NITER  .ALP.BETA.XO.H.I  .KflA  ) 

CALL  OiMECRR.SS.TT) 

WR I TEC  6 . 4006  ) KMA.SS 
210  CONTINUE 

CALL  OUTPCRPLT .CONFIG. DALPHA .DPHI .NITER) 

CALL  OIMECRR.SS.TT) 

999  WR I TEC  6 . 4005  ) TT 

1000  F0RMATC6F10.S) 

1001  FORMATC  61 10  J 

1005  FORMATC  6E12 .5 ) 

1006  FORMAT  CE12.6.5X.E12.6) 

1008  FORMAT  (E12-6.E12.6J 

1111  FORMAT! 3A4 ) 

1700  FORMAT  (E15.8) 

2002  FORMAT  ( 3X . '81 ,B2 .01 .02 */3X  . 4F 12  .6  ) 

2003  FORMAT ( 'O’.SX. 'M0DE'.5X. 'EIGENVALUE') 

2004  FORMAT ( 5X . 1 3 ,6X .E 12 .6 ) 

1112  FORMAT ( IX  . 'CONFIO'JRATION=  ' .3A4  ) 

2005  FORMATC  *0  * . 'RKOS  PARAMETERS : r/5X . ’START  TIME  =’.F8.3/ 

* SX . 'STOP  TIME  ='.F8.3/5X. 'INCREMENT  ='.F6.4) 

2006  FORMATC  'O'. 'ALPHAr ' ,F5 . 1 . IX . 'DEGREES '. 15X . 'PHI= ' ,F5 . 1 . IX . 

$ 'DEGREES ' ) 

2007  FORMATC'  '.’RECOIL  FORCE1) 

2222  F0RMATC5A4) 

2223  FORMATC  '1  ’.5A4) 

4000  FORMAT  (' 1 ’. 'PREPROCESSOR  CPU  TIME  =*.F8.3.'  SECONDS') 

4001  FORMATC  'O'. 'MODESHAPE  READ/WRITE  CPU  TIME  ='.F8.3.'  SECONDS') 

4003  FORMATC  'O'. 'PLOTTER  CPU  TIME  ='.F8.3.’  SECONDS’) 

4004  FORMATC '0*. 'PRINTER  CPU  TIME  ='.f8.3.*  SECONOS  ' ) 

4005  FORMATC  ’O’. ’TOTAL  CPU  TIME  =’.F8.3.’  SECONDS’) 

4006  FORMATC  'O’ . 'CPU  TIME  FOR  INTEGRATION  OF  MODAL  EQUATION '.  13 . 

• ’ ='.F8.3.  ' SECONDS') 

4007  FORMATC 'O' . 'STAT  CPU  TIME  ='.F8.3.'  SECONDS') 

STOP 

END 
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SUBROUTINE  DIME 


PURPOSE: 

CALCULATES  THE  ELAPSED  CPU  TIME  IN  SECONDS  AND  THE  CPU 
TIME  SINCE  THE  LAST  CALL.  SUBROUTINE  NOW  READS  THE 
DIGITAL  CLOCK  AND  RETURNS  THE  CPU  TIME  IN  MILLISECONDS. 

LISTINO  OF  PARAMETERS: 

RR  I/O  PARAMETER  - ON  INPUT  IS  THE  CPU  TIME  AT  THE  LAST 
CALL  OF  DIME:  ON  OUTPUT  IS  THE  PRESENT  CPU  TIME. 


SS  OUTPUT  PARAMETER  - THE  CPU  TIME  INCREMENT  SINCE 
THE  LAST  CALL  OF  DIME. 

TT  OUTPUT  PARAMETER  - TOTAL  CPU  TIME 

SUBROUTINES  CALLED: 

NOW 


SUBROUTINE  OIME(RR.SS.TT) 

CALL  NOW  (I TIME) 

TT=ITIME/lDOO .0 

SS=TT-RR 

RR=TT 

RETURN 

END 
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SUBROUTINE  DQSF 
PURPOSE 

TO  COMPUTE  THE  VECTOR  OF  INTEGRAL  VALUES  FOR  A GIVEN 
EQUIOISTANT  TABLE  OF  FUNCTION  VALUES. 

USAGE 

CALL  DQSF  (H.Y.Z.NOIM) 

DESCRIPTION  OF  PARAMETERS 

H - DOUBLE  PRECISION  INCREMENT  OF  ARGUMENT  VALUES. 

Y - DOUBLE  PRECISION  INPUT  VECTOR  OF  FUNCTION  VALUES. 

Z - RESULTING  DOUBLE  PRECISION  VECTOR  OF  INTEGRAL 

VALUES.  Z MAY  BE  IDENTICAL  WITH  Y. 

NOIM  - THE  DIMENSION  OF  VECTORS  Y AND  Z. 

REMARKS 

NO  ACTION  IN  CASE  NOIM  LESS  THAN  3. 

SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 
NONE 

METHOD 

BEGINNING  WITH  Z(n=0.  EVALUATION  OF  VECTOR  Z IS  DONE  BY 
MEANS  OF  SIMPSONS  RULE  TOGETHER  WITH  NEWTONS  3/8  RULE  OR  A 
COMBINATION  OF  THESE  TWO  RULES.  TRUNCATION  ERROR  IS  OF 
ORDER  H»»S  U.E.  FOURTH  ORDER  METHOD  I . ONLY  IN  CASE  NDIM=3 
TRUNCATION  ERROR  OF  2(2)  IS  OF  ORDER  H«*4. 

FOR  REFERENCE.  SEE 

(1)  F.B. HILDEBRAND.  INTRODUCTION  TO  NUMERICAL  ANALYSIS. 
MCGRAW-HILL.  NEW  YORK/TORONTJ/LONDON . 1956.  PP. 71-76. 

(2)  R.ZURMUEHL.  PRAKTISCHE  MAfHEMATIK  FUER  INGENIEURE  UNO 
PHYSIKER.  SPRINGER.  BERLIN/GOETTINGEN/HEIDELBERG.  1963, 
PP. 214-221 . 


SUBROUTINE  OQSF! H .Y  ,Z  .NDIM ) 


DIMENSION  YCl).Z(l) 

DOUBLE  PRECISION  Y .Z.H.HT .SUM1 .SUM2.AUX.AJX1 . RUX2 


DQSF  10 
.DQSF  20 
DQSF  30 
DQSF  40 
DQSF  50 
DQSF  60 
DQSF  70 
OQSF  80 
DQSF  90 
OQSF  100 
OQSF  110 
DQSF  120 
DQSF  130 
OQSF  140 
DQSF  150 
DQSF  160 
OQSF  170 
DQSF  180 
DQSF  190 
DQSF  200 
OQSF  210 
DQSF  220 
OQSF  230 
DQSF  240 
OQSF  250 
OQSF  260 
OQSF  270 
OQSF  280 
DQSF  290 
DQSF  300 
DQSF  310 
0QSF  320 
DQSF  330 
DQSF  340 
OQSF  350 
OQSF  360 
OQSF  370 
OQSF  380 
.DQSF  390 
DQSF  400 
DQSF  410 
DQSF  420 
DQSF  430 
DQSF  440 
DQSF  450 


i 


o o r>  ci  o o o 


C 


HT= -3333333333333333300»H 
IF( NDIN-6  )7 .8 . 1 


i m 


NDIH  IS  GREATER  THAN  6.  PREPARATIONS  OF  INTEGRATION  LOOP 

1 SUM1=Y( 2 ) ♦Y( 2 ) 

8Uril=8Uf1USUNl 
SUNl=HT*C Y( 1)+SUN1«-Y(3) ) 

AUX1=Y(  4 ) Y ( 43 
AUXUAUX1+AUX1 

AUXl=SUNUHTi«(  Y(  3 ) *AUX1>YC  5 ) ) 

AUX2=hT»(Y(  1 ) +3 .87500*  ( Y(  2 ) -t-YC  5 ) ) <-2 .62500*  ( Y(  3)+Y(4))+Y(6)) 
SUM2=Y( 5 ) ♦Y( 5 ) 

SUH2=SUN2*SUN2 

SUN2=AUX2-HT* ( Y ( 4 ) «-SUM2+Yf  6 ) ) 

Z( 1 JrO.OO 
AUX=Yt  3 ) -*-YC  3 ) 

AUX=AUX*AUX 

Z C 2 )sSUf12-HT»(  Y( 2 )+AUX+Y( 4 ) ) 

Z(3)=SUN1 

Z(4)=SUM2 

IF(N0IN-6)S.5.2 

INTEGRATION  LOOP 

2 DO  4 1=7 .N0IM.2 
SUH1=AUX1 
SUK2=AUX2 

AUX1=Y( I-1)+Y(I-1) 

AUX1=AUX1+AUX1 

AUX1=SUH1+HT»( YC I —2 ) +AUX 1+Y  ( 1 3 ) 

Z( I-2)=SUM1 
IF( I-N0IMJ3.6.6 

3 AUX2=Y(I)*Y( I) 

AUX2=AUX2*AUX2 

RUX2=SUN2+HT»(  Y(  1-1  )*AUX2*Y( 1*1 ) ) 

4 2(1-1 )=SUN2 

5 2(N0IM-1)=AUX1 
Z( NOIM )=AUX2 
RETURN 

6 Z( NOIN-1 )=SUM2 
Z(N0IN)=AUX1 
RETURN 

END  OF  INTEGRATION  LOOP 


DQSF  460 
DQSF  470 
OQSF  480 
DG8F  490 
OQSF  500 
DQSF  510 
DQSF  520 
DQSF  530 
DQSF  540 
DQSF  550 
DQSF  560 
DQSF  570 
OQSF  580 
DQSF  590 
DQSF  600 
DQSF  610 
DQSF  620 
DQSF  630 
DQSF  640 
OQSF  650 
DQSF  660 
OQSF  670 
OQSF  680 
DQSF  690 
DQSF  700 
DQSF  710 
OQSF  720 
DQSF  730 
OQSF  740 
DQSF  7S0 
DQSF  760 
DQSF  770 
DQSF  780 
OQSF  790 
OQSF  800 
DQSF  810 
DQSF  820 
DQSF  830 
DQSF  840 
DQSF  850 
DQSF  860 
DQSF  870 
DQSF  910 
OQSF  880 
DQSF  82° 
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7 IF( NOIM-3 ) 12 . 1 1 .d 

NOIM  IS  EQUAL  TO  4 OR  5 

8 SUf12=l  . 12500»HT* C Yt  1 ) *Y  ( 2 ) +YC  2 )*Y(  2 I ♦Vf  3 ) *-Y{  3 ) + Y(  3 )+Y(  4 )) 
SUMl=Y( 2 )+Y( 2 ) 

sumusumusuhi 

SUf11=HT»(  Y(  1 )+SUMl+Y(  3 ) ) 

zmso.oo 

RUX1=Y( 3 ) *YC3  3 
AUXUAUXWAUX1 

Z(  2 )=SUM2-HT» ( Y ( 2 ) -t-AUX  1 + Y(  4 ) ) 

IF( NQItt-5  ) 10  >9  • 9 

9 AUXt=Y( 4 ) *Y( 4 ) 

AUX1=AUX 1 +RUX1 

2(S)=SUIH*HT»(Y(3)+AUX1+Y(S)  ) 

10  2(3)=SUni 
Z14)=SUM2 
RETURN 

NOIN  IS  EQUAL  TO  3 

11  SUrU=HT»(l.2S00»Ym*Y(2)*Y(2)-.2SD0»Y(3P 
SUM2=Y ( 2 ) *Y( 2 ) 

SUn2=SUM2*SUri2 

Z(3J=H T»(Y( I ) +SUH2+YI 3 ) ) 

Z( 1 )=0.00 

zm=sum 

12  RETURN 
END 


DQSF  900 
OQSF  920 
OQSF  830 
DQSF  940 
OQSF  9S0 
DQSF  960 
OQSF  870 
OQSF  980 
DQSF  990 
DQSF1000 
0QSF1010 
DQSF 1 020 
0QSF1030 
0QSF1040 
0QSF10S0 
0QSF1060 
DQSF1070 
DQSF1080 
OQSF 1090 
DQSF1100 
DQSF 1 1 10 
0QSF1120 
DQSF 1130 
DQSF1140 
OQSF 1150 
0QSF1 160 
0QSF1170 
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SUBROUTINE  LORD 
PURPOSE  : 

CALCULATES  THE  THREE  FORCE  COMPONENTS  AND  THE  THREE 
MOMENT  COMPONENTS  OF  THE  RECOIL  FORCE  VECTOR  THAT  IS 
APPLIED  TO  THE  HEL ICOPTER-TURRET  ATTACH  f’OINT.  RESULTS 
C ARE  STORED  IN  THE  VECTOR  FOR. 

C 

C LISTING  OF  PARAMETERS: 

C X INPUT  PARAMETER  - TIME  IN  SECONDS 

C 

C COMMON  BLOCKS  AND  THE  PROGRAMS  THAT  REQUIRE  THEM: 

C /ANGLE/  LOAD. MAIN 

C /OIST/  LOAD. MAIN 

C /FORCE/  LOAD.  QMODE 

C 

c 

C SUBROUTINES  CALLED: 

C RECOIL 

C 

SUBROUTINE  LOAD(X) 

COMMON  /ANGLE/  ALPHA .PHI ,SA ,CA ,SP .CP 
COMMON  /DIST/  Bl.B2.01.D2 
COMMON  /FORCE/  FORI 6) 

CALL  RECOIL! X .RF  ) 

FORI  1 )=+RF*CA»CP 

FORI  2 j=  RF«CA»SP 

FORI  3 )=-RF»SA 

FX  r RF»C-Bl*SA+B2»CA+02) 

FOR! 4)  = -FX»SP 
FORIB)  = FX*CP 
FORI  6 ) = 0.0 
999  CONTINUE 
RETURN 
END 
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SUBROUTINE  OUTP 


C 
C 
C 

C PURPOSE: 

WRITES  RECOIL  FORCE  ORTR  RND  HELICOPTER  ACCELERArION  DRTR  ON 
— - UNIT  12 

LISTING  OK  PARAMETERS : 

RPLT  ARRAY  CONTAINING  T I TLE  OP  RECOIL  FORCE 
CONFIG  RRRRV  CONTAINING  T I Tl.E  OF  HELICOPTER 
CONFIGURATION 

DRLPHR  ELEVATION  RNQLE  IN  OEGREES 
OPH I RZimiTH  ANGLE  IN  OEGREES 
NITER  SIZE  OF  RRRRVS 

COMMON  8L0CKS  ANO  THE  PROGRAMS  THAT  REQUIRE  THEM: 

/MODE/  MAIN, OlJPr. RHODE. QURO 

/REG/  HAIN  .OUTP  .REDREC 

SUBROUTINES  CALLED: 

NONE 

SUBROUTINE  OUTP  I RnLT .CONF I G .OALPHA .OPH I . MI  TER j 

COMHON/MODE/BI  180 .30  ) ,GH(  30  ) .Of  ( 30  ) .NQ(  6 ) .NHOO.NflF  .MITER  .WORK!  9001 
« .NIN .NOUT 

COMMON  /RED/  X0.H.AI4C03) 

DIMENSION  RPLT(S) .CONFIG! 3) ,THETA!3) . SE C 3 ) ,AG(3) 

WRITE  I 12)  RPLT 
WRITE  112)  CONFIG 
WRITE  (12)  OALPHA. OPHI  .NITER. H.XO 
WR  T TE  (6.100)  RPLT. CONFIG 
WRITE  (6.200  OALPHA. OPHI  .NITER. H.XO 
WRITE  (6.400) 

TRANSFER  ORTA  FROM  THE  TEMPORARY  FILE  NOt'T  TO  THE  PERMANENT  FILE  12 

REWIND  NOL'T 
IPRT  = 0 
DO  20  1 = 1. MITER 
REAO(NOUT)  WORK 

IJMIN  = MINO( 100.NI  IER  - ( 1-1 i»lOO) 

DO  20  J= l. IJMIN 
INOX  = J ♦ C I - l )* ion 
00  ID  K--1.3 

SE  ( K ) = WORKf  K *-9»  ( J-  1 ) ) 


A-ll 
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10  THETA! K ) = P0RK(K*.i.9«f  J-l  < ) 

iprt  = rrR r » 1 

in  IPRT.EQ.il  t IPRTrl 
IFC  IN0X.EQ.Nf  TE«)  I^RTM 

IF(  IPRT.EQ.l  ) URI  TE'9.3001  JNOX .A' INOX ) . THE TA .Ft 
20  CONTINUE 

WRT  TEC  6 .401 ) 

REWINO  NOUT 
IPRTrO 

no  on  r=i .nr ter 

REAO  ' NOI'T  ) WORK 

JJMIN  = HI  NO!  100.N  r TER—  C f -l  J*  1001 

oo  oo  j=i  . r .»ri r n 

INOX  = J * ( I-t )»100 
00  25  Kr.1.3 

SE(K)  = H0RW!K*9-!.'-l  ) ) 

THETAfK)  = PORK  (K  *.■*►?»*(  J-:  )) 

25  ROCK)  = WORK ' <\ -f -9*  ( J-  . j) 

WRrrE'12j  Rf INOK ) .5E . 'META .AO 
IPRT=rPRT*l 
I F ( rPRT.EQ.il  ) IPR-  l 
IFC INOX.EQ. NITER)  rPR'  1 
r ft rpRr.EQ.i ) uRrTLffi.3nm  rNnx.Rt inpxj .ag 
30  CONTINUE 

tOO  FORMRT  ! ‘I ‘.5A4/IX.3A4) 

200  TORflRT  ('  ’ . 'ORLPflR  = * .FR  . 1 . 10X  . 'OPHl  s * .FS  . I . 10X  . "NI  TER  =’. 

» M.IOX.'H  r’.  F8.5.lOX.’XO  z'.IH.K) 

300  FORMAT  ! lX.l4.3X.F9.2.e! JX.EIS.9)) 

400  FORMAT!  *0  * ,T5 . 'I  * . T »2 . ‘RECO JL  ‘ . T2C  . ‘TOOl ‘ . '4  1 . ‘TOOL ’ . T62 . * TOO 
« T/9. ‘XIODl ‘ .T97. ‘ X 1 002  ’ . T 1 1 S . "XI 003  ' ) 

401  FORMAT!  ‘0‘  .TR.  ‘I  ‘ .T12.  *RECOP  * . TxS.  ' THE  TR  , ‘ . T 4-4 . TNETA2  ' .TO  l . 
* ‘THETAO’) 

RETURN 

END 


onnoononooo 
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c 

c 

c 

c 

c 

c 

r. 

c 
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c 

c 
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SUBROUTINE  QMODE 
PURPOSE: 

GENERATES  THE  TOTAL  LORO  VECTOR  FOR  THE  SYSTEM  OF  MOOPL 
VIBRATION  EQUATIONS 

LISTING  OF  PARAMETERS: 

INPUT: 

NITER  NUMBER  OF  DATA  POINTS  ON  THE  TIME  INTERVAL 
H TIME  INCREMENT 

N COLUMN  INflEX  Or  MATRIX  A CONTAINING  THE 
MOOESHAPE  WITH  INDEX  IMOGE(N) 

XO  STARTING  TIME 


v 


COMMON  BLOCKS  AND  THE  PROGRAMS  THAT  REQUIRE  THEM: 
/FORCE/  LOAO.  QMOOE 
/MODE/  MAIN .OUPT. QMOOE. QUAD 

/QR/  QMOOE. QUAD 

/WHY/  MAIN.  QMOOE.  QUAD 


SUBROUTINES  CALLED: 

LOAO 

SUBROUTINE  QMOOEI NI TER .H .N , XO ) 

COMMON/MODE/A ( 180 .30 ) .GMC 30 ) .OF  I 30 ) ,N0( 6 ) .NMOO .NDF 
COMMON  /QR/  QRHSI 4001) 

COMMON  /FORCE/  FORf  6 ) 

COMMON  /WHY/  Y( 90 ) . IMOOEI 30 ) 

DIMENSION  FEC 180) 

00  100  1 = 1. NITER 
QRHSI I )=0.0 
X=XO«-(  1-1  )»H 
CALL  LOAO(X) 

DO  11  J= l .6 
JJ=NO( J) 

11  QRHSI  I )=QRHS(  I ) *-A(  J.J . N )»F0R  I J ) 

100  CONTINUE 
RETURN 
ENO 
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SUBROUTINE  QUHD 


PURPOSE: 

INTEGRATES  UNCOUPLED  MODRL  EQUAT IONS  RNO  DETERMINES  TURRET 
RNOULRR  DISPLRCEMEN  rs 

FILES  13  RNO  14  ARE  USED  TO  TEMPORARILY  STORE  DRTR  BEFORE 
IT  IS  WRITTEN  ON  OUTPUT  FILE  12 


LISTING  OF  PARAMETERS: 

INPUT: 

NITER  NUMBER  Of  DATA  POINTS  ON  THE  TIME  INTERVAL 
ALP  EXPONEN  T I HI  DECAY  FACTOR 
BETA  DAMPED  MODAL  FREQUENCY 
XO  STARTING  TIME 
H TIME  INCREMENT 

N COLUMN  INDEX  Of  MATRIX  A CONTAINING  THE 
MODEShAPE  WITH  INDEX  IMX 
IMX  MODAL  COORDINATE  INDEX 

COMMON  BLOCKS  AND  THE  PROGRAMS  THAT  REQUIRE  THEM: 

/MODE''  MRIN. GUP T. QMODE. QUAD 

/OR/  QMODE .QUAD 

/WHY/  MAIN.  QMODE.  QUAD 


SUBROUTINES  CALLED: 

QMODE 

DQSF 

SUBROUTINE  QIJADINI  TER  .ALP  .BETA  .XO  .H  . N . IMX  ) 

COMMON/MODE/At 180 .30 ) .GM( 30 ) .GF( 30 ) .NGt 6 ) .NMOD .NDF .MI  TER .WORK! 900 ) 
•.NIN.NOUT 

COMMON  /QR/  QRHSI 4001 ) 

COMMON  /WHY/  Y( 90 ) . IMOOEC 30 ) 

DOUBLE  PRECISION  OH . DX .CB1 .CBI .SB  1 .SBI . E T J .QM I .3M0I . QMDDI 
DOUBLE  PRECISION  OQ' 4001) .QUAD1 ( 4001 ) 

generate  the  excitation  force  vector 

CALL  QMODEINI TER.H.N.XO) 

CHECK  IF  MODAL  EQUATION  REPRESENTS  A RIGID  BODY  MOTION 


c 


I F ( IMX.LE.6)  GO  TO  40 


C 

C — FLEXURRL  MODAL  EQUATION 
C 

c 

C FORM  THE  SOLUTION  INTEGRALS 

C 

El=EXPl-ALP»XO) 

OH=H 

DO  10  Ir.l. NITER 
OX=XO<-f  1-1  )*0H 
OE\'=DEXP(  -ALP*OX ) 

QUAOl ( I )=QRHS( I )»OCOS( BE  TAwDX )*DEK 
10  OQC  I )=QRHS(  I)*iOSIN(BETA»nX)*DEX 
CALL  nQSFIDH.QUADl . QUAD1 .NJTER) 

CALL  OQSFCDH.OQ. DO. NITER) 

cbi=be:ai»xo 

SBUDSINCCBl  ) 

CBl=DCOPCCBt  ) 

C 

C READY  THE  INPUT  AND  OUTPUT  FILES  ANO  READ  THE  PARTIAL  SOLUTION 
C FROH  THE  INPUT  FILE 
C 

NIN  = NOUT 
NO’J  T 13 

IFCNOUT.EQ.NIM)  NOUT  = 14 

REUINO  lo 

REWIND  14 

DO  30  J=l. MITER 

READ! NIN ) WORK 

IJMIN  = MINOCIOD. NITER-C  J-l  )*100) 

C 

C FORM  THE  SOLUTION  TO  THE  N TM  MODAL  EQUATION 
C 

DO  20  IJrl.I.ifllN 
I = (J-l)*100  ♦ IJ 
dx=xoi-(  I-;  )*dh 
SBI=0SIN(BE7A»0X) 

CBIcDCOSC  BET A* OX) 

E7Ir-DEXP(ALP»PX) 

QMI  =(  SBI  *aUAD..(  I )-CBl*OQI n»El»( ( Y(NMOP*N)-ALP«Y(N) )»(SBI»CBl- 
$SBl»CBr  )^YCN.*eErR»<(CBI«C81+FB:*S81 ) ) )»ET I/BE7A 
QMDT  =E7Ix(  ( ALP»SRI  ♦BETA»CBI  )■ QUAOl  ( I )- C ALP»CBI-BE  TA*SBI)*DQ(  I ) 

• *El»(  Y(  NMOO.-N  ) »BE  TA»  ( C8I*C8I,+SBI*SB1 ) •*•  I RLP»Y{  NMOO-.N  ) - ( ALP*»2>-BETA 
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$»»J;)»Y(N))»(Sar*CBl“SBl*CBr  ) ) J/BETA 
QHOOI  =QRHS(  r )-2.»RLP»QHDI  -CALP«>*2*-8ETA**2}*QriT 


HDD  THE  SOLUTION  TO  THE  N TH  HOOAL  EQUATION  TO  THE  PARTIAL  SUM 
AND  STORE  THE  RESULTS  ON  THE  OUTPUT  FILE 

DO  20  K=1  .3 

K3  = K*3 

KJ  = K+9«( IJ-l ) 

KJ3  = KJ*3 
Ko6  = KJ*6 

WORKCKJ)  = WORKCKJ)  ♦ A( NQ( K ) .N  )»QHDOI 
W0RKCKJ3)  = W0RKCKJ3)  ♦ A( NOC K3 ) .N )*QHOC I 
20  W0RKCKJ6)  = W0RKCK.I6)  ♦ AC NO( K3  ) ,N )*QHI 
WRITECNOUT)  WORK 
30  CONTINUE 
RETURN 

— - RIOID  BODY  HOOAL  EQUATION 

40  RETURN 
END 
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SUBROUTINE  REDREC 


PURPOSE: 

READS  RECOIL  DATA  AND  GENERATES  WEAPON  RECOIL 
LOADING  AT  THE  RECOIL  ADAPTER  ATTACH  POINTS. 


LISTING  OF  PARAMETERS: 

INPUT  : 

XO  STARTING  T I HE 
H TIME  INCREMENT 

NH  INTEGER  NUMBER  OF  INCREMENTS  OF  .0005  IN  THE 
TIME  STEP 

NITER  THE  OROER  OF  MATRIX  fl 
OUTPUT: 

F RECOIL  FORCE  AT  TIME  X.  A POSITIVE  RECOIL 
FORCE  IS  DIRECTED  TOWARDS  THE  TURRET. 

COMMON  BLOCKS  AND  T HE  PROGRAMS  THAT  REQUIRE  THEM: 

/RED/  MAIN .OUTP  .REOREC 

SUBROUTINES  CALLED: 

NONE 

SUBROUTINE  REORECt NH. NITER) 

COMMON  /RED/  X0.H.A(4003) 

DIMENSION  8(2600) 

EQUIVALENCE  ( All 400 ] .8( l ) ) 

ZERO  THE  RECOIL  FORCE  VECTOR 

DO  9 1=1.4003 

9 AT  I )=0 - 

00  10  1=1.2600 

10  8( I )=0. 

READ  THE  RECOIL  FORCE  DATA  FROM  THE  INPUT  CARO  FILE  II 
REAOC ll .100.EN0=2C)B 

ADJUST  THE  DATA  SO  THAT  IT  CONTAINS  THE  PROPER  TIME  STEP  BETWEEN 
AOJACENT  DATA  POINTS 


c 


20  00  50  N=l. MITER 
100  F0RMAT(8F9.2) 

00  60  N=t. NITER 
L=(N-U«NH/2 
XL=CN-l )»NH/2 . 

IF(L.EQ.XL)  00  TO  45 
Ll=( (M-l )*NH*l)/2 
L2=( ( N-l )»NH»3)/2 
A!  N )=- ( B( L 1 )+RCL2) ) / 2 . 

00  TO  50 
45  R(N!--B(L*n 
50  CONTINUE 
RETURN 

THE  ENTRY  RECOIL  ROUTINE  DETERMINES  THE  RECOIL  FORCE  F AT  AN  ARBITRARY 
POINT  X BY  LINEAR  INTERPOLATION 

ENTRY  RECOIL! X .F  ) 

Y=CX-XO)/H 

L=Y 

N=(L  * Y * 1.61/2  ♦ 1 
FrACN) 

RETURN 

END 
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SUBROUTINE  RLOT 
PURPOSE: 

— fl  GENERAL  PLOT  ROUTINE.  WHICH  PLOTS  A FUNCTION  Y = Y(X)  AT  N DATA 

— POINTS  (X(N).Y(NI) 

LISTINE  OF  PARAMETERS : 

N = NUMBER  OF  DATA  POINTS 
X = VECTOR  OF  ABSCISSA  DATA  POINTS 
Y = VECTOR  OF  ORDINATE  DATA  POINTS 

SUBROUTINES  CALLED: 

NONE 


SUBROUTINE  RLOT  (N.XA) 

DIMENSION  X(N).  Y(N).  AMARKC9) 

INTEOER  ORAPHnOl) 

OATA  I BLANK/ ' V.  IPLUS/*i-V.  iOASH/’-V.  IOOT/’.V.  ISTAR/'«V 

---  FIND  THE  MAXIMUM  AMAX  AND  THE  MINIMUM  AMIN  OF  THE  SET  OF  VALUES 
— • TO  BE  PLOTTED 

AMAX=Y( 1 ) 

AMIN=Y( 1 ) 

CO  20  Jul  .N 

IFCYC J) .LE.AMAX)  00  TO  10 
AMAX=Y( J) 

ID  I F I Y ( J ) .OE.AMIN)  GO  TO  20 
AMIN=Y( J) 

20  CONTINUE 

IF( AMAX  .NE. AMIN)  GO  TO  2S 
WRITE  (6.400)  AMAX 
RETURN 
25  CONTINUE 

WRITE:6.100i  AMIN. AMAX 

C LABEL  INCREMENTS  ON  THE  Y-AXIS 
C 

RANOE=AMAX-AM I N 
STEP=C  AMAX -AMIN )/l0 . 

R0=( ( 0 .O-AMIN ) /RANGE )■ 100 . 

1 0=R0 

I0=( I0*R0*1 .6 )/2 . 
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102=10*2 
DO  30  .1=1.9 

30  RMRRKC J)=RMIN*J»STEP 

WRITE. 6.200 ) ( RMRRKI J ) . J~1 .9 ) 

SCALE  RNO  PLOT  THE  POINTS 

DO  70  J=t.N 
YY=YI J) 

RY=CCYY-RMIN)/RRN0£)*100. 

IY=RY 

IY=(  IY-rRy+1  .5  )/2» 

I2=IY*2 
DO  40  1=1.101 
40  ORRPHI  I 3 = IBLANK. 

JFCIO.OE.P)  00  TO  48 
IF( IY.EQ.1 .OR.IY.EQ.O)  00  TO  44 
DO  42  1 = 1. IY 
42  ORRPHI I )=IDRSH 
44  DO  46  1=1.101.10 
46  ORRPHC I)=IDOT 

ORRPHI I Y*l ) = I PLUS 
00  TO  66 

48  /FI  IP. LE. 100)  00  TO  66 
IF! IY.0T.98)  00  TO  52 
00  50  1=12.101 
50  ORRPHI I )=IDRSH 
52  00  54  1=1.101.10 
54  ORRPHI I )=IOPT 

ORRPHI IY*l)=TPLUS 
00  TO  66 

66  IFl /Y.OT.IO-1 ) 00  TO  69 
00  58  I=IZ.IO 
SB  ORRPHI  n = IORSH 
00  TO  62 

59  IFIIY.LE.IO^l)  00  TO  62 
DO  60  1 = 102.  IY 

60  ORRPHC I )=IDRSH 
62  DO  64  1=1.101.10 
64  ORRPHI  n = ioor 

ORRPHI 10*1 )=ISTRR 
ORRPHI I Y*1 JrlPLUS 
66  WRI TEC  6 .300  ) XI J ) .Y( J ) .ORRPH 
70  CONTINUE 

100  FORMAT  I •0*.T2l.  •MINIMUM  = ' .E15 .8 . T97 . 'MAXIMUM  = \E15.8) 
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APPENDIX 

FORTRAN  IV  Source  Listing  of  the  Gun-Turret  Dynamic  Model 


TURRET  DYNAMIC  MODEL!  MAIN  PROGRAM 


PURPOSE: 

DETERMINES  THE  DYNAMIC  MOTIONS  OF  THE  MAJOR  COMPONENTS  OF  THE 
TURRET  AND  CALCULATES  THE  ERROR  IN  THE  GUN  POSITION  OURINO  THE 
WEAPON  FIRING  INTERVAL. 


COMMON  BLOCKS  AND  THE  PROGRAMS  THAT  REQUIRE  THEM: 
/ACCEL/  MAIN.  IFORCE 
/ANGLE/  MAIN.  IFORCE.  OUTP.  RFORCE 
/ANGROT/  MAIN.  OUTP 
/DIST/  MAIN.  IFORCE 
/MASS/  MAIN.  FCT.  IFORCE 

/MOTION/  MAIN.  CALCOM . IFORCE.  OUTP.  RFORCE 
/NPLOT/  MAIN.  CALCOM.  OUTP 

/NUMB/  MAIN.  CALCOM.  IFORCE.  OUTP.  POINT.  RFORCE 
/RED/  MAIN.  CALCOM.  RFORCE 
/STIFF/  MAIN.  FCT 


SUBROUTINES  CALLED: 
CALCOM 
PLOTS 
RKGS 
STATIS 


COMMON  /ACCEL/  T 1 ( 2600 ) . T2(  2600  ).  T3C  2600  ) .UDD(  2600 ) . 

$ VDO( 2600 ) . WOO! 2600 ) 

COMMON  /ANGLE/  PHI .ALPHA . SA .CA .SP .CP 
COMMON  /ANGROT/  AG1 ( 2600 1 .A02( 2600 ) .AG3( 2600 ) 

COMMON'/O IST/A 1 . A2 . A3 .81 .82 
COMMON/MASS/AM. All .AI2.GR 
COMMON/MOT  I ON/ 1 M ,N 

COMMON  /NPLOT/  NP. TIME (2600) .ANG( 2600 ) .WEAPON! 8 ) .OPHI .DALPHA . 
« RATE. ROUNDS 

COMMON  /NUMB/  NITER .HDAT ,X0 ,RPLT( 5 ) .CONFIG! 3 ) 

COMMON  /RED/  FREC12600) 

COMMON/STIFF/TK.TG 

DIMENSION  PRMT! 5 ) .DERY! 5).AUX(8»5) ,Y(5) 

EXTERNAL  FCT. OUTP 

CALL  PLOTS  ( IBUF.IDUM.14) 
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C 

C 


READ  PLOT  TITLE  PARAMETERS 


READ( 5 .1007 ) WEAPON 
READ! 5 .1005 ) RATE .R0UN0S 

SET  MOTION  TO  BE  FIRST. IM=1  AZIMUTH. IM=2  ELEVATION 
READ! 5 .1002 ) IM 

SET  NUMBER  OF  MOTIONS  TO  BE  ANALYZED. 1 OR  2 
READ! 5. 1002)  NM 

READ  CENTER  OF  ORAVITY  OFFSET  FEET 

READ! 5 . 1005 ) A1  .A2  .A3 ,B1 .82 

READ  DATA  FROM  UNIT  1 

READ  ( 1 ) RPLT 
READ  (1)  CONFIG 

REAO  (1)  DALPHA.OPHI .NI7ZR.H0AT.X0 

TURRET  INTERFACE  GRI0  POINT  PARAMETERS: 


UDD  = 

-XDD 

LINEAR 

ACCELERATIONS 

VOD  = 

ZDD 

LINEAR 

ACCELERATIONS 

WOO  = 

YDO 

LINEAR 

ACCELERATIONS 

T1  = - 

TODX 

ANGULAR 

ACCELERATIONS 

T2  = 

TDDZ 

ANGULAR 

ACCELERATIONS 

T3  = 

TDDY 

ANGULAR 

ACCELERATIONS 

AG1  = 

-TX 

ANGULAR 

POSITIONS 

AG2  = 

TZ 

ANGULAR 

POSITIONS 

AG3  = 

TY 

ANGULAR 

POSITIONS 

DO  4 1=1 .NITER 

REAOI 1 ) FRECI I ) .UDDC I ) .VDO( I ) .WDD( IJ.TICJ)  .T2( I ) ,T3( I ) . 
• AGim.AG2m.A03m 
4 CONTINUE 

PRINT  INPUT  DATA 

WR I TE ( 6 . 1 00 ) WEAPON .RPLT. RATE . ROUNDS . CONFIG 


J 
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WiRITEC  6 .200  ) DALPHA.DPHI  .NITER .HOAT ,X0 
WRI TEC  6 . 400 ) 

IPRT  = 0 

DO  1 1 = 1. NITER 

IPRT  = IPRT  ♦ 1 

IFC IPRT.EQ.il)  IPRT  = 1 

IFCI.EQ. NITER)  IPRT  = 1 

IFt IPRT -EQ - 1 ) WRITE! 6.300 ) I ,FREC( I).Tt(I).T2(I).T3(I). 

• UDO( I ) . VDD( I ) . WDDt I ) 

1 CONTINUE 
IPRT  = 0 

WR I TEC  6 . 401  ) 

DO  2 1=1. NITER 
IPRT  = IPRT  ♦ 1 
IFC IPRT.EQ.il  ) IPRT  = 1 
IFC I .EQ. NITER)  IPRT  = 1 

IFC IPRT. EQ.l ) WRITE! 6 .300 ) I .FRECC I ) .AG1 ( I) .AG2CI) .AG3C I ) 

2 CONTINUE 
C 

C INITIALIZE  PARAMETERS 
C 

PI=3. 1415926/180. 

PHI=DPHI*PI 
ALPHA=DALPHA»P  I 
SArSINCALPHA) 

CA=COS( ALPHA) 

SP=SIN( PHI ) 

CP=C0S( PHI ) 

PH  I =PH I*PI 
ALPHA=ALPHA*PI 
C 

C CONVERT  LINEAR  ACCELERATIONS  TO  FT/SEC»«2  UNITS 
r 

DO  5 1 = 1 .NITER 
UDDC I )=UOD( I ) / 1 2 - 
VDDC I )=VDD( I ) / 1 2 . 

5 WOOC I )=WDO( I )/ 12 . 

00  99  N=1.NM 
NP=0 

IFCIM.EQ.N)  GO  TO  12 
WR I TEC  6 . 1 600 ) 

GO  TO  13 

12  WR I TEC  6 . 1500 ) 

13  CONTINUE 
C 
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C RERO  MASS  IN  LBS  RNO  MOMENTS  OF  INERTIR  OF  LORO  RNO  MOTOR 
C MOMENTS  OF  INERTIR  IN  SLU0S-FT««2 
C 

RERD( 5 . 1000 1RM .RI 1 .RI2 
RM=RM/32.2 

RERD  OERR  BOX  RATIO 

RERDC 5 >1000 )GR 

RERO  TRANSMISSION  ST IFFNESS .FT-LBS . RNO  DAMPING  COEFFICIENT 

RERD( 5. 10031TK .TG 
WR I TEC  6 <2000 1 RM.RI1.RI2 
WRITE! 6.2001 ) Rl  .R2.R3.B1 .82 
WRI TEC  6 .2002 1 TK  .TG 
WR I TEC  6 .2003  ) GR 

CONVERT  MOTOR  MOMENT  OF  INERTIR  TO  TURRET  COORDINATES 
RI2=fiI2*CGR)»*2 
SET  INITIAL  DISPLACEMENTS 

00  10  1=1.5 

10  Y ( I 1=0.0 

SET  ERROR  WEIGHTS 

00  11  1=1.5 

11  DERYC I )=.2 

READ  THE  UPPER  ERROR  80UN0 

RERDC  5.1001 ) PRMTC4) 

PRMTC 1 )=X0 

PRMTC  2 )=X0+( NI TER-1  )»HDRT 
PRMTC  3)=HDRT 
WR I TEC  6 .2005 ) 

CALL  THE  RUNGE-KUTTR  EQURTION  SOLVER  RKGS. 

CRLL  THE  STATISTICAL  ERROR  PARAMETER  SUBROUTINE  STRTIS. 

C CRLL  THE  CRLCOMP  PLOT  SUBROUTINE  CRLCOM 
C 

20  CRLL  RKGSCPRMT.Y.0ERY.5.IHLF.FCT.0UTP.RUX) 


B-4 


I 


CALL  STATISC NITER. AND. ANOAV.ANOSO) 

CALL  CALCOMC  NM.ANOAV  .ANOSD .PRMT  ) 

IF( IM.EQ.N ) 00  TO  31 
WRITEC6.2008) 

WRI TEC 6.2007 ) 

00  TO  32 

31  CONTINUE 
WRI TEC  6 . 20C6  ) 

WRI TEC  6 .2007  ) 

32  CONTINUE 

WRI TEC 6 . 1 006 ) ANOAV. ANOSD 
99  CONTINUE 
100  FORMATC  '1  * .8A4/1X.5A4/ 

• ’ FIRINO  RATE  =’.F4.0.’  SPM;  ROUNDS  = ’ ,F4 .0/1X .3A4 ) 

200  FORMATC’  ’.’DALPHA  = ’ .F5 . 1 . 10X . *DPHI  =’ ,F5 . 1 . 10X . ’NITER  = ’. 

• I4.1DX.’H  = ’.F8.5.10X. 'XO  =’.F8.S) 

300  FORMATC IX . 1 4 .3X .F9 .2 ,6( 3X .E15  .8  )) 

400  FORMATC  ’C’.T5, ’I ’ .T12 . ’RECOIL ’ .T26 . ’TDD1 ’ .T44 . ’T0D2  ’ . T62 . ’T003  ’ . 

• T79.’XIDD1’.T97.’XIDD2’.T11S. ’X IDD3 ’ ) 

401  FORMATC 'O’ ,T5. ’I ’ . T12 . ’RECOIL ’ . T25  » ’ THETA 1 ’ . T 43 , ’THETA2  ’ .T61 . 

• ’THETA3  ’ 1 

1000  F0RMATCF10.5.2E1S.5) 

1001  FORMATC  2F10 .3 ) 

1002  FORMATC 110) 

1003  F0RMATCE10.3.F10.S) 

1005  F0RMATC5F10.3) 

1006  F0RMATC5X.E1S.6.10X.E15.6) 

1007  FORMAT (2004) 

1600  FORMATC '1 ’.SOX. 'AZIMUTH  VIBRATION’) 

1600  rORMRTC  ’1 ’.SOX. 'ELEVATION  VIBRATION’) 

2000  FORMATC 7X. ’MASS  TURRET  INERTIA  MOTOR  INERTIA ’/4X ,3E18 .6 ) 

2001  F0RMATC7X. 'CENTER  OF  GRAVITY  AND  PIVOT  POSI TIONS ’/4X .5F18 .6 ) 

2002  FORMATC 7X. ’GEAR  STIFFNESS  DAMPING  COEFFICIENT  ’/4X .2E18 .6 ) 

2003  FORMATC 4X. ’ GEARRATIO  ’/4X.F10.0) 

2005  FORMAT! 9X . ’TIME  . ANGULAR  DISPLACEMENT  IN  RADIANS’) 

2006  F0RMATC5X. 'AZIMUTH  MOTION') 

2007  FORMATC 5X. ’MEAN  VALUE  STANDARD  DEVIATION’) 

2008  FORMATCSX. 'ELEVATION  MOTION’) 

STOP 

END 
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C SUBROUTINE  CALCOM 

C 

C 

C PURPOSE: 

C PLOTS  THE  AZIMUTH  AND  ELEVATION  ANGLES  AND  THE  RECOIL  FORCE  VS.  TIME 
C ON  THE  CALCOMP  914  PLOTTER 

C 
C 

C LISTING  OF  PARAMETERS: 

C NM  = NUMBER  OF  MOTIONS  TO  BE  PLOTTED 

C ANOAV  = MEAN  VALUE  OF  THE  DYNAMIC  MOTION 

C ANGSD  = STANDARD  DEVIATION  OF  THE  DYNAMIC  MOTION 

C PRMT  = PARAMETER  ARRAY  CONTAINING  THE  START  TIME.  STOP  TIME. 

C TIME  INTERVAL.  AND  THE  UPPER  ERROR  BOUND 


COMMON  BLOCKS  AND  THE  PROGRAMS  THAT  REQUIRE  THEM: 

/MOTION/  MAIN.  CALCOM . IFORCE.  OUTP  . R» OR* E 
C /NPLOT/  MAIN.  CALCOM.  OUTP 

C /NUMB/  MAIN.  CALCOM.  IFORCE.  OUTP.  POINT.  RFORCE 

C /RED/  MAIN.  CALCOM.  RFORCE 

C 
C 

C SUBROUTINES  CALLED: 

C AXIS 

C FACTOR 

C GBDATE 

C LINE 

C NEWPEN 

C NUMBER 

C PLOT 

C RECT 

C SCHLE 

C SYMBOL 

C 
C 

SUBROUTINE  CALCOM!  NM. ANGAV .ANGSD .PRMT ) 

COMMON  /MOTION/  IM.N 

COMMON  /NPLOT/  NP . T IME(  2600 ) ,RNG( 2600 ) .WEAPON! 3 ) .DPHI .DALPHA . 
* RATE. ROUNDS 

COMMON  /NUMB/  NI TER ,HOAT ,X0 . RPLT! 5 ) .CONFIG! 3 ) 

COMMON  /REO/  FREC12600) 

DIMENSION  DAT(7 ) .PRMT! S ) 

IFINP.GT.2598)  NP=2598 
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INITIALIZE  THE  NEXT  PLOT 


CALL  PL0T10..-15..-3) 

CALL  PLOT! 0 . . 1 4 . .2 ) 

CALL  PLOT! 1 . .3 .37S ,-3 ) 

CALL  FACTOR (0. 75) 

CALL  RECT(0..0..  9..12..0..2) 
CALL  PL0TU..1..-3) 


SCALE  AND  DRAM  THE  AXES 


N1=NP+1 

N2=NPf2 

CALL  SCALE! ANG.7 . .NP , 1 ) 

TIHE(Nl)  = PRMT(l) 

TIME!  N2 ) = (PRMT(2)  - PRMT(1))/10. 

CALL  AXISC0..0.. ’TIME  ( SECONDS  )’  .-14 . 10  . .0 . .TIME! Nl) .TIME! N2 ) ) 
CALL  GBORTE(OAT) 

CALL  SYMBOL! 6. ,-0.8. 0.14. OAT .0. ,28) 

CALL  AXIS(0. .0. , ’ANGLE  ( MRADS  ) * . 13 .7 . .90 . .AND! Nl) .ANO! N2 )) 


DRAW  THE  ZERO  AXIS 


ZEROP=ABS( ANG! N l )/ANO( N2 ) ) 
CALL  PLOTIO. .ZER0P.3) 

CALL  PLOT! 10 . . ZEROP .2 ) 


WRITE  THE  HEADING  FOR  THE  AZIMUTH  OR  ELEVATION  ANOLE  PLOT 


IF!  IM 
CALL 
GO  TO 
20  CALL 
30  CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 


.EQ.N) 

SYMBOL! 

30 

SYMBOL! 

SYMBOL! 

SYMBOL! 

SYMBOL! 

SYMBOL! 

NUMBER! 

SYMBOL! 

NUMBER! 

SYMBOL! 

SYMBOL! 

NUMBER! 

SYMBOL! 


GO  TO  20 

2.5.7.38.0.21, ’ELEVATION  ’.0..10) 


,7.38.0.21. ‘AZIMUTH  ’.0..8) 
..999... 21. 'ANGLE  VS.  TIME*. 0.14) 
,7.  .0.14. ’RECOIL  ORTA:  ’.0..13) 


999. 

.0. 

14 

.RPLT.O. , 

20) 

.7.0 

.14 

« 1 

• 1 

ORIENTATION: 

999. 

.0. 

14 

.OPHI .0. . 

1 ) 

999. 

.0. 

14 

. ’ DEO  AZ 

; ’. 

999. 

.0. 

14 

.OALPHA.O 

..1) 

999. 

.0. 

14 

. ’ DEO  EL 

’.0. 

.4.0 

.14 

• 

MEAN  =’.0 

. .6) 

999. 

.0. 

14 

.ANGAV.O. 

.2) 

999. 

.0. 

14 

. ’S  STO. 

DEV 

CALL  NUMBER ( 999 . >999 . ,0 . 14 , ANDSO ,0 . ,2 ) 


PLOT  THE  AZIMUTH  OR  ELEVATION  ANOLE  VS.  TIME 


CRLL  NEHPEN! 2 ) 

CALL  LI NEC  TIME. AND. NP .1.0.1 ) 
CALL  NEHPEN! 1) 

CALL  FACTORC 1 .0) 

CALL  PL0TC9 .25 .-15 . .-3 1 
IFCN.NE.NM)  RETURN 


PLOT  RECOIL  FORCE 


CALL  PLOTC 0 . . 14 . .2 ) 

CALL  PLOTC 1.. 3. 375. -3) 

CALL  FACTORC 0.75) 

CALL  RECTCO..O..  9..12..0..2) 

CALL  PL0TC1..1..-3) 

CALL  SCALECFREC.7..NP.1 ) 

CALL  AXISC  0 . .0 . . ’TIME  ( SECONOS ) * .-14 . 10 . .0 . .TIMEC N1 ) .TIMEC N2 ) ) 
CALL  SYMBOL! 6.. -0.8. 0. 14. DAT. 0.. 28) 

CALL  AXISCO. .0. . ’RECOIL  FORCE  ( LBS ) ’ . 18 .7 . .90 . . 

« FRECC N1 ) .FRECC N2 ) ) 

ZEROP  = ABSC  FRECC  N1 J/FRECC  N2 ) ) 

CALL  PLOTC 0 . .ZEROP .3 ) 

CALL  PLOTC 10. .ZEROP. 2) 

CALL  SYMBOL ( 2 .0.7.38,0.21 . 

• ’TRANSMITTED  RECOIL  FORCE  VS.  TIME  ’ .0 . .33 ) 

CALL  SYMBOL! 2 .0.7.0,0.14 .MEAPON .0 . .32 ) 

CALL  SYMBOL! 2.0.6  .7 .0 . 14 .RPLT .0 . .20 ) 

CALL  SYMBOL! 2 .0.6.4.0.14, ’FIRINO  RATE  = ’.0..14) 

CALL  NUMBER! 999.. 999.. 0.1 4. RATE. 0., 1 ) 

CALL  SYMBGU999.  .999.  .0.14.’  SPM;  ROUNDS  = ’ .0  . . 16  ) 

CRLL  NUMBER! 999 . . 999 . . 0 . 1 4 . ROUNDS . 0 . . 1 ) 

CALL  NEWPENC2) 

CALL  LINECTIME.FREC.NP, 1,0.1) 

CALL  NEHPEN! 1) 

CALL  FACTOR! 1 .0) 

CALL  PLOT! 9. 2S.- 15. .-3) 

CALL  PLOTC 0 . . 12 . .2 ) 

CALL  PL0TC4..0. .999) 

RETURN 

END 
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SUBROUTINE  FCT 


I 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


PURPOSE: 

EVALUATES  THE  RIOHT  HAND  SIOE  OF  THE  QOVERNINO  EQUATIONS 


LISTINO  OF  PARAHETERS; 

X = TIME  INPUT 

Y = VECTOR  OF  SOLUTIONS  AT  TIME  X 

DERY  = VECTOR  OF  TIME  DERIVATIVES  OF  THE  SOLUTION  VECTOR  Y AS 

DEFINED  BY  THE  RIOHT  HAND  SIDE  OF  THE  OOVERNINO  EQUATIONS 


COMMON  BLOCKS  AND  THE  PROGRAMS  THAT  REQUIRE  THEM: 
/MASS/  MAIN.  FCT.  IFORCE 
/STIFF/  MAIN.  FCT 


SUBROUTINES  CALLED: 
IFORCE 
RFORCE 


SUBROUTINE  FCTIX.Y .DERY ) 

COHMON/MASS/AM . A 1 1 . A 1 2 . OR 
COMHON/STIFF/TK . TO 

DIMENSION  PRMT( 5 ) .DERY! 5 ) . AUX( 8 .5 ) .Y( 5 ) 

CALL  SUBROUTINE  RFORCE  AND  IFORCE  TO  OET  THE  RECOIL  AND  INERTIAL  FORCES 


CALL  RFORCE (X.TR) 

CALL  IFORCE! X . TH ) 

P=0.0 

AF=ABS( Y( 2 ) ) 

IF(AF.EQ.O.O)  AF=I.O 

T=4 .0125E*02«»(  P-Y!  1 ) — Y (51) 

TP=ABSm 

IF(TA.OT.2.6)  T=T/TA«2.6 
C=T«OR 


EVALUATE  RIOHT  HAND  SIDE  OF  OOVERNINO  EQUATIONS 


DERY! 1 )=Y! 2 ) 
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OERY( 2 )=-TK/RI 1»( Y( 1 ) — Y ( 3 ) ♦ Y C 2 )/RF»TO)*( TR*TH )/RI 1 
DERYC  3 )=YC  4 ) 

0ERY(4)=TK/RI2»(Y( 1 )-Y(3) ) +C/RI2 

0ERY(6J=(4.785047E-S»Y( 4)*1 .01U86E-4»0ERY(  4)  )»0R-1D  .■Y(S) 
RETURN 
END 


o o o o 


m 


f 


c 

C SUBROUTINE  IFORCE 
C 

c 

C PURPOSE: 

C CRLCULRTES  THE  INERTIRL  TORQUES  CAUSED  BY  THE  ACCELERATIONS  OF  THE 
C TURRET  RTTRCHMENT  PLATFORM 

C 
C 

C LISTING  OF  PARAMETERS: 

C X = TIME  INPUT 

C TH  = INERTIAL  TORQUE  ABOUT  THE  ORIGIN  OF  THE  TURRET 

C COORDINATE  SYSTEM 

C 

r 

C COMMON  BLOCKS  AND  THE  PROGRAMS  THAT  REQUIRE  THEM: 

C /ACCEL/  MAIN.  IFORCE 

C /ANGLE/  MAIN,  IFORCE.  OUTP.  RFORCE 

C /DIST/  MAIN.  IFORCE 

C /MASS/  MAIN.  FCT.  IFORCE 

C /MOTION/  MAIN.  CALCOM.  IFORCE.  OUTP.  RFORCE 

C /NUMB/  MAIN.  CALCOM.  IFORCE.  OUTP.  POINT.  RFORCE 

C /TORQUE/  IFORCE.  OUTP.  RFORCE 

C 
C 

C SUBROUTINES  CALLED: 

C POINT 

C 

c 

SUBROUTINE  IFORCE(X.TH) 

COMMON  /ACCEL/  T 1 ( 2600 ). T2( 2600  ).  T3C 2600  ) .UOOI 2600 ) . 

$ VQDI 2800 ) .WDD( 2600 ) 

COMMON  /ANGLE/  PHI .ALPHA .SA .CA  .SP .CP 
C0MM0N/DIST/A1 .A2.A3.B1 .82 
COMMON/MASS/AM. All .AI2.GR 
COMMON/MOT  I ON/ 1 M.N 

COMMON  /NUMB/  NITER. HDAT .XO.RPLT(S) .CONFIG! 3 ) 

COMMON/TORQUE/TR .THOUM 

CALCULATE  THE  ANGULAR  AND  LINEAR  ACCELERATIONS  OF  THE  TURRET 
ATTACHMENT  AT  TIME  X 

CALL  POINTCX.NUM.XNJ 

ThX=-Tl(NUM)  + (Tl(NUM)-Tl(NUM*l ) )■( X-XN J/HDAT 
THYrn(NUMW(T3(N'JM*l  )-T3(  NUM  ))■(  X-XN  J/HDAT 
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THZ=T2( NUM ) ♦( T2( NUM*l )-T2C NUM ) )■( X-XN )/HDRT 
RX=-UDD(  NUM  ) ♦( UDDC  NUfl ) -UDD(  NUM«-1 ) )»( X-XN )/H0RT 
flYrwoo(Nurn«-(woo(Nurt*i  )-udo(num)  )»(x-xn)/hdrt 
RZsVOOC  NUM ) ♦( V00( NUM+1  )-VOO(NUMn»(X-XN)/HORT 
IF( IM.EQ.N)  00  TO  9 

FORCE  FOR  ELEVATION  DIRECTION 

TH=-RM»(  ( R2+R3»CR  )■  RY-(RX*CP-flZ»SP)«(-B2fR3»Sfin-RIl» 
$ ( THX*SF ♦THZ»CP  ) 

00  TO  10 

9 CONTINUE 
C 

C FORCE  FOR  AZIMUTH  DIRECTION 
C 

TH=RM»fll ■( RX»SP*RZ*CP )-RI 1 »THY 

10  CONTINUE 
THDUM=Th 
RETURN 
END 


c 

C SUBROUTINE  OUTP 

C 

C 

C PURPOSE: 

C FILLS  THE  OUTPUT  ARRAY  ANG.  WHICH  IS  USED  IN  SUBROUTINE  CALCOM  TO 

C PLOT  THE  AZIMUTH  AND  ELEVATION  ANGLES.  ALSO  THE  TIME.  RELATIVE 

C TURRET  ANGLE.  AND  THE  DYNAMIC  MOTION  OF  THE  GUN  ARE  PRINTEO. 

C 

C 

C LISTING  OF  PARAMETERS : 

C X = TIME  INPUT 

C Y = SOLUTION  VECTOR 

C DERY  = TIME  DERIVATIVE  OF  THE  SOLUTION  VECTOR 

C IHLF  = OUTPUT  PARAMETER  FROM  SUBROUTINE  RKGS.  WHICH  SPECIFIES 

C THE  NUMBER  OF  BISECTIONS  OF  THE  INITIAL  INCREMENT 

C NDIM  = THE  NUMBER  OF  EQUATIONS  SOLVED  BY  SUBROUTINE  RKOS 

C PRMT  = PARAMETER  ARRAY  CONTAINING  THE  START  TIME.  STOP  TIME.  INITIAL 

C INCREMENT.  AND  THE  UPPER  ERROR  BOUNO 

C 
C 

C COMMON  BLOCKS  AND  THE  PROGRAMS  THAT  REQUIRE  THEM: 

C /ANGLE/  MAIN.  IFORCE.  OUTP.  RFORCE 

C /ANOROT/  MAIN.  OUTP 

C /MOTION/  MAIN.  CALCOM.  IFORCE,  OUTP.  RFORCE 

C /NPLOT/  MAIN,  CALCOM.  OUTP 

C /NUMB/  MAIN.  CALCOM.  IFORCE.  OUTP.  POINT.  RFORCE 

C /TORQUE/  IFORCE.  OUTP.  RFORCE 

C 
C 

C SUBROUTINES  CALLED: 

C POINT 

C 

C 

SUBROUTINE  OUTP ( X ,Y , DERY . IHLF .NDI M .PRMT ) 

COMMON/ANGLE/  PHI  , ALPHA  . SA  .CA  .SP .C? 

COMMON  /ANGROT/  AOl ( 2R00 ) . A02 I 2S0D ) , AG3 ( 2600 ) 

COMMON/MOTION/  IM.N 

COMMON  /NPLOT/  NP ,TI ME( 2600 ) .ANOC 2600 ) 

COMMON  /NUMB/  NITER. HDAT.XO.RPLT(S) .C0NFIGI3) 

COMMON/TORQUE/TR.TH 
DIMENSION  Y ( 5 ) .DERYIS)  .PRMT(S) 

NP=NP*1 

IFCNP.OT .2600)  GO  TO  9 
TIME! NP )=X 
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CALCULATE  THE  ANGULAR  POSITIONS  AT  THE  TURRET  ATTACHMENT  POINT 
AT  THE  TIME  X 

CALL  POINT(X.NUM.XN) 

G1  = AOl(NUM)  CAOHNUM  + 1)  - AOl(NUM))  a (X-XNJ/HOAT 

G2  = AG2CNUM)  + (AG2(NUM*n  - AG2CNUM))  a (X-XNJ/HDAT 

03  r A03CNUM)  + CA03(NUM+n  - A03(NUMJ)  a (X-XNl/HDAT 

IFCIM.EQ.N)  GO  TO  8 

ELEVATION  ANGLE 

ANO(NP)  = ( Y( 1 ) - GlaSP  ♦ G2aCP)al000. 

GO  TO  9 
8 CONTINUE 

AZIMUTH  ANGLE 

ANGINP)  = !(Y(1)  ♦ 03  JaCA  ♦ (GlaCP  + 02«SP  )»SA  )*1000 . 

9 CONTINUE 

IF  (NP.EQ.I)  GO  TO  10 

IF( (NP/10J-10.NE.NP)  RETURN 

10  CONTINUE 

WRITE!  6 . 1000  ) X.YU  ) .ANG1NP  ) 

1000  FORMAT! 3X .6E15 .6  ) 

RETURN 

END 


— 


I 


c 

C SUBROUTINE  POINT 
C 

C 

C PURPOSE  : 

C FINOS  THE  TIME  XN  AND  THE  CORRESPONOINO  INDEX  NUH  OF  AN  ARRAY , 

C SUCH  THAT  XN  IS  JUST  BELOW  X-  THESE  VHLUES  ARE  USEO  FOR  LINEAR 

C INTERPOLATION. 

C 

C 

C LISTING  OF  PARAMETERS: 

C X = TIME  INPUT 

C NUM  = CALCULATED  INDEX  OF  THE  ARRAY  SUCH  THAT  THE  TIME  XN  IS 

C JUST  BELOW  X 

C XN  = TIME  CORRESPONOINO  TO  THE  INDEX  NUM 

C 

C 

C COMMON  BLOCKS  AND  THE  PROGRAMS  THAT  REQUIRE  THEM: 

C /NUMB/  MAIN.  CRLCOM.  IFORCE.  OUTP.  POINT.  RFORCE 

C SUBROUTINES  CALLED: 

C NONE 

C 

C 

SUBROUTINE  POINT! X .NUM .XN ) 

COMMON  /NUMB/  NI TER .HQAT .XO ,RPLT( S ) .CONFIG! 3 ) 

AN=( X-XO ) /HDAT+l  . 

NUM=AN 

NUM=( AN+NUM+1 .S)/2 

IFtNUM.EQ. NITER)  NUM=NIT£R-1 

XN=XO+! NUM-1 )*HDAT 

RETURN 

END 


3 
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SUBROUTINE  RFORCE 


PURPOSE: 

CALCULATES  THE  APPLIED  TORQUES  GENERATED  BY  THE  WEAPON  RECOIL  FORCE 


LISTING  OF  PARAMETERS: 

X = TIME  INPUT 

TR  = CALCULATED  RECOIL  TORQUE 


COMMON  BLOCKS  AND  THE  PROGRAMS  THAT  REQUIRE  THEM: 
/ANGLE/  MAIN.  IFORCE.  OUTP.  RFORCE 
/MOTION/  MAIN.  CALCOM.  IFORCE.  OUTP.  R*  OR*  E 
/NUMB/  MAIN.  CALCOM.  IFORCE.  OUTP.  POINT.  RFORCE 
/RED/  MAIN.  CALCOM.  RFORCE 
/TORQUE/  IFORCE.  OUTP.  RFORCE 


SUBROUTINES  CALLED: 
POINT 


SUBROUTINE  RFORCE! X . TR  ) 

COMMON  /ANGLE/  PHI  .ALPHA .SR .CA . SP  .CP 
COMMON/MOT  I ON/ 1 M.N 

COMMON  /NUMB/  N I TER .HOAT ,XO . RPLT! 5 ) .CONFIG! 3 ) 
COMMON  /RED/  FRECC  2600 ) 

COMMON/TORQUE/TROUM ,TH 

DETERMINE  THE  WEAPON  RECOIL  FORCE  RT  THE  TIME  X 
CALL  POINT ! X ,NUM .XN ) 

F=FREC(NUM)  + !FREC!NUM+1  )-FREC(  NUM  ) )ti(  X-XN  J/HORT 
IF! IM.EQ.N)  GO  TO  9 


C 


9 

C 

C 


FORCE  FOR  ELEVATION  DIRECTION 


TR=F»0.16 
GO  TO  10 
CONTINUE 
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FORCE  FOR  AZIMUTH  DIRECTION 


TR=F»0 . 16»Cfl 

CONTINUE 

TRDUI1=TR 

RETURN 

END 
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SUBROUTINE  RKGS 

PURPOSE 

TO  SOLVE  fi  SYSTEM  OF  FIRST  ORDER  ORDINARY  DIFFERENTIAL 

EQUATIONS  WITH  GIVEN  INITIAL  VALUES. 

USAOE 

CALL  RKOS  (PRMT. Y.OERY.NOIM.IHLF.FCT. OUTP. AUX) 

PARAMETERS  FCT  AND  O'JTP  REQUIRE  AN  EXTERNAL  STATEMENT. 

DESCRIPTION  OF  PARAMETERS 

PRMT  - AN  INPUT  AND  OUTPUT  VECTOR  WITH  DIMENSION  GREATER 
OR  EQUAL  TO  5,  WHICH  SPECIFIES  THE  PARAMETERS  OF 
THE  INTERVAL  AND  OF  ACCURACY  AND  WHICH  SERVES  FOR 
COMMUNICATION  BETWEEN  OUTPUT  SUBROUTINE  (FURNISHED 
BY  THE  USER)  AND  SUBROUTINE  RKOS.  EXCEPT  PRMT(S) 
THE  COMPONENTS  ARE  NOT  DESTROYED  BY  SUBROUTINE 
RKOS  AND  THEY  ARE 

PRMTUI-  LOWER  BOUND  OF  THE  INTERVAL  (INPUT). 

PRMTl 2 )-  UPPER  BOUND  OF  THE  INTERVAL  (INPUT). 

PRMT( 3 )-  INITIAL  INCREMENT  OF  THE  INDEPENDENT  VARIABLE 
( INPUT) . 

PRMT( 4 )-  UPPER  ERROR  BOUND  (INPUT).  IF  ABSOLUTE  ERROR  IS 
GREATER  THAN  PRMT(4).  INCREMENT  GETS  HALVEO. 

IF  INCREMENT  IS  LESS  THAN  PRMT(3)  AND  ABSOLUTE 
ERROR  LESS  THAN  PRMT(4)/SO.  INCREMENT  GETS  OOUBLEO 
THE  USER  MAY  CHANGE  PRMT(4)  BY  MEANS  OF  HIS 
OUTPUT  SUBROUTINE. 

PRMT( 5 )-  NO  INPUT  PARAMETER.  SUBROUTINE  RKOS  INITIALIZES 
PRMT ( 5 ) =0 . IF  THE  USER  WANTS  TO  TERMINATE 
SUBROUTINE  RKGS  AT  ANY  OUTPUT  POINT.  HE  HAS  TO 
CHANGE  PRMT ( 5 ) TO  NON-ZERO  BY  MEANS  OF  SUBROUTINE 
OUTP.  FURTHER  COMPONENTS  OF  VECTOR  PRMT  ARE 
FEASIBLE  IF  ITS  DIMENSION  IS  DEFINED  GREATER 
THAN  5.  HOWEVER  SUBROUTINE  RKGS  DOES  NOT  REQUIRE 
AND  CHANGE  THEM.  NEVERTHELESS  THEY  MAY  BE  USEFUL 
FOR  HANDING  RESULT  VALUES  TO  THE  MAIN  PROORAM 
(CALLING  RKGS)  WHICH  ARE  OBTAIN'D  BY  SPECIAL 
MANIPULATIONS  WITH  OUTPUT  DATA  IN  SUBROUTINE  OUTP. 

Y - INPUT  VECTOR  OF  INITIAL  VALUES.  (DESTROYED) 

LATERON  Y IS  THE  RESULTING  VECTOR  OF  DEPENDENT 
VARIABLES  COMPUTED  AT  INTERMEDIATE  POINTS  X. 


RKOS  10 
•RKOS  20 
RKGS  30 
RKOS  40 
RKOS  SO 
RKGS  60 
RKGS  70 
RKOS  80 
RKGS  90 
RKGS  100 
RKGS  110 
RKOS  120 
RKGS  130 
RKOS  140 
RKGS  150 
RKGS  160 
RKGS  170 
RKGS  180 
RKGS  190 
RKGS  200 
RKGS  210 
RKOS  220 
RKGS  230 
RKGS  240 
RKGS  250 
RKGS  260 
RKGS  270 
RKOS  280 
.RKGS  290 
RKGS  300 
RKGS  310 
RKGS  320 
RKGS  330 
RKGS  340 
RKGS  350 
RKOS  360 
RKGS  370 
RKOS  380 
RKGS  390 
RKGS  400 
RKGS  410 
RKGS  420 
RKGS  430 
RKGS  440 
RKGS  450 
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OERY  - INPUT  VECTOR  OT  ERROR  WEIGHTS.  (DESTROYED)  RKOS 

THE  SUM  OF  ITS  COMPONENTS  MUST  BE  EQUAL  TO  1.  RKOP 

LATERON  DERY  IS  THE  VECTOR  OF  DERIVATIVES.  WHICH  RKOS 

BELONG  TO  FUNCTION  VALUES  Y AT  A POINT  X.  RKOS 

NDIM  - AN  INPUT  VALUE,  WHICH  SPECIFIES  THE  NUMBER  OF  RKOS 

EQUATIONS  IN  THE  SYSTEM.  RKOS 

IHLF  - AN  OUTPUT  VALUE.  WHICH  SPECIFIES  THE  NUMBER  OF  RKOS 
BISECTIONS  OF  THE  INITIAL  INCREMENT.  IF  IHLF  OETS  RKOS 
GREATER  THAN  10.  SUBROUTINE  RKOS  RETURNS  WITH  RKOS 

ERROR  MESSAGE  IHLF=11  INTO  MAIN  PROORAM.  ERROR  RKOS 

MESSHOE  IHLF=12  OR  IHLF=13  APPEARS  IN  CASE  RKOS 

PRMT( 3 ) =0  OR  IN  CASE  SI GN( PRMTt 3 ) ) .NE .SION( PRMT( 2 )-RKGS 
PRMTC 1)  ) RESPECTIVELY.  RKOS 

FCT  - THE  NAME  OF  AN  EXTERNAL  SUBROUTINE  USEO.  THIS  RKOS 

SUBROUTINE  COMPUTES  THE  RIOHT  HAND  SIDES  DERY  OF  RKOS 

THE  SYSTEM  TO  GIVEN  VALUES  X AND  Y.  ITS  PARAMETER  RKOS 
LIST  MUST  BE  X.Y.DERY.  SUBROUTINE  FCT  SHOULD  RKOS 

NOT  DESTROY  X AND  Y . RKOS 

OUTP  - THE  NAME  OF  AN  EXTERNAL  OUTPUT  SUBROUTINE  USED.  RKOS 

ITS  PARAMETER  LIST  MUST  BE  X .Y .DERY , IHLF .NDIM .PRMT .RKOS 
NONE  OF  THESE  PARAMETERS  (EXCEPT.  IF  NECESSARY,  RKOS 
PRMTt  4).  PRMT(  S ),...)-  SHOULD  BE  CHANGED  BY  RKOS 

SUBROUTINE  OUTP.  IF  PRMTt S)  IS  CHANGED  TO  NON-ZERO, RKOS 
SUBROUTINE  RKOS  IS  TERMINATED.  RKOS 

AUX  - AN  AUXILIARY  STORAGE  ARRAY  WITH  8 ROWS  AND  NDIM  RKOS 

COLUMNS.  RKOS 

RKOS 

REMARKS  RKOS 

THE  PROCEDURE  TERMINATES  AND  RETURNS  TO  CALLING  PROORAM.  IF  RKOS 

(1)  MORE  THAN  10  BISECTIONS  OF  THE  INITIAL  INCREMENT  ARE  RKOS 

NECESSARY  TO  OET  SATISFACTORY  ACCURACY  (ERROR  MESSAGE  RKGS 

IHLFsll).  RKOS 

(2)  INITIAL  INCREMENT  IS  EQUAL  TO  0 OR  HAS  WRONG  SION  RKOS 

(ERROR  MESSAGES  IHLF=12  OR  iHLF=13),  RKOS 

(3)  THE  WHOLE  INTEGRATION  INTERVRL  IS  WORKED  THROUGH.  RKOS 

(4)  SUBROUTINE  OUTP  HAS  CHANGED  PRMT(S)  TO  NON-ZERO . RKOS 

RKOS 

SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIREO  RKOS 

THE  EXTERNAL  SUBROUTINES  FCT( X . Y .OERY ) AND  RKOS 

OUTPtX.Y, OERY. IHLF. NDIM. PRMT)  MUST  BE  FURNISHED  BY  THE  USER. RKOS 

RKOS 

METHOO  RKOS 

EVALUATION  S DONE  BY  MEANS  OF  FOURTH  ORDER  RUNOE-KUTTA  RKOS 

FORMULAE  IN  THE  MODIFICATION  DUE  TO  OILL.  ACCURACY  IS  RKOS 

TESTED  COMPARING  THE  RESULTS  OF  THE  PROCEDURE  WITH  SINOLE  RKOS 


460 
470 
480 
490 
SOO 
S10 
520 
530 
540 
550 
560 
570 
580 
590 
600 
610 
620 
630 
640 
650 
660 
670 
680 
690 
700 
710 
720 
730 
740 
7 SO 
760 
770 
780 
790 
800 
810 
820 
830 
840 
850 
860 
870 
880 
890 
900 
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AND  DOUBLE  INCREMENT. 

SUBROUTINE  RKOS  AUTOMATICALLY  RDJUSTS  THE  INCREMENT  OURINO 
THE  WHOLE  COMPUTATION  BY  HALVING  OR  DOUBLINO.  IF  MORE  THAN 
10  BISECTIONS  OF  THE  INCREMENT  ARE  NECESSARY  TO  OCT 
SATISFACTORY  ACCURACY.  THE  SUBROUTINE  RETURNS  WITH 
ERROR  MESSAGE  IHLF=11  INTO  MAIN  PROGRAM. 

TO  GET  FULL  FLEXIBILITY  IN  OUTPUT.  AN  OUTPUT  SUBROUTINE 
MUST  BE  FURNISHED  BY  THE  USER. 

FOR  REFERENCE.  SEE 

RALSTON/WILF.  MATHEMATICAL  METHODS  FOR  DIGITAL  COMPUTERS. 
WILEY.  NEW  YORK/LONDON.  1960.  PP. 110-120. 


SUBROUTINE  RKGS( PRMT . Y .OERY .NDIM , I HLF.FCT .OUTP .AUX ) 


DIMENSION  Y( 5 ) .DERY( 5 ) ,AUX( 8 ,5 ) ,A(4),B(4).C(4) ,PRMT( 5 ) 

DO  1 1=1. NDIM 

AUX( 8.1  )=.06666667»0ERY( I ) 

X=PRMT( 1 ) 

XEN0=PRMT( 2 ) 

H=PRMT( 3 ) 

PRMT (5 )=0. 

CALL  FCT(X.Y.OERY) 

ERROR  TEST 

IFl H»(XEND-X ) 135.37.2 

PREPARATIONS  FOR  RUNGE-KUTTA  METHOO 

A( 1 )=.S 

A(2)=. 2928932 

A( 3 )=1 .707107 

A(4)=. 1866667 

Bl 1 1=2. 

B(2)=l . 

B(3)=l. 

BU)=2. 

C( 1)=.S 
C( 2 1= .2928932 
CT  3 1 = 1 .707107 
C(4)=.S 


PREPARATIONS  OF  FIRST  RUNOE-KUTTA  STEP 
DO  3 Irl.NOIM 


RKOS  910 
RKOS  920 
RKOS  930 
RKOS  940 
RKOS  950 
RKOS  960 
RKOS  970 
RKOS  980 
RKOS  990 
RKOSIOOO 
RKOSIOIO 
RKOS1020 
RKOS1030 
RKOS1040 
RKOS1050 
RKGS1060 
RKOS 1070 

RKOS  1090 
RKOS 11 00 
RKGSlllO 
RKGS1120 
RKOS 1130 
RKGS1140 
RKOS 11 50 
RKOS 11 60 
RKGSV170 
RK0S1180 
RKOS 1190 
RKOS 1200 
RKOS1210 
RKGS1220 
RKOS1230 
RKOS1240 
RKGS1250 
RKGS1260 
RKOS 1270 
RKOS 1280 
RK0S1290 
RKOS 1300 
RKOS 1310 
RKOS 1320 
RK0S1330 
RKOS1340 
RKOS1350 
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AUX( 1 .1 )=YC I ) 
AUX(  2 . I )=OERY( I ) 
AUX(3.I)=0. 

3  AUXC6.I )=0. 
IREC=0 
H-H+H 
IHLF=-1 
ISTEP=0 
IEND=0 


START  OF  A RUNOE-KUTTA  STEP 

4 IF( ( X+H-XEND )»H )7 .6 .5 

5 N=XENO-X 

6 IENO=l 

RECORDING  OF  INITIAL  VALUES  OF  THIS  STEP 

7 CALL  OUTPCX.Y.OERY.IREC.NOIM.PRMT) 
IF(PRMT(5) 140.8. 40 

8 ITEST=0 

9 ISTEP=ISTEP*L 


START  OF  INNERHOST  RUNOE-KUTTA  LOOP 
J=1 

10  AJ=h( J 1 
BJ=B( J) 

CJ=C( J ) 

DO  11  I=1.NDIH 
R1=H»0ERY( I ) 

R2=A J«l R1-BJ«AUX( 6.11 ) 

Y ( I )=Y(  I )>R2 
R2=R2+R2»R2 

11  AUX( 6 . I )=AUX ( 6 . I )*R2-CJnRl 
IF( J-4 112.15.15 

12  J=J*1 

IFC  J-3  >13,14.13 

13  X=X*.S»H 

14  CALL  FCT(X.Y.DERY) 

GOTO  10 

END  OF  INNERMOST  RUNOE-KUTTA  LOOP 


TEST  OF  ACCURACY 
IS  I F ( ITESTH6.16.20 


RKGS1360 

RKGS1370 

RKGS1380 

RKGS1390 

RKGS1400 

RKGS1410 

RKGS1420 

RKGS1430 

RK0S1440 

RKGS14S0 

RKGS1460 

RKGS1470 

RKGS1480 

RKGS1490 

RKGS1S00 

RKOSISIO 

RKGC1S20 

RKGS1S30 

RKGS1S40 

RKGS1550 

RKGS1S60 

RK0S1S70 

RKGS1SB0 

RK0S1590 

RKGS1600 

RK0S1610 

RKGS1620 

RK0S1630 

RKGS1640 

RK0S16S0 

RKGS1660 

RK0S1670 

RK0S1680 

RKGS1690 

RK0S1700 

RK0S1710 

RKGS1720 

RKGS1730 

RKGS1740 

RKGS1750 

RK0S1760 

RK0S1770 

RKGS1780 

RKOS1790 

RKGS1800 
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IN  CASE  ITEST=0  THERE  IS  NO  POSSIBILITY  FOR  TESTING  OF  ACCURACY 

RKGS1810 

RKGS1820 

16 

00  17  1=1. NOIM 

RKGSL830 

17 

AUX(4.I)=Ym 

RKGS1840 

ITEST=1 

RKGS1850 

ISTEP=I5TEP+ISTEP-2 

RKGS1860 

18 

IHLF=IHLF+1 

RKGS1870 

X=X-H 

RKGS1880 

H=.5»H 

RKGS1890 

00  19  1=1 .NDIM 

RKDS1900 

Y ( I )=RUX(  l.n 

RKGS1910 

DERY( I )=AUX( 2.1) 

RKGS1920 

19 

AUX(6.I )=AUX13.I ) 

RKGS1930 

GOTO  9 

RKGS1940 

TN  CASE  ITESTrl  TESTING  OF  ACCURACY  fS  POSSIBLE 

RKGS1950 

RKGS1960 

20 

If100=ISTEP/2 

RKGS1970 

IF(  ISTEP-rflOO-IMOO  121.23.21 

RK0S1980 

21 

CALL  FCT(X.Y.DERY) 

RKGS1990 

00  22  1=1 .NOIM 

RKGS2000 

flUXCS. I )=Y( I 1 

RKGS2010 

22 

AUX(7.I)=0ERY(  I ) 

RKGS2020 

GOTO  9 

RKGS2030 

COMPUTATION  OF  TEST  VALUE  OELT 

RKOS2040 

RKGS2050 

23 

0ELT=0 . 

RKGS2060 

00  24  1=1 .NDIM 

RKGS2070 

24 

OELT=DELT+flUX( 8. 1 luflBSt  RUXC 4.1 )-Y( I ) ) 

RKGS2080 

IF(DELT-PRMT( 41)28.28.25 

RKGS2090 

ERROR  IS  TOO  GREAT 

RKGS2100 
RKGS21 10 

25 

IFC IHLF-10 126.36.36 

RKGS2120 

26 

00  27  1=1 .NOIM 

RKGS2130 

27 

flux(4.n=fluxcs.n 

RKGS21 40 

ISTEP=ISTEP+ISTEP-4 

RK0S21S0 

X=X-H 

RKGS2160 

IEN0=O 

RKGS2170 

GOTO  18 

RKGS2180 

RESULT  VALUES  ARE  GOOO 

RKGS2190 

RKGS2200 

28 

CALL  FCT(X.Y.OERY) 

RKGS2210 

00  29  1 = 1 .NOIM 

RKGS2220 

flux(i.n=Y(i) 

RKGS2230 

flUXC  2 • I )=DERY( I ) 

RKGS2240 

AUX(3.n=AUXC6.I) 

RKGS2250 

E-22 


I 


o r>  o 


Y(I)=AUX(S.I) 

RKGS2260 

29 

OERY ( n =AUX ( 7 . n 

RKGS2270 

CALL  OUTPC X-H .Y ,OERY , IHLF , NO  HI  .PRMT ) 

RKGS2280 

i 

IF(PRMT( 51 ) 40 >30. 40 

RKGS2290 

30 

00  31  Irl.NOIM 

RKGS2300 

Y(I)sRUX(l.n 

RKGS2310 

31 

DERY( I )=RUX( 2 . 1 1 

RKGS2320 

IREC=IHLF 

RKGS2330 

IF( IEN0 )32 >32  >39 

RKGS2340 

RKGS2350 

INCREMENT  GETS  OOUBLEO 

RKGS2360 

32 

IHLFrlhLF-l 

RKGS2370 

ISTEP=ISTEP/2 

RKGS2380 

H=H«-H 

RKGS2390 

IFUHLF14.33.33 

RKGS2400 

33 

IM00=ISTEP/2 

RKGS2410 

IF( ISTEP-IM00-IM0D14.34.4 

RKGS2420 

34 

IFCOELT- .02»PRMT(  4) 135.35.4 

RKGS2430 

35 

IHLF=IHLF-1 

RKGS2440 

ISTEPrISTEP/2 

RKOS24SO 

H=H+H 

RKGS2460 

GOTO  4 

RKGS2470 

RKGS2480 

RKGS2490 

RKGS250C 

RKGS2510 

RKOS2S20 

RKOS2530 

RKGS2540 

RKGS2550 

RK082560 

RKOS2570 

RKGS2580 

RKOS2590 


RETURNS  TO  CALLING  PROGRAM 

36  IHCF=11 

CALL  FCT(X.Y,DERY1 
GOTO  39 

37  IHLF=12 
GOTO  39 

38  IHLF=13 

39  CALL  OUTPIX .Y ,DERY , IHLF.NDIM .PRHT ) 

40  RETURN 
END 


SUBROUTINE  STPTIS 


PURPOSE:  CALCULATES  THE  MEAN  AND  STANDARD  DEVIATION  OF  ARRAY  0 


LISTING  OF  PARAMETERS: 

INPUT: 

NITER  s NUMBER  IF  DATA  POINTS  ON  THE  TIME  INTERVAL 
Q = ARRAY  OF  INPUT  VALUES 
OUTPUT : 

QBAR  = MEAN  VALUE  OF  ARRAY  Q 

SOQ  s STANDARD  DEVIATION  OF  VALUES  IN  ARRAY  0 


COMMON  BLOCKS  ANO  THE  PROORAMS  THAT  REQUIRE  THEM: 
NONE 

SUBROUTINES  CALLED: 

NONE 


SUBROUTINE  STAPISt NITER .Q .QBAR .SOQ ) 

DIMENSION  O(NITER) 

QBAR  = 0. 

DO  10  1=1. NITER 
10  QBAR  = QBAR  + Q(  I ) 

QBAR  = (QBAR  - 0 .5-C Q( NITER ) ♦ Q(1)))/(NITER 
SDQ  = 0. 

DO  20  1=1 .NITER 
OQ  = (0(1)  - QBAR )«w2 
20  SOQ  = SOQ  + OQ 

SDQ  = SQRTISOQ/NITER) 

RETURN 

END 
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